CN106707338B - A kind of lower reservoir high-precision forecasting method of strong shielding - Google Patents
A kind of lower reservoir high-precision forecasting method of strong shielding Download PDFInfo
- Publication number
- CN106707338B CN106707338B CN201611025132.5A CN201611025132A CN106707338B CN 106707338 B CN106707338 B CN 106707338B CN 201611025132 A CN201611025132 A CN 201611025132A CN 106707338 B CN106707338 B CN 106707338B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- seismic
- strong
- shielded layer
- seismic data
- 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
- 238000013277 forecasting method Methods 0.000 title claims abstract description 17
- 230000003068 static effect Effects 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of strong lower reservoir high-precision forecasting methods of shielding, the described method comprises the following steps: picking up the seismic horizon for the strong reflection shielded layer being located above reservoir on seismic profile；Static time shift is carried out to seismic data using the seismic horizon of pickup；The seismic data after static time shift is pre-processed using median filtering, to remove strong reflection shielded layer；Seismic data after median filtering move static inverse time, the seismic profile after obtaining removal strong reflection shielded layer；Seismic properties are extracted to the seismic data after removal strong reflection shielded layer, and carve reservoir, obtain reservoir prediction result.According to the method for the present invention carry out it is strong shield lower reservoir prediction, significantly improve by the RESERVOIR RECOGNITION ability above covering strong shielding action and being influenced, reduce the multi-solution of seismic data RESERVOIR INTERPRETATION, improve reservoir prediction precision under high-strength shielding.
Description
Technical field
The present invention relates to geophysical signals to explain that field is more particularly related to a kind of lower reservoir high-precision of strong shielding
Prediction technique.
Background technique
With going deep into for exploration and development, gas reservoir geological condition becomes increasingly complex, and reservoir prediction required precision is higher and higher.When
Special formation or lithology are contained when showing as strong reflection on seismic profile in reservoir superstratum, can make underlying strata reservoir information
It is blanked, to interfere RESERVOIR RECOGNITION and influence precision of prediction, brings bigger difficulty to gas reservoir development.
Summary of the invention
It is above-mentioned existing in the prior art it is an object of the present invention to solving for the deficiencies in the prior art
One or more problems.For example, one of the objects of the present invention is to provide a kind of strong sides for shielding lower reservoir prediction precision of raising
Method.
To achieve the goals above, the present invention provides a kind of strong lower reservoir high-precision forecasting method of shielding, the methods
The following steps are included: picking up the seismic horizon for the strong reflection shielded layer being located above reservoir on seismic profile；Utilize the ground of pickup
It shakes layer position and static time shift is carried out to seismic data；The seismic data after static time shift is pre-processed using median filtering, with
Remove strong reflection shielded layer；Seismic data after median filtering move static inverse time, after obtaining removal strong reflection shielded layer
Seismic profile；Seismic properties are extracted to the seismic data after removal strong reflection shielded layer, and carve reservoir, obtain reservoir prediction
As a result.
Strong one embodiment for shielding lower reservoir high-precision forecasting method according to the present invention, the reservoir prediction result is carving
The reservoir development position at quarter.
Strong one embodiment for shielding lower reservoir high-precision forecasting method according to the present invention, the extraction seismic properties step
Extract earthquake instantaneous amplitude attribute.
Strong one embodiment for shielding lower reservoir high-precision forecasting method according to the present invention, the static time shift step is by ground
Data are shaken according to the seismic horizon time adjustment of the pickup, come into line lineups according to time orientation.
Compared with prior art, the beneficial effect comprise that significantly improving is influenced by strong shielding action is above covered
RESERVOIR RECOGNITION ability, reduce the multi-solution of seismic data RESERVOIR INTERPRETATION, improve reservoir prediction precision under high-strength shielding.
Detailed description of the invention
By the description carried out with reference to the accompanying drawing, above and other purpose of the invention and feature will become more clear
Chu, in which:
Fig. 1 is the strong flow chart for shielding lower reservoir high-precision forecasting method of exemplary embodiment of the present.
Fig. 2 is original seismic cross-section.
Fig. 3 is to be cutd open using the lower reservoir high-precision forecasting method of strong shielding of exemplary embodiment of the present to the earthquake of Fig. 2
Face carries out static time shift and the seismic cross-section after coming into line.
Fig. 4 is to be cutd open using the lower reservoir high-precision forecasting method of strong shielding of exemplary embodiment of the present to the earthquake of Fig. 3
Face carries out the seismic cross-section after median filter process.
Fig. 5 is to be cutd open using the lower reservoir high-precision forecasting method of strong shielding of exemplary embodiment of the present to the earthquake of Fig. 4
Face carries out static reversed time shift treated seismic cross-section.
Fig. 6 is the mathematical principle figure of median filtering.
Fig. 7 is original section earthquake instantaneous amplitude attributed graph.
Fig. 8 is after removing strong shielded layer using the lower reservoir high-precision forecasting method of strong shielding of exemplary embodiment of the present
Earthquake instantaneous amplitude attributed graph.
Description of symbols:
1- reservoir position, T1- pick up layer position, the reservoir position that T2- is carved out.
Specific embodiment
Hereinafter, it is high that reservoir under strong shielding according to the present invention will be described in detail in conjunction with attached drawing and exemplary embodiment
Accuracy prediction method.
Using wave field separation technology, strong shielded layer seismic reflection is peeled off from original seismic profile, makes the weak reflection of reservoir
Signal is highlighted, and is of great significance to reservoir prediction precision is improved.
In order to improve reservoir prediction precision under strong shield, the present invention provides one kind and shields lower reservoir high-precision forecast side by force
Method, comprising the following steps: firstly, picking up on seismic profile, strong reflection shielded layer seismic horizon above reservoir；Then, in utilization
The layer position that portion is picked up carries out static time shift to by seismic data, is pre-processed using median filter method to seismic data, shells
From falling (in other words remove) strongly reflecting layer；Secondly, move static inverse time to data after median filtering, removal strong reflection screen is obtained
Seismic profile after covering layer；In turn, seismic properties are extracted to the seismic data after removal strong reflection shielded layer, carves reservoir.Most
Afterwards, reservoir prediction result is exported.
Fig. 1 is the strong flow chart for shielding lower reservoir high-precision forecasting method of exemplary embodiment of the present.Such as Fig. 1 institute
Show, reservoir high-precision forecasting method under strong shielding according to an exemplary embodiment of the present invention the following steps are included:
(1) it picks up on seismic profile, strong reflection shielded layer seismic horizon above reservoir, as shown in Fig. 2, explaining line T1 in figure
The strong shielded layer in reservoir top is represented, polygon frame 1 is reservoir position, and T1 layers of presence has seriously affected the RESERVOIR INTERPRETATION that underlies.
(2) the layer position picked up using upper step static time shift and come into line to by seismic data, and static time shift is with coming into line
Lineups will be picked up on seismic data by time shift to be aligned in time.T1 is come into line into the same time.It static time shift and comes into line
Result afterwards is as shown in Figure 3.
(3) seismic data is pre-processed using median filter method, peels off strongly reflecting layer, median filter process knot
Fruit is as shown in Figure 4.
Median filtering: being a kind of nonlinear filtering, and mathematical principle is as shown in Figure 6, it is assumed that has data sequence an xi, i=
1,2,3 ..., m, if the span of median filtering is n, the median filtering process to jth point is: taking using jth point as the n at midpoint
A sample value is as input；To this n sample value order rearrangement by size；Taking the sample value of n data center after resetting to be used as should
The filtering output of point.
(4) data after median filtering move static inverse time, the seismic profile after obtaining removal strong reflection shielded layer is such as
Shown in Fig. 5, Fig. 7 is original section instantaneous amplitude attribute.
(5) earthquake instantaneous amplitude attribute is extracted to the seismic data after removal strong reflection shielded layer, carves reservoir, sees Fig. 8
It is shown.
Earthquake instantaneous amplitude attribute is shown in Fig. 7 and Fig. 8 after comparing original and removal strong reflection shielded layer, it can be seen that in Fig. 8
T2 reservoir position seismic amplitude value increases obviously, and reservoir signal is reinforced, and RESERVOIR RECOGNITION precision is improved.
(6) it exports, terminates.
Export the reservoir development position of engraving, comprising: abscissa X, ordinate Y can be directly used for reservoir development Favorable Areas
Cheng Tu.
As can be seen that carry out in this way it is strong shield lower reservoir prediction, significantly improve by above covering strong shielding action
The RESERVOIR RECOGNITION ability of influence, reduces the multi-solution of seismic data RESERVOIR INTERPRETATION, improves reservoir prediction essence under high-strength shielding
Degree.This method is mainly used in concealment lithologic deposit reservoir prediction.
In conclusion the present invention removes strong shielded layer using median filter process method, strong shielded layer is eliminated to storage of underliing
The influence of layer highlights reservoir characteristic with seismic properties Conjoint Analysis, preferably reflection reservoir distribution, realizes the lower storage of strong shielding
The high-precision forecast of layer, improves probing success rate.This method is efficient, strong applicability.
Although those skilled in the art should be clear above by combining exemplary embodiment to describe the present invention
Chu can carry out exemplary embodiment of the present invention each without departing from the spirit and scope defined by the claims
Kind modifications and changes.
Claims (3)
1. a kind of lower reservoir high-precision forecasting method of strong shielding, which is characterized in that the described method comprises the following steps:
Pick up the seismic horizon for the strong reflection shielded layer being located above reservoir on seismic profile；
Static time shift is carried out to seismic data using the seismic horizon of the pickup；
The seismic data after static time shift is pre-processed using median filtering, to remove strong reflection shielded layer；
Seismic data after median filtering move static inverse time, the seismic profile after obtaining removal strong reflection shielded layer；
Seismic properties are extracted to the seismic data after removal strong reflection shielded layer, and carve reservoir, obtain reservoir prediction as a result, its
In,
Seismic data according to the seismic horizon time adjustment of the pickup, is made lineups according to the time by the static time shift step
Direction comes into line.
2. reservoir high-precision forecasting method under strong shielding according to claim 1, which is characterized in that the reservoir prediction knot
Fruit is the reservoir development position of engraving.
3. reservoir high-precision forecasting method under strong shielding according to claim 1, which is characterized in that the extraction earthquake category
Property step extract earthquake instantaneous amplitude attribute.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611025132.5A CN106707338B (en) | 2016-11-18 | 2016-11-18 | A kind of lower reservoir high-precision forecasting method of strong shielding |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611025132.5A CN106707338B (en) | 2016-11-18 | 2016-11-18 | A kind of lower reservoir high-precision forecasting method of strong shielding |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106707338A CN106707338A (en) | 2017-05-24 |
| CN106707338B true CN106707338B (en) | 2018-12-07 |
Family
ID=58940130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611025132.5A Active CN106707338B (en) | 2016-11-18 | 2016-11-18 | A kind of lower reservoir high-precision forecasting method of strong shielding |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106707338B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107356967B (en) * | 2017-07-26 | 2019-04-12 | 西安交通大学 | A kind of compacting seismic data shields by force the sparse optimization method of interference |
| CN109343115B (en) | 2018-11-21 | 2019-12-03 | 成都理工大学 | A kind of gas-bearing reservoir depicting method based on Log-constrained |
| CN111694055B (en) * | 2019-03-14 | 2023-04-25 | 中国石油天然气股份有限公司 | Method and device for extracting and analyzing AVO attribute under strong reflection shielding |
| CN111505708B (en) | 2020-04-28 | 2021-04-20 | 西安交通大学 | Deep learning-based strong reflection layer stripping method |
| CN112965109B (en) * | 2021-02-05 | 2023-05-12 | 中国石油化工股份有限公司 | Method and device for predicting reservoir under strong reflection shielding |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8958267B2 (en) * | 2011-05-13 | 2015-02-17 | Conocophillips Company | Seismic true estimated wavelet |
| CN102809762B (en) * | 2012-08-13 | 2014-09-03 | 成都理工大学 | Reservoir imaging technique based on full-frequency-band seismic information mining |
| CN103412332B (en) * | 2013-01-22 | 2016-05-25 | 中国地质大学（北京） | A kind of method of determining thin reservoir thickness |
| CN104849758B (en) * | 2015-05-05 | 2017-10-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | For the drawing method of the rule interference in geological data |
| CN105259579A (en) * | 2015-11-10 | 2016-01-20 | 中国地质大学（北京） | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes |
| CN105467451B (en) * | 2016-01-13 | 2018-05-15 | 中国石油集团东方地球物理勘探有限责任公司 | Fractal inversion method based on total variation minimization constraint |
| CN105911585B (en) * | 2016-07-05 | 2018-05-15 | 中国石油集团东方地球物理勘探有限责任公司 | A kind of extracting method and device of earthquake record regular interference |
-
2016
- 2016-11-18 CN CN201611025132.5A patent/CN106707338B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106707338A (en) | 2017-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106707338B (en) | A kind of lower reservoir high-precision forecasting method of strong shielding | |
| CN105510976A (en) | Multiple-wave combined self-adaptive attenuation method | |
| Long et al. | Fast-AIC method for automatic first arrivals picking of microseismic event with multitrace energy stacking envelope summation | |
| CN106324669B (en) | Method of separating surface-related multiples of different orders in seismic exploration data | |
| Mints et al. | A three-dimensional model of the Early Precambrian crust under the southeastern Fennoscandian Shield: Karelia craton and Belomorian tectonic province | |
| Lu et al. | Instantaneous polarization filtering focused on suppression of surface waves | |
| CN108828664A (en) | Multiple wave identification method and device | |
| CN108427140B (en) | A method of for small scale fracture and cave reservoir seismic recognition | |
| Li et al. | Experimental analysis and application of sparsity constrained deconvolution | |
| Gómez Londoño et al. | Using the Karhunen-Loeve transform to suppress ground roll in seismic data | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section | |
| CN108254785A (en) | Reservoir Determination Method and Device | |
| CN104424393B (en) | A kind of geological data reservoir reflectance signature based on principal component analysis strengthens method | |
| CN106610506B (en) | Seismic exploration thin layer identification method | |
| CN108693558B (en) | Seismic data processing method and device | |
| Ma et al. | A Global and Multi-Scale Denoising Method Based on Generative Adversarial Network for DAS VSP Data | |
| CN113126157B (en) | Frequency wave number domain high-angle fracture extraction method and device, storage medium and equipment | |
| CN109212609A (en) | Near surface Noise Elimination method based on wave equation continuation | |
| CN111694055B (en) | Method and device for extracting and analyzing AVO attribute under strong reflection shielding | |
| CN109917459A (en) | Method, device and system for suppressing earthquake noise | |
| Cowan et al. | Screening kimberlite magnetic anomalies in magnetically active areas | |
| CN113126156A (en) | Method and device for extracting high-angle fracture in radon region, storage medium and equipment | |
| Gao et al. | Zero-offset VSP wavefield separation using two-step SVD method | |
| Ma et al. | An aeromagnetic denoising-decomposition-3D inversion approach for mineral exploration | |
| CN104297776B (en) | Method for improving excitation effect of desert area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180402 Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201113 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Patentee before: BGP Inc., China National Petroleum Corp. |