CN102650701A - Reservoir wave impedance prediction method based on dipole wave - Google Patents
Reservoir wave impedance prediction method based on dipole wave Download PDFInfo
- Publication number
- CN102650701A CN102650701A CN2012101291927A CN201210129192A CN102650701A CN 102650701 A CN102650701 A CN 102650701A CN 2012101291927 A CN2012101291927 A CN 2012101291927A CN 201210129192 A CN201210129192 A CN 201210129192A CN 102650701 A CN102650701 A CN 102650701A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- seismic
- wave
- impedance
- well
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a reservoir wave impedance prediction method based on dipole wave, which comprises the following steps that: logging data is analyzed to acquire the corresponding parameter data of the seismic response characteristics of a reservoir; well seismic calibration is carried out according to the logging and seismic data so as to extract wavelets and determine the time-depth relationship; the seismic response characteristic wave function of impedance changes is constructed according to the parameter data, the wavelets and the time-depth relationship; the way integral characteristics of the logging data are extracted according to the seismic response characteristic wave function to form a seismic response model database; the seismic data beside a well is analyzed according to the seismic response model database to work out an impedance interpretation spectrum of the seismic data beside the well; a seismic trace is analyzed according to the seismic response model database to work out a data body wave impedance spectrum; and the wave impedance is predicted according to the seismic response model database, the impedance interpretation spectrum and the data body wave impedance spectrum. According to the reservoir wave impedance prediction method based on the dipole wave, the wave impedance of the reservoir is predicted through the seismic characteristics and wave impedance change rules based on the analysis to the logging data, so that the precision in predicting the wave impedance of the reservoir and the resolution to thin reservoirs are improved.
Description
Technical field
The present invention is about the geophysical prospecting for oil field, particularly about the reservoirs exploration technical field, is a kind ofly to carry out reservoir wave impedance forecast method based on dipole wave specifically.
Background technology
The geophysical prospecting for oil technology is to use one of the oil of seeking underground reservoir the most widely, the space distribution of rock gas and method of inner structure rule; Rock acoustic feature through the seismic data record is confirmed the reservoir lithology interphase, oil-gas exploration and development is had the meaning of particular importance.
Usually said seismic inversion process promptly is to utilize the counter rock signature and the reservoir parameter of releasing zone of interest of numerical information of earthquake field acquisition.Seismic inversion is the important technology of seismic reservoir prediction, but problems such as complicated geological condition low owing to its precision of prediction, that face are restricting its application in reservoir exploration exploitation reservoir prediction.And how to effectively utilize data such as multiple data such as earthquake, drilling well, geology, well logging improve the reservoir modeling accuracy and simplify constrained parameters, seem particularly important.Organically combine the model seismic inversion technology of setting up according to well-log information and seismic data; Though can solve the resolution problems of some thin interbeds; But should technology have stronger multi-solution, diverse ways, different parameter, different model all possibly obtain different inversion results.
Current seismic inversion technical application is comparatively extensive, and its type mainly contains multiple modes such as recurrence inversion, well log constrained seismic inversion, and the wavelet transformation of seismic inversion, genetic algorithm and neural network algorithm all are non-linear.The general existence of these seismic reservoir wave impedance inversion techniques has the problem of multi-solution, has seriously reduced the accuracy of reservoir prediction, has increased the difficulty of petroleum prospecting.
Summary of the invention
Instance of the present invention provides a kind of and has carried out reservoir wave impedance forecast method based on dipole wave; On the basis of analyzing at well logging, seismic data; The reservoir top, the end seismic facies that utilize the variation of seismic signature and wave impedance to cause answer the waveform of dipole wave and energy road integration feature constraint interative computation to come the reservoir wave impedance is predicted; Improved the accuracy of reservoir wave impedance prediction, improved resolution thin reservoir.
The objective of the invention is, provide a kind of and carry out reservoir wave impedance forecast method based on dipole wave, described method comprises:
Analyze well-log information to obtain seismic reservoir response characteristic corresponding parameters data;
Carrying out well shake demarcation according to described well-log information concerns when wavelet is also confirmed to extract deeply;
Make up the response features wave function according to described supplemental characteristic, described wavelet in time to concern deeply;
According to described response features wave function described well-log information is carried out the integration feature extraction and form the seismic response model database;
Analyze the impedance explanation spectrum that the other geological data of well obtains the other geological data of well according to described seismic response model database;
Analyze seismic trace according to described seismic response model database and obtain data volume wave impedance spectrum;
Explain that according to described seismic response model database, impedance spectrum and data volume impedance spectrum carry out the wave impedance prediction.
Wherein, described supplemental characteristic comprises sound wave, density and the wave impedance reservoir parameter of reservoir thickness data, reservoir upper and lower medium.
According to described well-log information carry out well shake demarcate with extract wavelet and when confirming dark relation comprise: carry out according to described well-log information that big horizon calibration, frequency division are demarcated, wavelet is extracted in the meticulous demarcation of composite traces; Concern deeply when obtaining according to seismogeology layering, seismic interpretation layer position and waveform character.
Making up the response features wave function according to described supplemental characteristic, described wavelet with timely dark relation comprises: analyze corresponding waveform character and the described supplemental characteristic of dark relation when described; Bireflection coefficient and described wavelet are carried out convolution obtain composite traces; Analyze of the influence of seismic wavelet characteristic to reservoir top, the end and border waveform; Make up the response features wave function according to described waveform character, described composite traces and described influence.
According to described response features wave function described reservoir data being carried out the integration feature extraction forms the seismic response model database and comprise: the attenuation gradient fundamental function of take to push up, the end is different retrains described response features wave function; Response features wave function through after the constraint carries out the waveform correlation analysis to well-log information; Described well-log information is carried out the integration feature extraction; Described well-log information is carried out positive and negative interative computation form the seismic response model database.
The impedance that obtains the other geological data of well according to the other geological data of described seismic response model database analysis well explains that spectrum comprises: according to described seismic response model database well side geological data is analyzed; According to the waveform correlation coefficient of dipole wave the other geological data of described well is carried out the impedance that big impedance priority processing and road integral constraint function priority processing obtain the other geological data of well and explain spectrum; Described impedance is explained that spectrum makes an explanation obtain predicting the outcome of dipole wave.
Obtaining data volume wave impedance spectrum according to described seismic response model database analysis seismic trace comprises: according to described seismic response model database three-dimensional wave impedance data body is analyzed and obtained the reservoir impedance spectrum; Described seismic response model database is obtained data volume wave impedance spectrum to seismic trace analysis.
Preferably, described method comprises that also described data volume wave impedance spectrum is carried out 3 D stereo to be shown.
Advantageous effects of the present invention is; On the basis of analyzing at well logging, seismic data; The reservoir top, the end seismic facies that utilize the wave impedance variation to cause answer the waveform of dipole wave and energy road integration feature constraint interative computation to come the reservoir wave impedance is predicted, have improved the accuracy of reservoir wave impedance prediction, have reduced iterations; Improved resolution, particularly improved resolution heterogeneous reservoirs such as carbonatite solution cavity, algal limestone, volcanics to thin reservoir.
Description of drawings
In order to be illustrated more clearly in instance of the present invention or technical scheme of the prior art; To do to introduce simply to the accompanying drawing of required use in instance or the description of the Prior Art below; Obviously, the accompanying drawing in describing below only is instances more of the present invention, for those of ordinary skills; Under the prerequisite of not paying creative work property, can also obtain other accompanying drawing according to these accompanying drawings.
A kind of process flow diagram that carries out the embodiment one of reservoir wave impedance Forecasting Methodology based on dipole wave that Fig. 1 provides for instance of the present invention;
A kind of process flow diagram that carries out the embodiment two of reservoir wave impedance Forecasting Methodology based on dipole wave that Fig. 2 provides for instance of the present invention;
Fig. 3 is the particular flow sheet of the step S103 among Fig. 1;
Fig. 4 is the particular flow sheet of the step S104 among Fig. 1;
Fig. 5 is the step S105 among Fig. 1, the particular flow sheet of S106;
A kind of process flow diagram that carries out the embodiment three of reservoir wave impedance Forecasting Methodology based on dipole wave that Fig. 6 provides for instance of the present invention;
Fig. 7 is the reservoir section well logging well logging sectional view in LTM carbonate reservoir work area;
Fig. 8 is the impedance variations theoretical model figure in LTM carbonate reservoir work area;
Fig. 9 is the theogram illustraton of model in LTM carbonate reservoir work area;
Figure 10 is that the spectrum synoptic diagram is explained in the reservoir impedance in LTM carbonate reservoir work area;
Figure 11 is the sensitive parameter and the effect analysis figure in LTM carbonate reservoir work area.
Embodiment
To combine the accompanying drawing in the instance of the present invention below, the technical scheme in the embodiment of the invention is carried out clear, intactly description, obviously, described instance only is the present invention's part instance, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills are not making the every other embodiment that is obtained under the creative work prerequisite, all belong to the scope of the present invention's protection.
Fig. 1 a kind ofly carries out the process flow diagram of reservoir wave impedance Forecasting Methodology based on dipole wave for what instance of the present invention provided, can be known that by Fig. 1 this method specifically comprises:
S101: analyze well-log information to obtain seismic reservoir response characteristic corresponding parameters data.Well-log information generally comprises the explanation of well-log information individual well, geological data data, logging data data, drilling data etc.The supplemental characteristic that obtains comprises the reservoir thickness data of " in pairs " dipolar character of the end, top reflection, sound wave, density and the wave impedance reservoir parameter of reservoir upper and lower medium.Analyzing the get parms concrete steps of data of well-log information is: at first explain the data parameters of obtaining accurate reflection seismic reservoir response characteristic from the well-log information individual well; Then on the basis that makes full use of earthquake, well logging and logging data data; Analyze zone of interest heterogeneous reservoir characteristic parameter data and logging trace Changing Pattern in conjunction with drilling well, well-log information, explain and analyze sound wave, density and the wave impedance reservoir parameter information of extracting reservoir thickness data and reservoir upper and lower medium simultaneously.
S102: carry out well shake demarcation according to described well-log information and concern deeply when wavelet is also confirmed to extract.Fig. 2 a kind ofly carries out the process flow diagram of the embodiment two of reservoir wave impedance Forecasting Methodology based on dipole wave for what the embodiment of the invention provided, can be known that by Fig. 2 this step specifically comprises:
S202: carry out big horizon calibration, frequency division demarcation, the meticulous demarcation extraction of composite traces wavelet according to described well-log information; Promptly make full use of known well logging, earthquake information, adopt big horizon calibration, frequency division demarcation, the meticulous calibration technique of composite traces to extract wavelet and well shake coupling A accurately
2The coefficient factor.
When wavelet extraction, at first will carry out the structure of objective function: the variation of seismic wave field characteristic is corresponding to the variation of wave impedance; Between objective interval phase neighboring trace, has similarity preferably.Putting before this, make up iterative inversion minimize objective function J (Z, W):
J(Z，W)＝||S-S
1||+A
1||Z-Z
1||+A
2||W-W
0|| (1)
Each meaning of parameters is in the formula (1): S is expressed as seismologic record, A
1, A
2Be expressed as weighting coefficient; In reality extrapolation refutation process, dipole wave W
0Be that ψ and wave impedance Z (t) ask for respectively:
From (1) formula, can know: when being positioned at the well point, wave impedance is known, and core missions are to find the solution wavelet, so objective function is reduced to:
J(Z，W)＝||S-S
1||+A
2||W-W
0|| (2)
In the following formula, can confirm initial dipole wave W with experimental formula according to the spectral amplitude of well lie
0And A
2The coefficient factor.
S203: concern deeply when obtaining according to seismogeology layering, seismic interpretation layer position and waveform character; Promptly utilize multiple information such as seismogeology layering, seismic interpretation layer position and well logging impedance data composite traces and the contrast of seismologic record waveform character to confirm each other, obtain the best the time concern deeply.
S103: make up the response features wave function in time to concern deeply according to described supplemental characteristic, described wavelet.Fig. 3 is the particular flow sheet of step S103, can be known by Fig. 3, and this step specifically comprises:
S301: analyze corresponding waveform character and the described supplemental characteristic of dark relation when described;
S302: the wavelet of bireflection coefficient and described dipole is carried out convolution obtain composite traces, this step is carried out based on convolution principle algorithm, and the composite traces that convolution obtains can be used to analyze wavelet and the synthetic waveform character of reflection coefficient convolution.
S303: the seismic wavelet characteristic of the wavelet correspondence of extracting among the analytical procedure S102 is to the synthetic superimposed influence of reservoir top, the end and the border waveform of different-thickness;
S304: the response features wave function that makes up dipole according to described waveform character, described composite traces and described influence.In concrete embodiment, the response features wave function of the different impedance variations of a series of reservoirs of structure can abbreviate the basic wavelet function of reservoir as, and reflection coefficient is different.
Make up the basic wavelet function of the energy and the response features wave function of the different impedance variations of reservoir of road integral constraint characteristic;
Through small echo ψ and reflection coefficient constraint convolution model ψ
S(t)＝R(t)*W(t) (3)
S (t) seismologic record wherein, R (t) is a reflection coefficient, W (t) is dipole wave or small echo.
S104: according to described response features wave function described well-log information is carried out the integration feature extraction and form the seismic response model database, Fig. 4 is the particular flow sheet of this step S104, can know that by Fig. 4 step S104 specifically comprises:
S401: the attenuation gradient fundamental function (c, d) of take to push up, the end is different retrains described response features wave function.The locality rings when answering the characteristic wave function to change from small to large with impedance parameter, and the attenuation gradient fundamental function of take simultaneously to push up, the end is different retrains, and just can obtain one group of constraint wave function of representing the ascending variation characteristic of the different wave impedance of reservoir.
Each parameter can be expressed as in the formula (4): ψ is a wavelet function; T is the time, and b is the time difference of top bottom boundary, and a is the scale parameter with the different time of change of frequency; C, d are cosine square fringing attenuation function (being cos (x) * * 2), and x is far and near up and down relevant with reservoir.C (i) or d (j-i)=cos (3.14159265* (i/j) * * 2, i=1, j.From formula (4) to recently seeing, the end fringing attenuation gradient function constraint on the seismic response top, more can embody the reservoir country rock up and down waveform the influence of seismic reservoir response characteristic is improved the resolution of earthquake prediction reservoir.On the basis of wavelet and bireflection coefficient convolution, increase epimere constraint function c and hypomere constraint function d; Window function when having adopted different attenuation gradient methods to decay with different fringings; Help contrasting the predicting reliability analysis; To recently seeing, the difference that predicts the outcome is big from the decay of linear attenuation, cosine, three kinds of different attenuation functions of cosine square, explains that to receive the upper and lower waveform influence in border bigger.When the window parameter can keep the characteristic of seismic waveshape effectively when the cosine square decay decayed with the 10ms fringing, also overcome the superimposed influence of border waveform, thereby improve reservoir wave impedance accuracy of predicting to reservoir section seismic response waveform.
S402: well-log information is carried out waveform constraint correlation analysis through the dipole response features wave function after the constraint;
S403: described well-log information is carried out the integration feature extraction;
S404: described well-log information is carried out positive and negative interative computation form each individual well seismic reservoir response model database.Response features wave function through different impedance variations carries out waveform correlation analysis, the extraction of road integration characteristic attribute and FORWARD AND INVERSE PROBLEMS interative computation to the seismic data that comprises abundant information; Promptly predictably ring and answer model database, thereby form each individual well seismic reservoir response model database.
S105: analyze the impedance explanation spectrum that the other geological data of well obtains the other geological data of well according to described seismic response model database;
S106: analyze seismic trace according to described seismic response model database and obtain data volume wave impedance spectrum;
Fig. 5 is the particular flow sheet of step S105, S106, can be known by Fig. 5, and this step specifically comprises:
S501: according to described seismic response model database the other geological data of well is analyzed, the seismic reservoir response model database that promptly utilizes the place, well point to make up comes that the other geological data of well is carried out the different wave impedance of certain thickness reservoir and changes the analysis of reservoir response characteristic.
S502: according to the waveform correlation coefficient of dipole wave the other geological data of described well is carried out the impedance that big impedance priority processing and road integral constraint function priority processing obtain the other geological data of well and explain spectrum; Promptly take big impedance to be preferably principle; Consider waveform constraint related coefficient (one of standard of explaining as the reservoir impedance) and road integral constraint function (be seismic amplitude with energy response, it is not both the data operand different) the multiparameter principle of priority of dipole seismic wavelet (bireflection coefficient seismic wavelet) simultaneously; Spectrum is explained in the impedance that obtains the other geological data of well, can be described as the reservoir impedance and explains spectrum.
Take into full account the influence that the seismic impedance prediction receives waveform; On the basis of space constraint extrapolation; Need to add more geological information; Realize that the thickness error of explaining on seismic amplitude energy and road integral constraint reservoir impedance spectrum and the individual well is less, can organically combine multiple message reflection reservoir impedance spectrum accuracy; Improve the precision that initial model is set up, can find out that from formula analysis the foundation of primary wave impedance model Z1 seems particularly important.The Z1 function can be expressed as:
Z
1＝f(T_length，t_length，TI，coef，con，h，var，dif) (5)
Wherein to be respectively T_length be the big window length factor to meaning of parameters; T_length is the little window length factor, and TI is an integral function, and coef is a related coefficient; Con is stratum geology schema constrains (integration, toplap, a baselap); H is a reservoir thickness, and var is the impedance variations scope, and djf is an integration; Z
1Be the primary wave impedance model;
S503: described impedance is explained that composing the impedance that obtains dipole wave that makes an explanation predicts the outcome.The reservoir impedance explains that the demonstration numerical value of spectrum is predicting the outcome of dipole wave, also is simultaneously related coefficient and the ratio of road integral constraint energy.Fig. 7 is the reservoir section well logging well logging sectional view in LTM carbonate reservoir work area, reproduces reservoir rock signature up and down, is reservoir in the frame of broken lines among Fig. 7; Fig. 8 is the impedance variations theoretical model figure in LTM carbonate reservoir work area, is the influence of research impedance variations to waveform character, and among Fig. 8, ordinate is the time, and horizontal ordinate is a resistance value, and 50% among the figure is meant the variation of reservoir impedance phase for standard value; Fig. 9 is the theogram illustraton of model in LTM carbonate reservoir work area, and the variation that can find out resistance value has the variation of obvious waveform character, and among Fig. 9, horizontal ordinate is the time, and ordinate is a frequency, and 50% among the figure is meant the variation of reservoir impedance; Figure 10 is the reservoir impedance energy group synoptic diagram in LTM carbonate reservoir work area; Concentration of energy is explained a little less than the reservoir thickness multi-solution; Spectrum is explained in the reservoir impedance in carbonate reservoir work area; Among Figure 10, the left side is that resistance value is 2500 o'clock a seismologic record, and the right side is that resistance value is 12500 o'clock a reservoir wave impedance variation prediction spectrum; Can know the result who expresses the variation of reservoir wave impedance by Figure 11, convergence is comparatively stable relatively, the thickness interpretation good reliability.
S504: according to described seismic response model database three-dimensional wave impedance data body is analyzed and to be obtained the reservoir impedance spectrum;
S505: seismic trace analysis is obtained data volume wave impedance spectrum according to described seismic response model database.According to each individual well seismic reservoir response model database, whole three-dimensional wave impedance data body is analyzed, thereby can be obtained the reservoir impedance spectrum; Promptly obtain whole data volume wave impedance spectrum through each seismic trace being carried out analysis interpretation.The corresponding impedance value sign that adopts numerical value to write is taked transparent display mode, improves geological meaning and characteristic that impedance data shows, improves the three-dimensional effect of explaining of body space prediction between purpose section reservoir.
S107: spectrum and data volume impedance spectrum are explained in described seismic response model database, impedance; Through under prerequisites such as wavelet extraction, initial model structure and impedance interative computation; The prediction 3-d seismic data set carries out integrated interpretation, obtains impedance and predicts the outcome.
The ultimate principle of spectrum interative computation is explained in impedance: outside the well point, in the timing window, wavelet is stable in (1) formula.Fixing wavelet is so that find the solution wave impedance.So objective function is reduced to:
J(Z，W)＝||S-S
1||+A
1||Z-Z
1|| (6)
From above formula analysis, from (1) formula, can find out primary wave impedance model Z
1Foundation seem particularly important.Z
1Function can be expressed as:
Z
1＝f(T_length，t_length，TI，coef，con，h，djf) (7)
Wherein to be respectively T_length be the big window length factor to meaning of parameters, and t_length is the little window length factor, and TI is an integral function; Coef is a related coefficient; Con is stratum geology schema constrains (integration, toplap, a baselap), and h is a reservoir thickness, and djf is an integration; Z
1Be the primary wave impedance model, Z is final inversion result, can find out Z net result and Z from (6) formula
1With S and S
1All there is relation very close.
A kind of process flow diagram that carries out the embodiment three of reservoir wave impedance Forecasting Methodology based on dipole wave that Fig. 6 provides for the embodiment of the invention; Can know that by Fig. 6 the embodiment of the invention provides a kind ofly carries out reservoir wave impedance Forecasting Methodology based on dipole wave and comprise also that except above-mentioned steps described data volume wave impedance spectrum is carried out 3 D stereo to be shown.The data volume Acoustic Impedance Data that obtains is carried out 3 D stereo show, improve the sign of Acoustic Impedance Data speed of displaying and geological meaning, adopt the spatial prediction of interlayer purpose section, the Acoustic Impedance Data body space is made an explanation.
Below in conjunction with concrete embodiment, introduce in detail and provided by the inventionly a kind ofly carry out reservoir wave impedance forecast method based on dipole wave.Be applied to TLM carbonate reservoir work area, the performing step that it is concrete:
1. the finishing analysis of well-log information
Draw through analyzing drilling datas such as this carbonate reservoir work area H9, H7, H6c, H10, H11, H12, H13; This area has many mouthfuls of wells to obtain industry oil stream in the reservoir section; The reservoir type is seam hole type; Its logging trace shows as the reservoir characteristic of praetersonic, and H9 proves the existence in this reservoir section seam hole 1 meter of well section 6689.79-6690.79 rice drilling well emptying.This is participated in prediction and has only used a bite H9 well, and other wells are as the checking well.Obtain that sound wave (wave impedance) parameter attribute of country rock is that subsequent step improves parameter at the bottom of the top of reservoir section.Wave impedance 16946, the bottom that can obtain the reservoir top according to statistics is 17226, the reservoir section is 7546.(unit of its wave impedance is: sound wave m/s* density g/cc)
2. the well shake is demarcated
The H9 key well is carried out meticulous demarcation, utilize well shake associating extraction wavelet, the related coefficient of composite traces reaches 97.9%, obtains dark relation and optimum wavelet when best, is familiar with zone of interest thickness seismic reservoir response characteristic better.
3, make up different reservoir impedance seismic response model bank
Utilize wave impedance value, the country rock wave impedance value up and down of the reservoir section of log analysis to set up the wedge set model; With desirable Ricker wavelet of 25Hz or optimum wavelet and the model reflection coefficient convolution that extracts; Obtain the seismic reservoir response model of different impedance variations; Add cosine square attenuation gradient function constraint on this basis, formed final reflection seismic reservoir response characteristic model, form the integration data storehouse simultaneously; Set up the response features database, for next step actual earthquake data three-dimensional prediction provides standard prediction.Impedance variations-50% ,-25%, 0%, 25%, 50% basic model, compartment thickness is 20ms.
4, the reservoir impedance spectrum is set up
At first utilize response features that seismic trace near well is carried out the prediction of reservoir impedance magnitude at the individual well place; Make seismic trace and composite traces that comparatively desirable energy matching relationship arranged through self-adaptation adjustment parameter; Can find out that through the parameter tentative calculation analysis of place, well point the individual well thickness prediction is better with aboveground corresponding thickness.Time is the reservoir thickness of 2ms.Figure 11 is the sensitive parameter and the effect analysis figure in LTM carbonate reservoir work area; Can know by Figure 11; Related coefficient is comparatively responsive to reservoir thickness prediction; Related coefficient and the prediction of seismic amplitude attribute are comparatively responsive to thickness and impedance variations, improve reservoir impedance spectrum prediction effect through these parameter analyses, and the parameter of reliable constraint is provided for next step whole three-dimensional data prediction.
5,3-d seismic data set reservoir wave impedance prediction
On the basis that meticulous well shake is demarcated, objective interval is carried out 3-D seismics layer position T
03Follow the trail of to explain that the locking objective interval helps saving the machine arithmetic time, the method requires not too high to layer position, have only layer position near zone of interest just, window moves up and down the time through layer position, same measurablely goes out reservoir thickness and distributes.
Explain that from the reservoir impedance spectrum predicts the outcome, can better predict the different reservoir resistance values that change, effect is comparatively desirable.
In sum; Useful achievement of the present invention is: provide a kind of and carried out reservoir wave impedance forecast method based on dipole wave; On the basis of analyzing at well logging, seismic data; The reservoir top, the end seismic facies that utilize the wave impedance variation to cause answer the waveform of dipole wave and energy road integration feature constraint interative computation to come the reservoir wave impedance is predicted, have improved the accuracy of reservoir wave impedance prediction, have improved the resolution to thin reservoir.
Advantage of the present invention is:
1. creatively cross well logging, seismic data analysis and obtain model key parameter data such as the physical characterization data such as sound wave, density and wave impedance that information makes up reservoir section and reservoir top, end characteristic; Earthquake information is demarcated (being that the well shake is demarcated), the well shake is united extraction and optimized technology such as wavelet under the constraint of fringing attenuation function, to adopt theogram to make, log well subsequently; Utilize the reservoir thickness of known well logging interpretation to make up the reservoir model of different-thickness and resistance value, make up the identification that the seismic reservoir response characteristic model bank of dipole wave of the bireflection coefficient of different thick and resistance value degree improves reservoir wave impedance modeling and forecasting precision through forward modeling.This method has improved reservoir wave impedance accuracy of predicting, reduces iterations, improves the resolution to thin reservoir.
2. creatively adopt dipole wave attenuation gradient fringing function constraint; Set up and to reflect seismic reservoir response wave shape characteristic model database; The precision that raising is set up the heterogeneous reservoir initial model; Under 90 degree phase transition road integrations and the constraint of seismic waveshape energy feature, realize the prediction of 3D seismic data reservoir wave impedance; It is different from conventional method and is that with well logging information constraint raising modeling precision the quantification that iterative manner can directly utilize the information realization reservoir wave impedance of seismic amplitude, energy, road integration to change is effectively predicted.
Used specific embodiment among the present invention principle of the present invention and embodiment are set forth, the explanation of above embodiment just is used for helping to understand method of the present invention and core concept thereof; Simultaneously, for one of ordinary skill in the art, according to thought of the present invention, the part that on embodiment and range of application, all can change, in sum, this description should not be construed as limitation of the present invention.
Claims (8)
1. one kind is carried out reservoir wave impedance forecast method based on dipole wave, it is characterized in that described method comprises:
Analyze well-log information to obtain seismic reservoir response characteristic corresponding parameters data;
Carrying out well shake demarcation according to described well-log information concerns when wavelet is also confirmed to extract deeply;
Make up the response features wave function according to described supplemental characteristic, described wavelet in time to concern deeply;
According to described response features wave function described well-log information is carried out the integration feature extraction and form the seismic response model database;
Analyze the impedance explanation spectrum that the other geological data of well obtains the other geological data of well according to described seismic response model database;
Analyze seismic trace according to described seismic response model database and obtain data volume wave impedance spectrum;
Explain that according to described seismic response model database, impedance spectrum and data volume impedance spectrum carry out the wave impedance prediction.
2. method according to claim 1 is characterized in that, described supplemental characteristic comprises sound wave, density and the wave impedance reservoir parameter of reservoir thickness data, reservoir upper and lower medium.
3. method according to claim 1 is characterized in that, according to described well-log information carry out well shake demarcate with extract wavelet and when confirming dark relation comprise:
Carry out big horizon calibration, frequency division demarcation, the meticulous demarcation extraction of composite traces wavelet according to described well-log information;
Concern deeply when obtaining according to seismogeology layering, seismic interpretation layer position and waveform character.
4. method according to claim 1 is characterized in that, concerns that with in time dark making up the response features wave function comprises according to described supplemental characteristic, described wavelet:
Analyze corresponding waveform character and the described supplemental characteristic of dark relation when described;
Bireflection coefficient and described wavelet are carried out reservoir dipole wave convolution obtain composite traces;
Analyze the influence synthetic of seismic wavelet characteristic to reservoir top, the end and border waveform;
Make up the response features wave function according to described waveform character, described composite traces and described influence.
5. method according to claim 1 is characterized in that, according to described response features wave function described reservoir data is carried out integration feature extraction formation seismic response model database and comprises:
Take top, bottom boundary differential declines gradient fundamental function that described response features wave function is retrained;
Response features wave function through after the constraint carries out the waveform correlation analysis to well-log information;
Described well-log information is carried out the integration feature extraction to be analyzed;
Described well-log information is carried out positive and negative interative computation form each individual well seismic reservoir response model database.
6. method according to claim 1 is characterized in that, the impedance that obtains the other geological data of well according to the other geological data of described seismic response model database analysis well explains that spectrum comprises:
According to described seismic response model database the other geological data of well is analyzed;
According to the waveform correlation coefficient of dipole wave the other geological data of described well is carried out the impedance that big impedance priority processing and road integral constraint function priority processing obtain the other geological data of well and explain spectrum;
Described impedance is explained that spectrum makes an explanation obtain predicting the outcome of dipole wave.
7. method according to claim 1 is characterized in that, obtains data volume wave impedance spectrum according to described seismic response model database analysis seismic trace and comprises:
According to described seismic response model database three-dimensional wave impedance data body is analyzed and to be obtained the reservoir impedance spectrum;
According to described seismic response model database seismic trace analysis is obtained data volume wave impedance spectrum.
8. method according to claim 7 is characterized in that, described method comprises that also described data volume wave impedance spectrum is carried out 3 D stereo to be shown.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210129192.7A CN102650701B (en) | 2012-04-27 | 2012-04-27 | Reservoir wave impedance prediction method based on dipole wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210129192.7A CN102650701B (en) | 2012-04-27 | 2012-04-27 | Reservoir wave impedance prediction method based on dipole wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102650701A true CN102650701A (en) | 2012-08-29 |
| CN102650701B CN102650701B (en) | 2014-02-05 |
Family
ID=46692741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210129192.7A Active CN102650701B (en) | 2012-04-27 | 2012-04-27 | Reservoir wave impedance prediction method based on dipole wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102650701B (en) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967884A (en) * | 2012-11-13 | 2013-03-13 | 中国石油大学(北京) | Method and device for evaluating reliability of wave impedance inversion data |
| CN103149589A (en) * | 2013-02-22 | 2013-06-12 | 中国石油天然气股份有限公司 | Igneous rock oil and gas exploration method and device |
| CN103439740A (en) * | 2013-07-24 | 2013-12-11 | 中国石油天然气股份有限公司 | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integrals |
| CN104297791A (en) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | Inversion method and system based on earthquake dominant frequency |
| CN105093312A (en) * | 2015-07-06 | 2015-11-25 | 中国石油天然气股份有限公司 | Seismic relative wave impedance predicting method and apparatus based on frequency domain multi-level differentiations |
| CN106842316A (en) * | 2016-12-20 | 2017-06-13 | 中国石油天然气股份有限公司 | Crack determines method and apparatus |
| CN106842318A (en) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | Microcosmic erosion hole geophysics two-dimensional characterization determines method and apparatus |
| CN108594304A (en) * | 2018-07-25 | 2018-09-28 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | Based on the multipole Cooley impedance inversion approach of linear programming for solution L1 norms |
| CN110941271A (en) * | 2019-11-26 | 2020-03-31 | 北方工业大学 | Automatic driving lane changing model, automatic driving lane changing method and system thereof |
| CN111158048A (en) * | 2020-01-04 | 2020-05-15 | 杨林海 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
| CN111427086A (en) * | 2020-04-02 | 2020-07-17 | 同济大学 | Method for predicting maximum attenuation parameter of stratum by utilizing amplitude frequency component of longitudinal wave and transverse wave |
| CN111751877A (en) * | 2019-03-26 | 2020-10-09 | 中国石油天然气股份有限公司 | Seismic data multiple-integration coherent fracture prediction method and device |
| CN111897004A (en) * | 2020-08-31 | 2020-11-06 | 北京瑞码恒杰科技有限公司 | Logging prediction method based on big data analysis technology |
| CN112289385A (en) * | 2020-09-17 | 2021-01-29 | 西南交通大学 | Electrochemical impedance spectrum prediction method for high-power proton exchange membrane fuel cell stack |
| CN114994758A (en) * | 2022-08-02 | 2022-09-02 | 北京京鲁聚源能源科技有限公司 | Wave impedance extraction and structure characterization method and system for carbonate fracture control reservoir |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107422384B (en) * | 2017-03-20 | 2019-03-01 | 中石化石油工程技术服务有限公司 | A method of vertical seismic profile data is imitated using well-log information |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
| CN1501097A (en) * | 1999-01-15 | 2004-06-02 | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method | |
| CN101158724A (en) * | 2007-09-14 | 2008-04-09 | 中国石油集团西北地质研究所 | Reservoir thickness prediction method based on dipolar wavelet |
| CN102736107A (en) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | Energy constraint heterogeneous reservoir thickness identification system |
-
2012
- 2012-04-27 CN CN201210129192.7A patent/CN102650701B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
| CN1501097A (en) * | 1999-01-15 | 2004-06-02 | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method | |
| CN101158724A (en) * | 2007-09-14 | 2008-04-09 | 中国石油集团西北地质研究所 | Reservoir thickness prediction method based on dipolar wavelet |
| CN102736107A (en) * | 2011-04-07 | 2012-10-17 | 中国石油天然气股份有限公司 | Energy constraint heterogeneous reservoir thickness identification system |
Non-Patent Citations (5)
| Title |
|---|
| 刘彦君等: "沉积规律控制下的测井约束波阻抗反演及其应用", 《大庆石油地质与开发》 * |
| 段云卿等: "储层地震反演在辽河油田大民屯凹陷的应用", 《地球科学(中国地质大学学报)》 * |
| 王西文: "提取地震信号高频信息方法的若干问题", 《石油地球物理勘探》 * |
| 王西文: "相对波阻抗数据体约束下的多井测井参数反演方法及应用", 《石油地球物理勘探》 * |
| 薄永德等: "苏北盆地高邮凹陷隐蔽油气藏勘探方法", 《天然气工业》 * |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967884B (en) * | 2012-11-13 | 2016-01-20 | 中国石油大学(北京) | Reliability of wave impedance inversion data evaluation method and device |
| CN102967884A (en) * | 2012-11-13 | 2013-03-13 | 中国石油大学(北京) | Method and device for evaluating reliability of wave impedance inversion data |
| CN103149589A (en) * | 2013-02-22 | 2013-06-12 | 中国石油天然气股份有限公司 | Igneous rock oil and gas exploration method and device |
| CN103149589B (en) * | 2013-02-22 | 2016-01-06 | 中国石油天然气股份有限公司 | Pyrogenic rock Petroleum Exploration Methods and device |
| CN103439740A (en) * | 2013-07-24 | 2013-12-11 | 中国石油天然气股份有限公司 | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integrals |
| CN104297791A (en) * | 2014-09-25 | 2015-01-21 | 中国石油天然气股份有限公司 | Inversion method and system based on earthquake dominant frequency |
| CN105093312A (en) * | 2015-07-06 | 2015-11-25 | 中国石油天然气股份有限公司 | Seismic relative wave impedance predicting method and apparatus based on frequency domain multi-level differentiations |
| CN105093312B (en) * | 2015-07-06 | 2017-06-13 | 中国石油天然气股份有限公司 | Earthquake Relative Wave Impedance Forecasting Methodology and device based on frequency domain pluridifferentiation |
| CN106842316A (en) * | 2016-12-20 | 2017-06-13 | 中国石油天然气股份有限公司 | Crack determines method and apparatus |
| CN106842318B (en) * | 2016-12-30 | 2019-01-18 | 中国石油天然气股份有限公司 | Microcosmic erosion hole geophysics two-dimensional characterization determines method and apparatus |
| CN106842318A (en) * | 2016-12-30 | 2017-06-13 | 中国石油天然气股份有限公司 | Microcosmic erosion hole geophysics two-dimensional characterization determines method and apparatus |
| CN108594304A (en) * | 2018-07-25 | 2018-09-28 | 中国石油化工股份有限公司胜利油田分公司勘探开发研究院 | Based on the multipole Cooley impedance inversion approach of linear programming for solution L1 norms |
| CN111751877A (en) * | 2019-03-26 | 2020-10-09 | 中国石油天然气股份有限公司 | Seismic data multiple-integration coherent fracture prediction method and device |
| CN110941271A (en) * | 2019-11-26 | 2020-03-31 | 北方工业大学 | Automatic driving lane changing model, automatic driving lane changing method and system thereof |
| CN110941271B (en) * | 2019-11-26 | 2023-04-18 | 北方工业大学 | Automatic driving lane changing model, automatic driving lane changing method and system thereof |
| CN111158048A (en) * | 2020-01-04 | 2020-05-15 | 杨林海 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
| CN111158048B (en) * | 2020-01-04 | 2022-10-04 | 杨林海 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
| CN111427086A (en) * | 2020-04-02 | 2020-07-17 | 同济大学 | Method for predicting maximum attenuation parameter of stratum by utilizing amplitude frequency component of longitudinal wave and transverse wave |
| CN111897004A (en) * | 2020-08-31 | 2020-11-06 | 北京瑞码恒杰科技有限公司 | Logging prediction method based on big data analysis technology |
| CN112289385A (en) * | 2020-09-17 | 2021-01-29 | 西南交通大学 | Electrochemical impedance spectrum prediction method for high-power proton exchange membrane fuel cell stack |
| CN112289385B (en) * | 2020-09-17 | 2022-08-09 | 西南交通大学 | Electrochemical impedance spectrum prediction method for high-power proton exchange membrane fuel cell stack |
| CN114994758A (en) * | 2022-08-02 | 2022-09-02 | 北京京鲁聚源能源科技有限公司 | Wave impedance extraction and structure characterization method and system for carbonate fracture control reservoir |
| CN114994758B (en) * | 2022-08-02 | 2022-10-28 | 北京京鲁聚源能源科技有限公司 | Wave impedance extraction and structure characterization method and system for carbonate rock fracture control reservoir |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102650701B (en) | 2014-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102650701B (en) | Reservoir wave impedance prediction method based on dipole wave | |
| CN101158724B (en) | Reservoir thickness prediction method based on dipolar wavelet | |
| CN102736107B (en) | Energy constraint heterogeneous reservoir thickness identification system | |
| CN102759745B (en) | Method for forecasting carbonate reservoir based on forward modeling of digital geological outcrop model | |
| CN102053270B (en) | Sedimentary formation unit-based seismic facies analysis method | |
| CN104516018B (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| CN109143361B (en) | Method for compiling carbonate rock stratum paleogeology map based on sequence stratigraphy | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| CN109541685B (en) | River channel sand body identification method | |
| CN102466815A (en) | Triassic clastic rock petroleum reservoir identification method | |
| Aminzadeh et al. | Geophysics for petroleum engineers | |
| CN106501854A (en) | Oil gas drilling target geological quantifying risk evaluation methodology | |
| CN101872024A (en) | Method for carrying out well design by using time-lapse seismic | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| CN106353806B (en) | It is a kind of to detect the connective method of reef reservoir | |
| US20210396897A1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
| Roy et al. | Integrated characterization of hydraulic fracture treatments in the Barnett Shale: The Stocker geophysical experiment | |
| CN106526697A (en) | Method for identifying main transform plane of delta construction and destruction | |
| Venieri et al. | Vertical and lateral facies variability in organic-rich mudstones at the reservoir scale: A case study from the Devonian Duvernay formation of Alberta, Canada | |
| Khan et al. | Three-dimensional Structural Modeling (3D SM) and Joint Geophysical Characterization (JGC) of Hydrocarbon Reservoir: A Case Study of the Kadanwari field in Middle Indus Basin (MIB), Southeastern Pakistan | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
| Naseer | Application of instantaneous spectral decomposition-based porosity simulations for imaging shallow-marine stratigraphic traps of Lower-Eocene carbonates sequences of Indus Basin, Onshore Pakistan | |
| CN109283577A (en) | A kind of seismic layer labeling method | |
| CN106990433A (en) | A kind of recognition methods of the small erosion channel in massif |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |