CN103176211B - Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter - Google Patents
Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter Download PDFInfo
- Publication number
- CN103176211B CN103176211B CN201310072499.2A CN201310072499A CN103176211B CN 103176211 B CN103176211 B CN 103176211B CN 201310072499 A CN201310072499 A CN 201310072499A CN 103176211 B CN103176211 B CN 103176211B
- Authority
- CN
- China
- Prior art keywords
- gas
- parameter
- intersection
- data
- wave
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Embodiments provide a kind of gas-bearing reservoir prediction method based on many sensibility elasticities parameter and device, described method comprises: according to the geophysical character of Effective Reservoirs acoustic wave train logging data, determines the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir; Utilize Prestack seismic data to carry out prestack Simultaneous Inversion, obtain elastic parameter section; Project on described seismic elastic parameter inverting section obtain sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir; According to described sensibility elasticity parameter intersection section delineation gas-bearing reservoir scope, and pick up gas sand thickness.Method of the present invention improves the ability of application earthquake prestack inversion prediction Effective Reservoirs, can reduce the risk of explanation, improve the precision of reservoir prediction.
Description
Technical field
The invention belongs to geophysical prospecting for oil interpretation technique, is that one is applicable to based on acoustic wave train logging data and Prestack seismic data use in conjunction in seismic prospecting, the method for effective (gassiness) reservoir thickness of quantitative interpretation that crossed by sensibility elasticity parameter.
Background technology
Ago seismic exploration is with exploitation, the prediction of reservoir thickness adopts post-stack inversion technology often, it take well-log information as constraint, with seismic interpretation layer position for controlling, from well, set up impedance initial value model by interpolation extrapolation, by the error analysis of composite traces and actual seismic data, adopt method of conjugate gradient constantly to revise impedance initial value model, until composite traces and real seismic record the best are approached, surge impedance model is now inversion result.The time thickness of sand body is explained in sand body wave impedance scope pickup according to well logging statistics on inverting section, then is multiplied by sandstone speed, the thickness of sandstone, to reach the object of predicting reservoir thickness.After Soviet Union's Sulige gas field Sandstone Gas Bearing, sandstone velocity of longitudinal wave obviously reduces, cause the wave impedance of gas sand and mud stone or Sandy Silt impedance close or stacked, adopt post-stack inversion in the past to obtain single p-wave impedance and there is multi-solution, thus can not Accurate Prediction reservoir thickness.Therefore, this just needs a kind of new method to predict effectively (gassiness) reservoir thickness.
In recent years, prestack inversion grows up just and solves a kind of stacked reservoir thickness prediction method of p-wave impedance, it is method qualitative forecasting reservoir on inverting section by determining certain single sensibility elasticity parameter threshold value that current prestack inversion method is studied mostly, precision is lower, there is instability and multi-solution.
Summary of the invention
The present invention seeks to, by to Effective Reservoirs acoustic wave train logging geophysical character and elastic parameter analysis, the best sensibility elasticity parameter of preferred gas-bearing reservoir, obtains elastic parameter section by carrying out prestack inversion to seismic data, finally carries out elasticity intersection interpretation prediction gas-bearing reservoir thickness.
For reaching above-mentioned purpose, embodiments provide a kind of gas-bearing reservoir prediction method based on many sensibility elasticities parameter, described method comprises:
According to the geophysical character of Effective Reservoirs acoustic wave train logging data, determine the sensibility elasticity parameter can distinguishing lithology and gas-bearing reservoir;
Utilize Prestack seismic data to carry out prestack Simultaneous Inversion, obtain seismic elastic parameter inverting section;
Project on described seismic elastic parameter inverting section obtain sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir;
According to described sensibility elasticity parameter intersection section delineation gas-bearing reservoir scope, and pick up gas sand thickness.
On the other hand, the embodiment of the present invention additionally provides a kind of gas-bearing reservoir prediction device based on many sensibility elasticities parameter, and described device comprises:
Well logging sensibility elasticity parameter determination unit, for the geophysical character according to Effective Reservoirs acoustic wave train logging data, determines the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir;
Seismic elastic parameter inverting section acquiring unit, for utilizing Prestack seismic data to carry out prestack joint inversion, obtains seismic elastic parameter inverting section;
Sensibility elasticity parameter intersection section acquiring unit, for being projected on described seismic elastic parameter inverting section obtained sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir;
Gas-bearing reservoir delineation unit, for drawing a circle to approve gas-bearing reservoir scope according to described sensibility elasticity parameter intersection section, and picks up gas sand thickness.
The Advantageous Effects of technique scheme of the present invention is:
The gas-bearing reservoir prediction method crossed based on sensibility elasticity parameter of the embodiment of the present invention and device, employ multiple elastic parameter intersection, more stable than unitary elasticity parameter prediction reservoir, decrease multi-solution, for reservoir thickness sxemiquantitative and quantitative interpretation provide good instrument.Well logging Rock physical analysis higher for longitudinal frame and the higher earthquake prestack inversion of lateral resolution combine by the intersection of sensibility elasticity parameter, improve the ability of application earthquake prestack inversion prediction Effective Reservoirs, the risk of explanation can be reduced, improve the precision of reservoir prediction.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, introduce doing one to the accompanying drawing used required in embodiment or description of the prior art simply below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the basic flow sheet of the method for the embodiment of the present invention;
Fig. 2 is that the S-wave impedance SI of the embodiment of the present invention and P-S wave velocity ratio Vp/Vs cross figure;
Fig. 3 is the prestack Simultaneous Inversion sectional view of the embodiment of the present invention;
Fig. 4 is the sensibility elasticity parameter intersection sectional view of the embodiment of the present invention;
Fig. 5 is the Effective Reservoirs figure of the embodiment of the present invention;
Fig. 6 is the allomeric function block diagram of the device of the embodiment of the present invention;
Fig. 7 is the concrete function block diagram of sensibility elasticity parameter determination unit in embodiment of the present invention Fig. 6;
Fig. 8 is the concrete function block diagram of seismic elastic parameter inverting section acquiring unit in embodiment of the present invention Fig. 6.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
For solving Problems existing in background technology, the gas-bearing reservoir prediction method based on many sensibility elasticities parameter of the embodiment of the present invention is from the geophysical character of gas-bearing reservoir acoustic wave train logging data, preferably can distinguish the sensibility elasticity parameter of reservoir and gas-bearing reservoir, prestack Simultaneous Inversion is carried out with Prestack seismic data, obtain seismic elastic parameter inverting section, elastic parameter intersection section is obtained again by the intersection of gas-bearing reservoir sensibility elasticity parameter, finally according to intersection section delineation gas-bearing reservoir scope, pickup gas sand thickness.Concrete steps are as follows:
Refer to the overall flow figure of a kind of gas-bearing reservoir prediction method based on many sensibility elasticities parameter of the embodiment of the present invention shown in Fig. 1:
(1) preferred sensibility elasticity parameter, i.e. step 110: according to the geophysical character of Effective Reservoirs acoustic wave train logging data, determine the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir.
Study area and 20 mouthfuls of full wave train log objective interval reservoir properties are around added up, (formula 1.-Shi 4.) calculates the elastic curve (as p-wave impedance AI, S-wave impedance SI, P-S wave velocity ratio γ, Poisson ratioσ etc.) of acoustic wave train logging data respectively according to the following formula, and crossplot analysis and its is carried out to these elastic curve data, result shows: P-S wave velocity ratio γ, Poisson ratioσ convection cell are comparatively responsive, utilize the intersection of S-wave impedance SI and P-S wave velocity ratio γ to explain and can distinguish gas-bearing horizon (Effective Reservoirs) preferably.Utilize the intersection of elastic curve S-wave impedance SI and P-S wave velocity ratio γ to set up the X plot reflecting gas-bearing formation, determine the codomain scope of gas-bearing reservoir.
AI＝Vp*ρ①
SI＝Vs*ρ②
γ＝Vp/Vs③
σ＝(Vp
2-2Vs
2)/2(Vp
2-Vs
2)④
In formula: Vp is velocity of longitudinal wave, Vs is shear wave velocity, and ρ is density, and AI is p-wave impedance, SI is S-wave impedance, and γ is P-S wave velocity ratio, Poisson ratioσ, represent the scale-up factor of object transverse strain and longitudinal strain, also known as Poisson ratio, instruction change of fluid.
(2) Earthquake Resilient data volume is asked in prestack inversion, i.e. step 120: utilize Prestack seismic data to carry out prestack Simultaneous Inversion, obtains seismic elastic parameter inverting section.
1. according to zone of interest buried depth and geophone offset seismic data state of signal-to-noise far away, the geological data choosing maximum offset scope does common offset superposition, obtain at least three partial stack seismic data volumes (nearly geophone offset superposition, the superposition of middle geophone offset, geophone offset superposition far away), partial stack data volume carries out relative amplitude preserved processing, and require that its frequency, phase equalization are good, keep the prestack dynamic characteristic of seismic trace.
2. near to earthquake, in, geophone offset partial stack data volume far away utilizes acoustic wave train logging data to carry out repeatedly well shake respectively and demarcates, each partial stack data volume is extracted to the wavelet of corresponding angular range, wavelet lengths is between 100-200ms, wavelet half length can be divided exactly by earthquake sampling rate, during wavelet, window is 2-5 wavelet lengths doubly, near zone of interest, avoid tomography, good seismic data is selected to extract, different wavelet ensures that within the scope of dominant frequency sub-wave phase is more stable, and the simple shape of extraction wavelet, concentration of energy, secondary lobe are few.
3. the earthquake inputting the superposition of three different angles scopes is near, in, the wavelet of geophone offset superposition of data body far away and correspondence, solve full Knott-zoeppritz equation (formula 5.) with acoustic wave train logging data constraint employing and carry out prestack Simultaneous Inversion, obtain p-wave impedance data volume AI, S-wave impedance data volume SI and the elastic data body such as density data body ρ and P-S wave velocity ratio Vp/Vs.
Full Knott-zoeppritz equation:
⑤
In formula: r
ppfor longitudinal wave reflection coefficient; r
psfor converted shear wave reflection coefficient; r
spfor compressional wave transmission coefficient; r
ssfor converted shear wave transmission coefficient; ρ
1and ρ
2be respectively the density of upper and lower rock stratum, interface; v
p1and v
p2be respectively the velocity of longitudinal wave of upper and lower rock stratum, interface; v
s1and v
s2be respectively the shear wave velocity of upper and lower rock stratum, interface; θ
1for incident compressional angle and reflection angle; θ
2for the refraction angle through compressional wave;
the reflection angle of converted shear wave;
through the refraction angle of converted shear wave.
(3) intersection of sensibility elasticity parameter, i.e. step 130: project on described seismic elastic parameter inverting section obtain sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir.
Utilize intersection function, by the intersection section that gas-bearing reservoir sensibility elasticity parameter S-wave impedance and P-S wave velocity ratio intersection obtain, the region that intersection identifies meets petrophysical parameter feature, finally obtains elastic parameter intersection section, thus realizes the lateral prediction of Effective Reservoirs.
(4) Effective Reservoirs thickness is calculated, i.e. step 140: according to described sensibility elasticity parameter intersection section delineation gas-bearing reservoir scope, and pick up gas sand thickness.
On elastic parameter intersection section, draw a circle to approve gas-bearing reservoir scope according to intersection region, pickup gas sand time thick T, (formula 6.) calculates Effective Reservoirs thickness H according to the following formula finally, and V is gassiness (effectively) reservoir velocities.
H＝V*T/2⑥
By to intersection Interpretation of profile can quantitatively, sxemiquantitative prediction Effective Reservoirs thickness.
The gas-bearing reservoir prediction method crossed based on sensibility elasticity parameter of the embodiment of the present invention employs multiple elastic parameter intersection, more stable than unitary elasticity parameter prediction reservoir, decrease multi-solution, for reservoir thickness sxemiquantitative and quantitative interpretation provide good instrument.Well logging Rock physical analysis higher for longitudinal frame and the higher earthquake prestack inversion of lateral resolution combine by the intersection of sensibility elasticity parameter, improve the ability of application earthquake prestack inversion prediction Effective Reservoirs, the risk of explanation can be reduced, improve the precision of reservoir prediction.
Embodiment 1:
Utilize the intersection of acoustic wave train logging data S-wave impedance SI and P-S wave velocity ratio γ to explain and distinguish gassiness (effectively) reservoir, make gassiness (effectively) reservoir X plot (Fig. 2) with this.
According to zone of interest buried depth and seismic data signal to noise ratio (S/N ratio) feature, angular stack altogether made by the geological data choosing three angular ranges, obtain earthquake near, in, geophone offset partial stack data volume far away, to closely, in, geophone offset partial stack data volume far away utilizes acoustic wave train logging data to carry out repeatedly well shake respectively and demarcates, each partial stack data volume is extracted to the wavelet of corresponding angular range, the earthquake inputting three different angles superpositions is near, in, superposition of data body and corresponding wavelet a long way, solve full Knott-zoeppritz equation and carry out prestack inversion, obtain p-wave impedance data volume AI, the elastic data bodies (Fig. 3) such as S-wave impedance data volume SI and P-S wave velocity ratio γ.
Gassiness (effectively) the reservoir masterplate of acoustic wave train logging data make is mapped on the elastic parameter data volume such as p-wave impedance, S-wave impedance, P-S wave velocity ratio that prestack inversion obtains, obtain elastic parameter to cross explanation of seismic section (Fig. 4), the region that intersection identifies meets gas-bearing reservoir well logging intersection region petrophysical parameter feature.
Make an explanation to elastic parameter intersection section intersection region, obtain Effective Reservoirs time thickness profile figure (Fig. 5), pickup gas sand time thickness, finally calculates Effective Reservoirs thickness.
The advantage of the inventive method embodiment is: it is preferred to the present invention is based on sensibility elasticity parameter, Prestack seismic data is carried out prestack Simultaneous Inversion, and made an explanation by the intersection of sensibility elasticity parameter and can predict gassiness (effectively) reservoir thickness quantitatively, the method can be predicted Effective Reservoirs thickness easily and fast, reduces the multi-solution of the parameter of prestack inversion unitary elasticity in the past qualitative forecasting reservoir.
Further, the embodiment of the present invention additionally provides a kind of gas-bearing reservoir prediction device based on many sensibility elasticities parameter, and as shown in Figure 6, described device 600 comprises:
Well logging sensibility elasticity parameter determination unit 610, for the geophysical character according to Effective Reservoirs acoustic wave train logging data, determines the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir;
Seismic elastic parameter inverting section acquiring unit 620, for utilizing Prestack seismic data to carry out prestack Simultaneous Inversion, obtains seismic elastic parameter inverting section;
Sensibility elasticity parameter intersection section acquiring unit 630, for being projected on described seismic elastic parameter inverting section obtained sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir;
Gas-bearing reservoir delineation unit 640, for drawing a circle to approve gas-bearing reservoir scope according to described sensibility elasticity parameter intersection section, and picks up gas sand thickness.
Preferably, consult Fig. 7, described well logging sensibility elasticity parameter determination unit 610, can comprise:
Elastic parameter computing module 6101, for adding up study area and many mouthfuls of full wave train log objective interval reservoir properties around, calculate the elastic parameter of acoustic wave train logging data respectively to the 4th relational expression according to the first relational expression, described elastic parameter comprises: p-wave impedance AI, S-wave impedance SI, P-S wave velocity ratio γ and Poisson ratioσ;
Sensibility elasticity parameter determination module 6102, for carrying out crossplot analysis and its to described elastic parameter, determines that P-S wave velocity ratio γ, Poisson ratioσ are the sensibility elasticity parameter for fluid;
Gas-bearing horizon explanation module 6103, distinguishes gas-bearing horizon for utilizing the intersection of S-wave impedance SI and P-S wave velocity ratio γ to explain;
The codomain range determination module 6104 of gas-bearing reservoir, for utilizing the intersection of S-wave impedance SI and P-S wave velocity ratio γ to set up the X plot reflecting gas-bearing formation, determines the codomain scope of gas-bearing reservoir;
Wherein, the first relational expression is as follows respectively to the 4th relational expression:
AI＝Vp*ρ；
SI＝Vs*ρ；
γ＝Vp/Vs；
σ＝(Vp
2-2Vs
2)/2(Vp
2-Vs
2)；
In formula: Vp is velocity of longitudinal wave, Vs is shear wave velocity, and ρ is density, and AI is p-wave impedance, and SI is S-wave impedance, and γ is P-S wave velocity ratio, and σ is Poisson ratio; Described Poisson ratioσ represents the scale-up factor of object transverse strain and longitudinal strain, instruction change of fluid.
Preferably, consult Fig. 8, described seismic elastic parameter inverting section acquiring unit 620, specifically can comprise:
Partial stack seismic data volume acquisition module 6201, for according to zone of interest buried depth and geophone offset seismic data state of signal-to-noise far away, the geological data choosing maximum offset scope does common offset superposition, obtain at least three partial stack seismic data volumes, described three partial stack seismic data volumes comprise: nearly geophone offset superposition of data body, middle geophone offset superposition of data body, geophone offset superposition of data body far away;
Wavelet extraction module 6202, demarcating for utilizing acoustic wave train logging data to carry out repeatedly well shake respectively to earthquake nearly geophone offset superposition of data body, middle geophone offset superposition of data body and geophone offset superposition of data body far away, each superposition of data body being extracted to the wavelet of respective angles scope; The length of described wavelet is 100-200ms, and half length of wavelet can be divided exactly by earthquake sampling rate, during wavelet, window is 2-5 wavelet lengths doubly, avoids tomography near zone of interest;
Inverting module 6203, for inputting the earthquake nearly geophone offset superposition of data body of three different angles superpositions, middle geophone offset superposition of data body and geophone offset superposition of data body far away, and the wavelet of correspondence, solve full Knott-zoeppritz equation with acoustic wave train logging data constraint employing and carry out prestack Simultaneous Inversion, obtain p-wave impedance data volume AI, S-wave impedance data volume SI and density data body ρ and P-S wave velocity ratio Vp/Vs.
Preferably, described sensibility elasticity parameter intersection section acquiring unit 630, specifically may be used for utilizing intersection function, by the elastic parameter intersection section that gas-bearing reservoir sensibility elasticity parameter S-wave impedance and P-S wave velocity ratio intersection obtain, the region that intersection identifies meets well logging intersection region petrophysical parameter feature, finally obtains elastic parameter intersection section.
Preferably, described gas-bearing reservoir delineation unit 640, specifically may be used on sensibility elasticity parameter intersection section, drawing a circle to approve gas-bearing reservoir scope according to intersection region, pickup gas sand time thick T, calculates Effective Reservoirs thickness H according to the 6th relational expression; Described 6th relational expression is: H=V*T/2, and wherein, V is gas-bearing reservoir speed.
The advantage of apparatus of the present invention embodiment is: it is preferred to the present invention is based on sensibility elasticity parameter, the intersection of sensibility elasticity parameter is made an explanation and can predict gassiness (effectively) reservoir thickness quantitatively, this device can be predicted Effective Reservoirs thickness easily and fast, reduces the multi-solution of the parameter of prestack inversion unitary elasticity in the past qualitative forecasting reservoir.
In one or more exemplary design, the above-mentioned functions described by the embodiment of the present invention can realize in the combination in any of hardware, software, firmware or this three.If realized in software, these functions can store on the medium with computer-readable, or are transmitted on the medium of computer-readable with one or more instruction or code form.Computer readable medium comprises computer storage medium and is convenient to make to allow computer program transfer to the telecommunication media in other place from a place.Storage medium can be that any general or special computer can the useable medium of access.Such as, such computer readable media can include but not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other anyly may be used for carrying or store the medium that can be read the program code of form with instruction or data structure and other by general or special computer or general or special processor.In addition, any connection can be properly termed computer readable medium, such as, if software is by a concentric cable, optical fiber computer, twisted-pair feeder, Digital Subscriber Line (DSL) or being also comprised in defined computer readable medium with wireless way for transmittings such as such as infrared, wireless and microwaves from a web-site, server or other remote resource.Described video disc (disk) and disk (disc) comprise Zip disk, radium-shine dish, CD, DVD, floppy disk and Blu-ray Disc, and disk is usually with magnetic duplication data, and video disc carries out optical reproduction data with laser usually.Above-mentioned combination also can be included in computer readable medium.
Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (4)
1. based on a gas-bearing reservoir prediction method for many sensibility elasticities parameter, it is characterized in that, described method comprises:
According to the geophysical character of Effective Reservoirs acoustic wave train logging data, determine the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir;
Utilize Prestack seismic data to carry out prestack Simultaneous Inversion, obtain seismic elastic parameter inverting section;
Project on described seismic elastic parameter inverting section obtain sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir;
According to described sensibility elasticity parameter intersection section delineation gas-bearing reservoir scope, and pick up gas sand thickness;
Wherein, the described geophysical character according to Effective Reservoirs acoustic wave train logging data, determine that the sensibility elasticity parameter can distinguishing lithology and gas-bearing reservoir comprises:
Study area and many mouthfuls of full wave train log objective interval reservoir properties are around added up, calculate the elastic parameter of acoustic wave train logging data respectively to the 4th relational expression according to the first relational expression, described elastic parameter comprises: p-wave impedance AI, S-wave impedance SI, P-S wave velocity ratio γ and Poisson ratioσ;
Intersection analysis is carried out to described elastic parameter, determines that P-S wave velocity ratio γ, Poisson ratioσ are the sensibility elasticity parameter for fluid;
Utilize the intersection of S-wave impedance SI and P-S wave velocity ratio γ to explain and distinguish gas-bearing reservoir;
Utilize the intersection of S-wave impedance SI and P-S wave velocity ratio γ to set up the X plot reflecting gas-bearing formation, determine the codomain scope of gas-bearing reservoir;
Wherein, the first relational expression is as follows respectively to the 4th relational expression:
AI＝Vp*ρ；
SI＝Vs*ρ；
γ＝Vp/Vs；
σ＝(Vp
2-2Vs
2)/2(Vp
2-Vs
2)；
In formula: Vp is velocity of longitudinal wave, Vs is shear wave velocity, and ρ is density, and AI is p-wave impedance, and SI is S-wave impedance, and γ is P-S wave velocity ratio, and σ is Poisson ratio; Described Poisson ratioσ represents the scale-up factor of object transverse strain and longitudinal strain, instruction change of fluid;
The described Prestack seismic data that utilizes carries out prestack Simultaneous Inversion, obtains seismic elastic parameter inverting section and comprises:
According to zone of interest buried depth and geophone offset seismic data state of signal-to-noise far away, the geological data choosing maximum offset scope does common offset superposition, obtain at least three partial stack seismic data volumes, described three partial stack seismic data volumes comprise: nearly geophone offset superposition of data body, middle geophone offset superposition of data body, geophone offset superposition of data body far away;
Utilize acoustic wave train logging data to carry out repeatedly well shake respectively to earthquake nearly geophone offset superposition of data body, middle geophone offset superposition of data body and geophone offset superposition of data body far away to demarcate, each superposition of data body is extracted to the wavelet of respective angles scope; The length of described wavelet is 100-200ms, and half length of wavelet can be divided exactly by earthquake sampling rate, during wavelet, window is 2-5 wavelet lengths doubly, avoids tomography near zone of interest;
Input the earthquake nearly geophone offset superposition of data body of three different angles superpositions, middle geophone offset superposition of data body and geophone offset superposition of data body far away, and each self-corresponding wavelet, solve full Knott-zoeppritz equation with acoustic wave train logging data constraint employing and carry out prestack Simultaneous Inversion, obtain p-wave impedance data volume AI, S-wave impedance data volume SI and density data body ρ and P-S wave velocity ratio Vp/Vs;
The described sensibility elasticity parameter intersection by gas-bearing reservoir obtains sensibility elasticity parameter intersection section and comprises:
Utilize intersection function, elastic parameter S-wave impedance, P-S wave velocity ratio intersection section is obtained by gas-bearing reservoir sensibility elasticity parameter S-wave impedance and P-S wave velocity ratio intersection, the region that intersection identifies meets well logging intersection region petrophysical parameter feature, finally obtains sensibility elasticity parameter intersection section.
2. method according to claim 1, is characterized in that, described according to described sensibility elasticity parameter intersection section delineation gas-bearing reservoir scope, and picks up gas sand thickness and comprise:
On sensibility elasticity parameter intersection section, draw a circle to approve gas-bearing reservoir scope according to intersection region, pickup gas sand time thickness T, calculate Effective Reservoirs thickness H according to the 6th relational expression; Described 6th relational expression is: H=V*T/2, and wherein, V is gas-bearing reservoir speed.
3. based on a gas-bearing reservoir prediction device for many sensibility elasticities parameter, it is characterized in that, described device comprises:
Well logging sensibility elasticity parameter determination unit, for the geophysical character according to Effective Reservoirs acoustic wave train logging data, determines the multiple sensibility elasticity parameters can distinguishing lithology and gas-bearing reservoir;
Seismic elastic parameter inverting section acquiring unit, for utilizing Prestack seismic data to carry out prestack Simultaneous Inversion, obtains seismic elastic parameter inverting section;
Sensibility elasticity parameter intersection section acquiring unit, for being projected on described seismic elastic parameter inverting section obtained sensibility elasticity parameter intersection section by the sensibility elasticity parameter intersection of gas-bearing reservoir;
Gas-bearing reservoir delineation unit, for drawing a circle to approve gas-bearing reservoir scope according to described sensibility elasticity parameter intersection section, and picks up gas sand thickness;
Wherein, described well logging sensibility elasticity parameter determination unit, comprising:
Elastic parameter computing module, for adding up study area and many mouthfuls of full wave train log objective interval reservoir properties around, calculate the elastic parameter of acoustic wave train logging data respectively to the 4th relational expression according to the first relational expression, described elastic parameter comprises: p-wave impedance AI, S-wave impedance SI, P-S wave velocity ratio γ and Poisson ratioσ;
Sensibility elasticity parameter determination module, for carrying out crossplot analysis and its to described elastic parameter, determines that P-S wave velocity ratio γ, Poisson ratioσ are the sensibility elasticity parameter for fluid;
Gas-bearing horizon explanation module, distinguishes gas-bearing horizon for utilizing the intersection of S-wave impedance SI and P-S wave velocity ratio γ to explain;
The codomain range determination module of gas-bearing reservoir, for utilizing the intersection of S-wave impedance SI and P-S wave velocity ratio γ to set up the X plot reflecting gas-bearing formation, determines the codomain scope of gas-bearing reservoir;
Wherein, the first relational expression is as follows respectively to the 4th relational expression:
AI＝Vp*ρ；
SI＝Vs*ρ；
γ＝Vp/Vs；
σ＝(Vp
2-2Vs
2)/2(Vp
2-Vs
2)；
In formula: Vp is velocity of longitudinal wave, Vs is shear wave velocity, and ρ is density, and AI is p-wave impedance, and SI is S-wave impedance, and γ is P-S wave velocity ratio, and σ is Poisson ratio; Described Poisson ratioσ represents the scale-up factor of object transverse strain and longitudinal strain, instruction change of fluid;
Described seismic elastic parameter inverting section acquiring unit, comprising:
Partial stack seismic data volume acquisition module, for according to zone of interest buried depth and geophone offset seismic data state of signal-to-noise far away, the geological data choosing maximum offset scope does common offset superposition, obtain at least three partial stack seismic data volumes, described three partial stack seismic data volumes comprise: nearly geophone offset superposition of data body, middle geophone offset superposition of data body, geophone offset superposition of data body far away;
Wavelet extraction module, demarcating for utilizing acoustic wave train logging data to carry out repeatedly well shake respectively to earthquake nearly geophone offset superposition of data body, middle geophone offset superposition of data body and geophone offset superposition of data body far away, each superposition of data body being extracted to the wavelet of respective angles scope; The length of described wavelet is 100-200ms, and half length of wavelet can be divided exactly by earthquake sampling rate, during wavelet, window is 2-5 wavelet lengths doubly, avoids tomography near zone of interest;
Inverting module, for inputting the earthquake nearly geophone offset superposition of data body of three different angles superpositions, middle geophone offset superposition of data body and geophone offset superposition of data body far away, and each self-corresponding wavelet, solve full Knott-zoeppritz equation with acoustic wave train logging data constraint employing and carry out prestack Simultaneous Inversion, obtain p-wave impedance data volume AI, S-wave impedance data volume SI and density data body ρ and P-S wave velocity ratio Vp/Vs;
Described sensibility elasticity parameter intersection section acquiring unit, specifically for utilizing intersection function, S-wave impedance, the intersection of P-S wave velocity ratio section is obtained by the sensibility elasticity parameter S-wave impedance of gas-bearing reservoir and P-S wave velocity ratio intersection, the region that intersection identifies meets well logging intersection region petrophysical parameter feature, finally obtains elastic parameter intersection section.
4. device according to claim 3, it is characterized in that, described gas-bearing reservoir delineation unit, specifically for drawing a circle to approve gas-bearing reservoir scope according to intersection region on sensibility elasticity parameter intersection section, pickup gas sand time thickness T, calculates Effective Reservoirs thickness H according to the 6th relational expression; Described 6th relational expression is: H=V*T/2, and wherein, V is gas-bearing reservoir speed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310072499.2A CN103176211B (en) | 2013-03-07 | 2013-03-07 | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310072499.2A CN103176211B (en) | 2013-03-07 | 2013-03-07 | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103176211A CN103176211A (en) | 2013-06-26 |
| CN103176211B true CN103176211B (en) | 2016-04-20 |
Family
ID=48636190
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310072499.2A Active CN103176211B (en) | 2013-03-07 | 2013-03-07 | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103176211B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104675392B (en) * | 2013-12-02 | 2017-11-24 | 中国石油化工股份有限公司 | A kind of reservoir lithology recognition methods based on prestack multi-parameter dimensionality reduction |
| CN104991272A (en) * | 2015-07-02 | 2015-10-21 | 河海大学 | Earthquake speed disturbance modeling method for well-free earthquake reversion |
| CN107462927A (en) * | 2016-06-03 | 2017-12-12 | 中国石油化工股份有限公司 | Seismic facies Forecasting Methodology and device based on Naive Bayes Classification |
| CN106353800B (en) * | 2016-08-12 | 2018-11-09 | 北京天工新源科技有限公司 | A kind of information extracting method of multiple wavelet atom dictionary interaction analysis |
| CN106707345B (en) * | 2016-12-13 | 2018-10-16 | 中国石油天然气股份有限公司 | A kind of angle elastic parameter Lithology Identification Methods and device |
| CN107831537B (en) * | 2017-01-05 | 2020-07-14 | 中国石油化工股份有限公司 | Method for describing high sand-to-land ratio nepheloid rock reservoir |
| CN108629459B (en) * | 2018-05-10 | 2022-05-10 | 中国石油天然气股份有限公司 | Method and device for detecting hydrocarbon-containing pore of reservoir |
| CN112014879A (en) * | 2019-05-31 | 2020-12-01 | 中国石油天然气集团有限公司 | Reservoir stratum prediction method and device based on mixed rock |
| CN113311483B (en) * | 2020-02-27 | 2024-04-30 | 中国石油天然气集团有限公司 | Pre-stack elastic parameter combination inversion method and system based on shale oil |
| CN114509812B (en) * | 2022-01-29 | 2024-05-10 | 中海石油(中国)有限公司上海分公司 | Clastic rock reservoir dessert area prediction method, device, medium and electronic equipment |
Citations (3)
| 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 |
| CN101329405A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Simple method of multi-parameter seismic inversion |
| CN101872024A (en) * | 2010-06-02 | 2010-10-27 | 中国海洋石油总公司 | Method for carrying out well design by using time-lapse seismic |
-
2013
- 2013-03-07 CN CN201310072499.2A patent/CN103176211B/en active Active
Patent Citations (3)
| 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 |
| CN101329405A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Simple method of multi-parameter seismic inversion |
| CN101872024A (en) * | 2010-06-02 | 2010-10-27 | 中国海洋石油总公司 | Method for carrying out well design by using time-lapse seismic |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103176211A (en) | 2013-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103176211B (en) | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter | |
| US8830788B2 (en) | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media | |
| Tang et al. | Processing array acoustic-logging data to image near-borehole geologic structures | |
| CN107526101B (en) | A kind of acquisition and processing method obtaining earthquake reflected wave | |
| US8406081B2 (en) | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers | |
| CA1256975A (en) | Geophysical exploration by interpretation of variations in seismic amplitudes | |
| NO337875B1 (en) | Method and medium for predicting free-surface multiples in 3D seismic data | |
| CN103576200B (en) | Low signal-to-noise ratio zone shallow wave impedance interface static correction | |
| MX2014003060A (en) | Systems and methods for frequency-domain filtering and space-time domain discrimination of seismic data. | |
| CN105388518A (en) | Centroid frequency and spectral ratio integrated borehole seismic quality factor inversion method | |
| CN101551466A (en) | Method for improving prediction precision of oil and gas reservoir by using seismic attribute related to offset distance | |
| CN104459777A (en) | Fluid identification method and system based on fluid bulk modulus AVO inversion | |
| WO2015092542A2 (en) | Time-lapse simultaneous inversion of amplitudes and time shifts constrained by pre-computed input maps | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN111722284A (en) | Method for establishing speed depth model based on gather data | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| Pafeng et al. | Prestack waveform inversion of three-dimensional seismic data—An example from the Rock Springs Uplift, Wyoming, USA | |
| CN104570071B (en) | A kind of linear AVA and AVF inversion methods of viscous acoustic medium Bayes | |
| CN113589375A (en) | VSP layer velocity inversion method based on inclined layer constraint travel time calculation | |
| Fei et al. | Interactive Multimedia Data Coscattering Point Imaging for Low Signal‐to‐Noise Ratio 3D Seismic Data Processing | |
| US20230184974A1 (en) | Systems and methods for reservoir characterization | |
| US20230184982A1 (en) | Full-wavefield angle gather for high-contrast inter thin-bed models | |
| US20240184007A1 (en) | Method and system for determining attenuated seismic time | |
| Yang et al. | The Inverse Fresnel Beam XSP-CDP Stack Imaging in Crosswell Seismic | |
| EP3087417B1 (en) | Methods and systems of joint inversion determination of surface-consistent amplitude corrections |
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 |