CN104820241B - A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background - Google Patents
A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background Download PDFInfo
- Publication number
- CN104820241B CN104820241B CN201510107455.8A CN201510107455A CN104820241B CN 104820241 B CN104820241 B CN 104820241B CN 201510107455 A CN201510107455 A CN 201510107455A CN 104820241 B CN104820241 B CN 104820241B
- Authority
- CN
- China
- Prior art keywords
- frequency
- target zone
- thickness
- seismic data
- compact reservoir
- 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
Abstract
The invention belongs to geophysical prospecting for oil technical field, and in particular to the unconventional compact reservoir Forecasting Methodology under a kind of strong reflection background, by forward modeling, establish target zone different-thickness sand body corresponding to different frequency feature, establish forward model storehouse；Time frequency analysis is carried out to the seismic data of target zone, obtains the frequency domain data of target zone；Frequency domain data has been obtained using previous step, has been contrasted with forward model, has determined target zone sand thickness, utilize above-mentioned steps, seismic data and model library are contrasted to obtain final Reservoir Analysis situation, time-frequency analysis technology is converted using GS, can preferably identify lengthSection more than 8m compact reservoirs, improve compact reservoir and predict that coincidence rate up to 76.6%, plays a significant role in fine and close oily well site deployment and horizontal well region are preferred.
Description
Technical field
The invention belongs to geophysical prospecting for oil technical field, and in particular to the unconventional cause under a kind of strong reflection background
Close method for predicting reservoir.
Background technology
The rock stratum big with country rock resistance difference such as the coal seam developed in subsurface formations, hydrocarbon source rock are formed by force on seismic profile
Reflection, reservoir adjacent thereto are influenceed by its strong reflection, and ground rings weak, prediction more difficulty.Forefathers are attempted at seismic data
Strong reflection is eliminated in reason, or carries frequency and handles, but effect is all undesirable.
It is particularly long for example the fine and close oily major developmental of Ordos Basin is in Triassic Yanchang Formation Chang 7 memberIt is fine and close
Reservoir is grown close to hydrocarbon source rock, reservoir top reflective by its underpartThe strong reflection of oil shale interface influences, and is presented on seismic profile
Trough feature, Yi JichangSand body reflection is covered in lengthIn the secondary lobe of oil shale reflection, conventional method or time-domain inverting side
Method is difficult to.
It is longSand body is comparatively developed in main abrasive band, shows as block thick sand body（" box-like " or " trapezoidal shape "）, frequency
Domain is shown as based on low frequency；The argillaceous siltstoue of edge development high impedance in river course, with sandstone, mud stone composition thin interbed（" sawtooth
Shape "）, frequency domain shown as based on high frequency.Therefore, trough feature is presented on time-domain seismic profile, it is longSand body thickness without
Method is differentiated.
The content of the invention
The purpose of the present invention is：In time-domain, what seismic amplitude represented is the overlapping response on a set of stratum, therefore, is differentiated
Rate is relatively low, but is corresponded on frequency domain, the stratum of different-thickness with the frequency of response.According to this feature, using high-resolution
Seismic data, is transformed into frequency domain by rate GS converter techniques, and compact reservoir thickness is predicted by the analysis to frequecy characteristic.
Therefore, the invention provides the unconventional compact reservoir Forecasting Methodology under a kind of strong reflection background, it is characterised in that
Comprise the following steps：
（1）By forward modeling, establish target zone different-thickness sand body corresponding to different frequency feature, establish forward model storehouse；
（2）Time frequency analysis is carried out to the seismic data of target zone, obtains the frequency domain data of target zone；
（3）Frequency domain data has been obtained using previous step, has been contrasted with forward model, has determined target zone sand body thickness
Degree.
Described step（1）In forward modeling, be accomplished by the following way：Multiple wells are selected to carry out GS time-frequency transformation analysis,
The thickness of the target zone corresponding to well selected by these is different.
The well selected during described forward modeling includes three classes, and I class well target zone sand body is thick massive feature, II class well
Target zone sand body is thick massive+thin interbed feature, and III class well target zone sand body is thin layer or thin interbed feature, passes through contrast
The frequecy characteristic of the interval of interest of this three classes well, establish spectrum signature corresponding to different compact reservoir thickness.
Described step（2）Time frequency analysis process it is as follows：GS time frequency analysis is carried out to the seismic data of target zone, it is determined that
CDP trace gathers, certain one of selection CDP trace gathers, initial frequency value is determined, the road seismic data is transformed into frequency domain, Ran Houyi
The secondary frequency domain data that all seismic datas are converted into target zone.
The determination of described initial frequency value is determined that the frequency range of actual seismic data rises by the frequency range of actual seismic data
Frequency values corresponding at beginning are initial frequency value.
Beneficial effects of the present invention：Unconventional compact reservoir prediction side under this strong reflection background provided by the invention
Method, by establishing forward model storehouse, seismic data and model library are contrasted to obtain final Reservoir Analysis situation, when converting using GS
Frequency analysis technology, it can preferably identify lengthSection more than 8m compact reservoirs, improve compact reservoir and predict coincidence rate up to 76.6%,
Played a significant role in fine and close oily well site deployment and horizontal well region are preferred.
The present invention is described in further details below with reference to accompanying drawing.
Brief description of the drawings
Fig. 1 is X1 wells lengthCompact reservoir predicts instance graph.
Fig. 2 is the Seismic forward section spectral contrast analysis chart of western 3 mouthfuls of 233 wellblock well.
Embodiment
In time-domain, what seismic amplitude represented is the overlapping response on a set of stratum, and therefore, resolution ratio is relatively low, but in frequency
Domain, the stratum of different-thickness correspond with the frequency of response.According to this feature, using high-resolution GS converter techniques,
Seismic data is transformed into frequency domain, by lengthCompact reservoir thickness is predicted in the analysis of band frequency feature.
High-resolution GS conversion is the extension using Morlet small echos as the continuous wavelet transform thought of wavelet.In high score
In resolution GS conversion, wavelet is made up of the product of monochromatic wave and Gaussian functions, the monochromatic wave in wavelet
Only make stretching in time-domain, and Gaussian functions are then stretched and translated.This point and continuous wavelet transform are not
With.In continuous wavelet transform, monochromatic wave carries out same flexible and translation with Gaussian functions.Become with continuous wavelet
Change, the time-frequency domain analysis method such as Instant Fourier Transform is compared, high-resolution GS conversion has the advantages of its is unique, such as：
The resolution ratio and frequency of the S-transformation of signal（That is yardstick）Relevant, at the same time, S-transformation result and its Fourier spectrums of signal are protected
Direct contact is held, wavelet need not meet admissibility condition etc..Therefore, high-resolution GS conversion can preferably reflect Hubei Province
The development characteristics of 7 sections of compact reservoirs of length under the hydrocarbon source rock strong reflection background of Er Duosi basins.Carried out in the past using this technology
Sequence is explained, by analyzing the frequecy characteristic of big jacket layer position, determines that the set stratum is progradation or the deposition for moving back the long-pending phase, according to this
To judge the development degree of sand body, the prediction that sand thickness in a certain axis of reflection is carried out with technique has not been reported.
Embodiment 1：
The present embodiment provides the unconventional compact reservoir Forecasting Methodology under a kind of strong reflection background, it is characterised in that including
Following steps：
（1）By forward modeling, establish target zone different-thickness sand body corresponding to different frequency feature, establish forward model storehouse；
Step（1）In forward modeling, be accomplished by the following way：Multiple wells are selected to carry out GS time-frequency transformation analysis, these institutes
The thickness of target zone corresponding to the well of selection is different.
The well selected during described forward modeling includes three classes, and I class well target zone sand body is thick massive feature, II class well
Target zone sand body is thick massive+thin interbed feature, and III class well target zone sand body is thin layer or thin interbed feature, passes through contrast
The frequecy characteristic of the interval of interest of this three classes well, establish spectrum signature corresponding to different compact reservoir thickness.
（2）Time frequency analysis is carried out to the seismic data of target zone, obtains the frequency domain data of target zone；
Step（2）Time frequency analysis process it is as follows：GS time frequency analysis is carried out to the seismic data of target zone, determines CDP roads
Collection（CDP trace gathers：In seismic data acquisition, when reflecting interface level, on survey line in different common-shot-gathers, it can look for
To different roads, they both are from some common ground on subsurface interface, and the point is referred to as common depth point or common reflection point, have altogether
Corresponding each recorded trace composition common-depth-point gather with depth pip is referred to as CDP trace gathers）, certain one of selection CDP trace gathers, really
Determine initial frequency value, the road seismic data is transformed into frequency domain, for the purpose of then all seismic datas are converted successively
The frequency domain data of layer.
The determination of initial frequency value determines by the frequency range of actual seismic data, the section start pair of the frequency range of actual seismic data
The frequency values answered are initial frequency value.
（3）Frequency domain data has been obtained using previous step, has been contrasted with forward model, has determined target zone sand body thickness
Degree.
When determining frequecy characteristic and sand thickness relation using forward modeling, to consider the data situation in actual area, some when
It is consistent in trend to wait result and the actual conditions of forward modeling, such as, low frequency component for it is main when, represent sand body comparative development, it is high
Thin interbed feature is represented based on frequency composition, but high and low frequency here is relative, the result and actual seismic data of forward modeling
Can difference in codomain scope.Therefore it is determined that during radio-frequency component corresponding to different sand thickness, forward modeling and actual money
Suspect combination.The area that this method is applied under strong reflection background, amplitude generic attribute is insensitive to reservoir thickness, in reservoir and hydrocarbon
Source rock and coal bed is not in same reflection lineups, the obvious regional application effect in sand shale impedance area point is more preferable.
Embodiment 2：
As shown in figure 1, it is X1 wells lengthCompact reservoir predicts example, and the present embodiment is said with reference to specific implementation result
It is bright, time-frequency analysis technology is converted using GS, can preferably identify lengthSection more than 8m compact reservoirs.Using this of the present invention
Unconventional compact reservoir Forecasting Methodology under strong reflection background, Ordos Basin lengthGround compact reservoir prediction coincidence rate reaches
76.6%, played a significant role in fine and close oily well site deployment and horizontal well region are preferred.
In Fig. 1, x1 wells are grown in figureSection advantage dominant frequency is long on 35-45Hz, 35Hz single-frequency sectionsResponse is obvious, bores
Preceding prediction lengthFine and close sand thickness is 10m, and real bore is 9.2m, and earthquake prediction is that real brill is identical preferable.
Embodiment 3：
When carrying out unconventional compact reservoir prediction using time frequency analysis, to be passed through according to the specific geology characteristic in different regions
Seismic forward simultaneously combines actual seismic data situation, determines the relation between different sand thickness and seismic spectrum.
Fig. 2 is the Seismic forward section spectral contrast analysis chart of western 3 mouthfuls of 233 wellblock well, it can be seen that X2 wells are grownDuan Zhimi sand thickness is 25.7m, and sand body is self-contained, long in forward modeling sectionThe dominant frequency of section is 25-45Hz；X3 well sitesDuan Zhimi sand thickness is 12.3m, there is block sand body, is also had long in thin interbed forward modeling sectionThe dominant frequency of section is 25-
55Hz；X4 wells are grownSand body agensis, fine and close sand thickness are 7.2m, are grown in forward modeling sectionThe dominant frequency of section is 30-
55Hz；Analyzed by the forward modeling of multiple wells, with reference to actual regional data situation, sum up the head of districtFine and close sand thickness is more than
It is long during 15mSection dominant frequency is 25-45Hz, and 30Hz single-shot section seismic responses are obvious；When compact reservoir thickness is 8-15m,
It is longSection dominant frequency is 25-55Hz, and the response of 35Hz-40Hz single-frequency section is obvious；It is longDuan Zhimi sand thickness is less than 8m
When, it is longSection dominant frequency is 30-55Hz, and the response of 40-45Hz single-frequency section is obvious.According to frequency spectrum corresponding to different fine and close sand bodies
Feature, realize the earthquake that area's hydrocarbon source rock radiates by force more than 8m compact reservoirs under background and effectively predict.
It is exemplified as above be only to the present invention for example, do not form the limitation to protection scope of the present invention, it is all
It is to be belonged to the same or analogous design of the present invention within protection scope of the present invention.
Claims (2)
1. the unconventional compact reservoir Forecasting Methodology under a kind of strong reflection background, it is characterised in that comprise the following steps：
(1) by forward modeling, establish target zone different-thickness sand body corresponding to different frequency feature, establish forward model storehouse；
(2) time frequency analysis is carried out to the seismic data of target zone, obtains the frequency domain data of target zone；
(3) the frequency domain data obtained using previous step, is contrasted with forward model, determines target zone sand thickness；
Forward modeling in described step (1), is accomplished by the following way：Multiple wells are selected to carry out GS time-frequency transformation analysis, these
The thickness of target zone corresponding to selected well is different；
The well selected during described forward modeling includes three classes, and I class well target zone sand body is thick massive feature, II class well purpose
Layer sand body is thick massive+thin interbed feature, and III class well target zone sand body is thin layer or thin interbed feature, by contrast this three
The frequecy characteristic of the interval of interest of class well, establish spectrum signature corresponding to different compact reservoir thickness；
The time frequency analysis process of described step (2) is as follows：GS time frequency analysis is carried out to the seismic data of target zone, determines CDP
Trace gather, certain one of selection CDP trace gathers, initial frequency value is determined, the road seismic data is transformed into frequency domain, then successively will
All seismic datas are converted into the frequency domain data of target zone.
2. the unconventional compact reservoir Forecasting Methodology under strong reflection background as claimed in claim 1, it is characterised in that：Described
The determination of initial frequency value determines by the frequency range of actual seismic data, frequency corresponding to the section start of the frequency range of actual seismic data
Value is initial frequency value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510107455.8A CN104820241B (en) | 2015-03-12 | 2015-03-12 | A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510107455.8A CN104820241B (en) | 2015-03-12 | 2015-03-12 | A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104820241A CN104820241A (en) | 2015-08-05 |
| CN104820241B true CN104820241B (en) | 2017-12-05 |
Family
ID=53730582
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510107455.8A Active CN104820241B (en) | 2015-03-12 | 2015-03-12 | A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104820241B (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007096850A2 (en) * | 2006-02-24 | 2007-08-30 | National University Of Ireland, Galway | Use of single-channel seismic-reflection profiling data usable for the quantitative determination of sediment physical properties of the seabed and the layered subbottom |
| CN102109613B (en) * | 2009-12-23 | 2012-11-14 | 中国石油天然气股份有限公司 | Method for defining effective thickness of target reservoir bed under complex geological conditions |
| BR112015011142A8 (en) * | 2012-11-21 | 2019-10-01 | Halliburton Energy Services Inc | time domain compensation processing method and apparatus and non-transient computer readable storage medium |
| CN104280770B (en) * | 2014-09-28 | 2017-01-11 | 中国石油大港油田勘探开发研究院 | Prediction method of compact transition rock reservoir stratum |
-
2015
- 2015-03-12 CN CN201510107455.8A patent/CN104820241B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104820241A (en) | 2015-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Castellaro | The complementarity of H/V and dispersion curves | |
| Sullivan et al. | Architectural analysis of deep-water outcrops: Implications for exploration and development of the Diana sub-basin, western Gulf of Mexico | |
| Burnett et al. | Application of spectral decomposition to gas basins in Mexico | |
| CN105445800A (en) | Thick sand body top differentiation lithologic reservoir identification method | |
| Naseer et al. | Porosity prediction of lower cretaceous unconventional resource play, south Indus Basin, Pakistan, using the seismic spectral decomposition technique | |
| CN104698496A (en) | Small sand body boundary identification method and small sand body space quantitative description method | |
| Naseer et al. | Continuous wavelet transforms of spectral decomposition analyses for fluvial reservoir characterization of Miano Gas Field, Indus Platform, Pakistan | |
| CN105116449A (en) | Method for identifying weak reflection reservoir | |
| Naseer | Spectral decomposition'application for stratigraphic-based quantitative controls on Lower-Cretaceous deltaic systems, Pakistan: significances for hydrocarbon exploration | |
| Martorana et al. | Seismo-stratigraphic model of “La Bandita” area in the Palermo Plain (Sicily, Italy) through HVSR inversion constrained by stratigraphic data | |
| Hanafy et al. | New insights on structure and stratigraphic interpretation for assessing the hydrocarbon potentiality of the offshore Nile Delta basin, Egypt | |
| Ray et al. | Attribute-assisted interpretation of deltaic channel system using enhanced 3D seismic data, offshore Nova Scotia | |
| Naseer et al. | Simulating the stratigraphy of meandering channels and point bars of Cretaceous system using spectral decomposition tool, Southwest Pakistan: Implications for petroleum exploration | |
| Naseer | Spectral decomposition-based quantitative inverted velocity dynamical simulations of early cretaceous shaly-sandstone natural gas system, Indus Basin, Pakistan: Implications for low-velocity anomalous zones for gas exploration | |
| CN104820241B (en) | A kind of unconventional compact reservoir Forecasting Methodology under strong reflection background | |
| 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 | |
| Asim et al. | Spectral decomposition application for analyzing the structure and the reservoir potential: a case study of Penobscot, Nova Scotia offshore, Canada | |
| Ahmad et al. | Detection of fluvial sand systems using seismic attributes and continuous wavelet transform spectral decomposition: case study from the Gulf of Thailand | |
| Assis et al. | On the applicability of seismic attributes in the description and interpretation of mass-transport deposits | |
| Soeparjadi et al. | PD 7 (2) Exploration Play Concepts in Indonesia | |
| Elkhodary et al. | Integrated potential field study on the subsurface structural characterization of the area North Bahariya Oasis, Western Desert, Egypt | |
| Zdanowski et al. | Dim and bright spots as indicators of the Zechstein Main Dolomite hydrocarbon reservoir in Poland | |
| Lottaroli et al. | West Africa Presalt Exploration: Background and Perspective of the West African “Tupi” myth | |
| Isaac et al. | Improvement of seismic data quality and recognition of fault discontinuities through seismic data conditioning applications: a case study of Issaran oil field, Gulf of Suez, Egypt | |
| Geerdes et al. | Spectral decomposition of 3D ground-penetrating radar data from an alluvial environment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |