CN109343116A - A kind of stratum deformation earthquake detection method of non-structural ge nesis - Google Patents
A kind of stratum deformation earthquake detection method of non-structural ge nesis Download PDFInfo
- Publication number
- CN109343116A CN109343116A CN201811508361.1A CN201811508361A CN109343116A CN 109343116 A CN109343116 A CN 109343116A CN 201811508361 A CN201811508361 A CN 201811508361A CN 109343116 A CN109343116 A CN 109343116A
- Authority
- CN
- China
- Prior art keywords
- deformation
- stratum
- structural
- stratum deformation
- nesis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 22
- 238000010276 construction Methods 0.000 claims abstract description 32
- 238000003708 edge detection Methods 0.000 claims abstract description 31
- 239000003027 oil sand Substances 0.000 claims abstract description 11
- 238000012512 characterization method Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 2
- 239000002131 composite material Substances 0.000 abstract description 2
- 238000005056 compaction Methods 0.000 description 7
- 238000000034 method Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010796 Steam-assisted gravity drainage Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 150000003839 salts Chemical class 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
- 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 invention belongs to stratum deformation earthquake detection methods, in particular to the stratum deformation earthquake detection method of a kind of non-structural ge nesis, characterized by the following steps: 1) tectonic level explain: being based on seismic data, carries out the explanation of layer position to destination layer, obtain destination layer structural map；2) edge detection: carrying out edge detection to the destination layer structural map, calculates construction variable gradient, obtains the edge detection results of destination layer tectonic level；3) qualitative, quantitatively characterizing: qualitative and quantitatively characterizing is carried out, the distribution of stratum deformation and the size of deformation quantity are obtained.It the composite can be widely applied to Canadian oil-sand reservoir and domestic and international offshore oilfield, there is certain promotion value.
Description
Technical field
The invention belongs to stratum deformation earthquake detection method, in particular to the stratum deformation earthquake of a kind of non-structural ge nesis is examined
Survey method.
Background technique
Canadian oil-sand assets have become the important overseas assets of CNOOC.Currently, the oil sands reserves of Canada 1/5 are big
Mostly use in-situ extraction technology.Steam assisted gravity drainage (SAGD) is at present using more one of technology.For this
Special oil-sand development scheme, many geologic(al) factors will affect its development effectiveness, and stratum deformation is one of them.Stratum deformation
Presence will lead to vapor chamber heat along stratum loss, substantially reduce development efficiency, therefore oil-sand is stored up in the prediction of stratum deformation
The Efficient Development of layer is most important.And the stratum deformation of oil-sand reservoir is mainly caused by karst collapse and differential compaction,
Discontinuously there is larger difference with the stratum of structural ge nesis.
Stratum deformation is in the nature a kind of noncoherent boundary on stratum.And the discontinuity on stratum, as tomography, reefs, salt dome,
In the presence of the exception boundaries such as river, crack, seismic data will appear seismic reflection discontinuity zone, and reflection is on seismic image
For marginal information.
Currently, the detection of information discontinuous for the stratum of structural ge nesis, such as tomography, crack, there are many more at
Ripe seismic properties, such as attribute of coherent data volume, differential body attribute, dip and azimuth attribute and edge enhanced properties.These attributes
Potential tectonic information in extractable seismic data, is sliced by attribute, seismic interpretation personnel is instructed to carry out fault plane rapidly
The work such as combination and Fractured Zone identification.
The detection technique of the discontinuous information in the existing stratum based on seismic data, such as attribute of coherent data volume, incline at differential body attribute
Angular range angle attribute and edge enhanced properties etc. are proposed for the stratum noncoherent boundary of structural ge nesis, for non-construction
Miniature scale stratum deformation effects caused by moving are undesirable.
The caused stratum deformation of the effects of for non-structural ge nesis, especially karst collapse and differential compaction, due to it
Planar distribution is not strong, and earthquake reflection line-ups bad break phenomenon is unobvious on section, and traditional seismic properties cannot be preferable
Portray this kind of stratum plane of deformation spread situation in ground.
How it to be based on seismic data, carry out non-structural ge nesis, in particular for karst collapse and the differential compaction origin cause of formation
Small scale, faint stratum deformation prediction are of great significance to oil-sand reservoir SAGD development efficiency is improved.
It can be seen that the detection for the discontinuous information in stratum of non-structural ge nesis needs a kind of detection method.
Summary of the invention
In view of the above-mentioned problems, the object of the present invention is to provide a kind of stratum deformation earthquake detection method of non-structural ge nesis,
This method is particular for small scale stratum deformation earthquake detection, especially with regard to a kind of karst collapse suitable for non-structural ge nesis
With the stratum deformation detecting method of the differential compaction origin cause of formation.
To achieve the above object, the present invention takes following technical scheme:
A kind of stratum deformation earthquake detection method of non-structural ge nesis comprising following steps:
1) tectonic level is explained: being based on seismic data, is carried out the explanation of layer position to destination layer, obtain destination layer structural map；
2) edge detection: carrying out edge detection to the destination layer structural map, calculates construction variable gradient, obtains destination layer
The edge detection results of tectonic level；
3) qualitative, quantitatively characterizing: carry out qualitative and quantitatively characterizing, obtain stratum deformation distribution and deformation quantity it is big
It is small.
Further, in the step 3), the qualitative characterization refers to by planar distribution and seismic profile analysis, over the ground
The deformation quantity of layer deformation carries out qualitative characterization.
Further, the qualitative characterization method the following steps are included:
Planar distribution: the structural map of the edge detection results of step 2) and step 1) is overlapped, and is developed according to oil-sand is influenced
Stratum deformation tolerance given threshold lower limit, extract target zone in meet construction Grads threshold stratum deformation, obtain stratum
Plane of deformation spread result；
It eliminates the false and retains the true: according to stratum plane of deformation spread as a result, further being fallen in conjunction with seismic profile deformation size and feature
The reliability of the shape changing detection result obtained in fact rejects false information, obtains the distribution of stratum deformation.
Further, the quantitatively characterizing refers to step 2) edge detection results through construction gradient magnitude analysis, right
The deformation quantity of stratum deformation carries out quantitatively characterizing, obtains the size of stratum deformation quantity.
The invention adopts the above technical scheme, which has the following advantages:
1, the present invention identifies the stratum deformation of non-structural ge nesis using the tectonic level variation of plane, solves and does not advise
Rule, the not analysis of the stratum deformation of the non-structural ge nesis of tomography substantially.
2, the present invention reflects that tectonic level changes speed using construction edge detection, identifies the stratum shape of non-structural ge nesis
The distribution of change and deformation quantity, it is intuitive and reliable.
3, the present invention is due to only carrying out edge detection to tectonic level, without directly carrying out attribute calculating to seismic data,
Computational efficiency greatly improves, and can obtain the planar distribution information of stratum discontinuity rapidly.
4, the present invention is not due to carrying out operation directly to seismic data, and calculated result is not by seismic data resolution
Limitation.
5, the present invention finds construction and changes relatively violent position by calculating construction gradient, for small scale, faint
Discontinuous information has stronger directive significance.
6, the present invention, which combines edge detection results with tectonic information, portrays stratum deformation data, can reduce stratum deformation
The multi-solution of explanation.
7, the present invention is by the not strong small scale of planar distribution caused by non-structural ge nesis, faint stratum deformation
Tectonic level mutation analysis is carried out, and quantitatively characterizing is carried out to deformation quantity, not only the distribution of available stratum deformation, but also can
To obtain the deformation quantity on stratum, computational efficiency is high, and result credibility is high.
8, Canadian oil-sand reservoir and domestic and international offshore oilfield be the composite can be widely applied to, there is certain promotion value.
Detailed description of the invention
Fig. 1 is the schematic diagram of the method for the present invention；
The step of Fig. 2 is the method for the present invention figure；
Fig. 3 is the destination layer structural map of embodiment；
Fig. 4 is the destination layer structural map edge detection results figure of embodiment；
Fig. 5 is the structural map and edge detection results overlapping figure of embodiment；
Fig. 6 is the stratum shape changing detection result figure of embodiment；
Fig. 7 is the quantitatively characterizing figure of the shape changing detection result of embodiment.
Specific embodiment
The present invention is described in detail below with reference to the accompanying drawings and embodiments.
As shown in Figure 1, a kind of stratum deformation earthquake detection method of non-structural ge nesis, non-using tectonic level variation identification
The stratum deformation of structural ge nesis.
Using the speed of construction edge detection reflection tectonic level variation, the stratum deformation of non-structural ge nesis is identified.
It is combined by destination layer tectonic level with construction edge detection, identifies the miniature scale stratum shape of non-structural ge nesis
Become, including qualitative and quantitatively characterizing is carried out to stratum deformation, qualitative characterization is carried out for the distribution of stratum deformation, calmly
Scale sign is the deformation quantity size progress for stratum deformation.
The following steps are included:
1) explanation of layer position is carried out to destination layer, obtains destination layer structural map；
2) edge detection is carried out to structural map, calculates construction variable gradient, obtains the edge detection results of tectonic level；
3) qualitative and quantitatively characterizing is carried out, the distribution of stratum deformation and the size of deformation quantity are obtained.
(1) by planar distribution and seismic profile analysis, qualitative characterization qualitative characterization: is carried out to deformation quantity.
(2) quantitatively characterizing: to step 2) edge detection results by construction gradient magnitude analysis, deformation quantity is quantified
Characterization.
As shown in Fig. 2, specifically includes the following steps:
S1. layer position is explained: being based on seismic data, is carried out the explanation of layer position to destination layer, obtain destination layer structural map；
S2. edge detection: carrying out edge detection to structural map, calculates construction variable gradient, obtains the edge of tectonic level
Testing result；
S3. qualitative characterization:
(1) planar distribution: the structural map of the edge detection results of step 2) and step 1) is overlapped, and is opened according to oil-sand is influenced
The stratum deformation tolerance given threshold lower limit of hair, obtains stratum plane of deformation spread result；
(2) it eliminates the false and retains the true: is obtained as a result, further implementing (1) in conjunction with seismic profile according to stratum plane of deformation spread
The reliability of the shape changing detection result；
S4. according to construction gradient magnitude, quantitative table quantitatively characterizing: is carried out to the deformation quantity of the edge detection results of step 2)
Sign.
Embodiment
The present embodiment is based on 3D seismic data, right according to the deformation of seismic profile lineups and tectonic level delta data pair
Non- structural ge nesis, the especially caused stratum deformation of karst collapse and the effects of differential compaction is tested and analyzed, including with
Lower step:
1) the layer position for carrying out certain density to destination layer based on 3D seismic data is explained, and carries out grid to interpretation horizon
Change, obtains the structural map of destination layer, obtain the high low value of tectonic level, as shown in figure 3, different colors represents different depths
It spends, stratum height is represented with color [RGB, (255,0,0)] (red), [RGB, (255,225,0)] (yellow) in the present embodiment,
[RGB, (160,32,240)] (purple), [RGB, (0,0,255)] (blue) represent that stratum is low, and the faster place of color change is said
Bright construction variation is faster.
2) edge detection is carried out to tectonic level data, calculates construction gradient, obtains destination layer tectonic level edge detection
As a result, obtain construction gradient value, the speed of construction variation is represented, as shown in figure 4, to represent different variations fast for different colors
Slowly, it is that construction variation is big that [RGB, (160,32,240)] (purple), [RGB, (0,0,255)] (blue), which represent construction gradient greatly,
Construction variation is represented than faster region；[RGB, (255,0,0)] (red), [RGB, (0,255,0)] (green), [RGB,
(255,225,0)] (yellow) represent that construction change of gradient is small, and the relatively slow region of construction variation is represented, construction is paid close attention to and becomes
Change fast place, is equivalent to the edge detection results of construction.
3) edge detection results are superimposed together with structural map, i.e., two figures stack display, according to research area
Matter condition and the stratum deformation tolerance lower limit for influencing oil-sand exploitation, setting construction Grads threshold, extraction target zone has certain
The stratum deformation of scale obtains stratum plane of deformation spread as a result, possible stratum deformation distribution is obtained, as shown in figure 5, emphasis
Concern construction changes fast place, possible karst deformation or stratum deformation region.
4) according to stratum plane of deformation spread as a result, further implementing deformation inspection in conjunction with seismic profile deformation size and feature
The reliability for surveying result rejects false information, as shown in fig. 6, the range of White curves delineation is the karst collapse identified
And the effects of differential compaction caused stratum deformation.
5) according to construction gradient size values, quantitatively characterizing is carried out to deformation quantity, obtains the size of stratum deformation quantity, such as Fig. 7
Shown, it is larger that [RGB, (255,0,0)] (red) and [RGB, (255,225,0)] (yellow) represents deformation quantity, and [RGB, (0,
255,0)] (green) takes second place, and [RGB, (0,0,255)] (blue) and [RGB, (160,32,240)] (purple) represents deformation quantity phase
To smaller.Deformation quantity is bigger, and stratum deformation is faster, and formation variation is than comparatively fast, obtaining the size of stratum deformation quantity.
The deformation origin cause of formation is analyzed according to deformation quantity size 5) and planar distribution rule 4) as a result, based on difference
The origin cause of formation classifies to deformation, obtains non-structural ge nesis, especially karst collapse and differential compaction etc., caused stratum deformation
Distribution and layer position deformation quantity size.
The various embodiments described above are merely to illustrate the present invention, wherein the structure of each component, connection type and manufacture craft etc. are all
It can be varied, all equivalents and improvement carried out based on the technical solution of the present invention should not exclude
Except protection scope of the present invention.
Claims (4)
1. a kind of stratum deformation earthquake detection method of non-structural ge nesis, it is characterised in that the following steps are included:
1) tectonic level is explained: being based on seismic data, is carried out the explanation of layer position to destination layer, obtain destination layer structural map；
2) edge detection: carrying out edge detection to the destination layer structural map, calculates construction variable gradient, obtains destination layer construction
The edge detection results of layer position；
3) qualitative, quantitatively characterizing: qualitative and quantitatively characterizing is carried out, the distribution of stratum deformation and the size of deformation quantity are obtained.
2. a kind of stratum deformation earthquake detection method of non-structural ge nesis as described in claim 1, it is characterised in that: the step
It is rapid 3) in, the qualitative characterization refers to through planar distribution and seismic profile analysis, carries out to the deformation quantity of stratum deformation qualitative
Characterization.
3. a kind of stratum deformation earthquake detection method of non-structural ge nesis as claimed in claim 2, it is characterised in that: described fixed
Property characterization method the following steps are included:
Planar distribution: the structural map of the edge detection results of step 2) and step 1) is overlapped, according to the ground for influencing oil-sand exploitation
Layer deformation tolerance given threshold lower limit extracts the stratum deformation for meeting construction Grads threshold in target zone, obtains stratum deformation
Planar distribution result；
It eliminates the false and retains the true: according to stratum plane of deformation spread as a result, further being implemented in conjunction with seismic profile deformation size and feature
The reliability of the shape changing detection result arrived rejects false information, obtains the distribution of stratum deformation.
4. a kind of stratum deformation earthquake detection method of non-structural ge nesis as described in claim 1, it is characterised in that: described fixed
Scale sign refers to step 2) edge detection results through construction gradient magnitude analysis, quantifies to the deformation quantity of stratum deformation
Characterization, obtains the size of stratum deformation quantity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811508361.1A CN109343116B (en) | 2018-12-11 | 2018-12-11 | Stratum deformation earthquake detection method for non-structural cause |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811508361.1A CN109343116B (en) | 2018-12-11 | 2018-12-11 | Stratum deformation earthquake detection method for non-structural cause |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109343116A true CN109343116A (en) | 2019-02-15 |
| CN109343116B CN109343116B (en) | 2020-08-04 |
Family
ID=65303735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811508361.1A Active CN109343116B (en) | 2018-12-11 | 2018-12-11 | Stratum deformation earthquake detection method for non-structural cause |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109343116B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114002739A (en) * | 2021-11-11 | 2022-02-01 | 中海石油(中国)有限公司 | Edge detection method, device and medium based on geometric non-parallel statistical attributes |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104658037A (en) * | 2015-03-04 | 2015-05-27 | 中国地质科学院矿产资源研究所 | Automatic extraction method for potential field structural framework |
| CN104991269A (en) * | 2015-06-04 | 2015-10-21 | 中国科学技术大学 | Quick full-waveform inversion method for edge guide and structural constraint |
| CN105510964A (en) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | Seismic recognition method of low-order strike-slip faults in complex structural areas |
-
2018
- 2018-12-11 CN CN201811508361.1A patent/CN109343116B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104658037A (en) * | 2015-03-04 | 2015-05-27 | 中国地质科学院矿产资源研究所 | Automatic extraction method for potential field structural framework |
| CN104991269A (en) * | 2015-06-04 | 2015-10-21 | 中国科学技术大学 | Quick full-waveform inversion method for edge guide and structural constraint |
| CN105510964A (en) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | Seismic recognition method of low-order strike-slip faults in complex structural areas |
Non-Patent Citations (1)
| Title |
|---|
| 宋建国 等: ""基于结构导向的梯度属性边缘检测技术"", 《地球物理学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114002739A (en) * | 2021-11-11 | 2022-02-01 | 中海石油(中国)有限公司 | Edge detection method, device and medium based on geometric non-parallel statistical attributes |
| CN114002739B (en) * | 2021-11-11 | 2024-01-26 | 中海石油(中国)有限公司 | Edge detection method, device and medium based on geometric non-parallel statistical attribute |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109343116B (en) | 2020-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109655896B (en) | Method and system for identifying and depicting karst ancient river channel | |
| NO318302B1 (en) | Method of identifying large-scale / small-scale structures and migration pathways in hydrocarbon-containing areas / reservoirs | |
| CN105866832B (en) | The recognition methods for the low-grade fault that Deep Tight Sandstone is hidden and device | |
| CN104977617A (en) | Reservoir fracture identification method and imaging logging reservoir fracture identification method | |
| ES2927319T3 (en) | Method and system of measurement of fracturability | |
| CN109425899A (en) | A kind of prediction technique and device of the distribution of carbonate rock fault belt | |
| Alzubaidi et al. | Automatic fracture detection and characterization from unwrapped drill-core images using mask R–CNN | |
| CN109655905A (en) | Based on AVAZ inverting crack prediction method and system before full azimuth stacked data | |
| Basso et al. | Three‐dimensional seismic geomorphology of paleokarst in the cretaceous macaé group carbonates, campos basin, Brazil | |
| Nicolaevich et al. | Applying full–azimuth angle domain pre-stack migration and AVAZ inversion to study fractures in carbonate reservoirs in the Russian Middle Volga region | |
| CN109343116A (en) | A kind of stratum deformation earthquake detection method of non-structural ge nesis | |
| Liang et al. | Fracture characterization using diffraction attributes in tight sandstone reservoirs: A case study from Keshen Gas Field, Tarim Basin | |
| Ibanez et al. | An integrated remote sensing and biomarker maturity parameter-based approach to trace petroleum migration in a complex fault zone lineament of the Brazilian Equatorial Margin | |
| Nworie et al. | Paleokarst features in the aptian carbonates of the barra velha formation, santos basin, brazil | |
| Cheng et al. | Multiscale fracture prediction technique via deep learning, seismic gradient disorder, and aberrance: Applied to tight sandstone reservoirs in the Hutubi block, southern Junggar Basin | |
| Meng et al. | A case study of complex carbonate reservoir connectivity analysis, Tarim Basin, China | |
| Ozkaya | Use of exclusion zones in mapping and modeling fracture corridors | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Refayee et al. | Fault and fracture detection in unconventional reservoirs: A Utica shale study | |
| Monier et al. | Delineation of reservoir channels by different seismic attributes and geobody extractions for robust volumetric estimation, Saffron Field, offshore Nile Delta, Egypt | |
| Bourdon et al. | Karst identification and impact on development plan | |
| Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data | |
| Zhang et al. | Remote sensing detection of heavy oil through spectral enhancement techniques in the western slope zone of Songliao Basin, China | |
| Heggland | Hydrocarbon migration and accumulation above salt domes—Risking of prospects by the use of gas chimneys | |
| Sun et al. | Complex carbonate reservoir characterization using diffraction-imaging data set |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |