CN109799531B - Fracture reservoir prediction method based on seismic frequency division coherence attribute - Google Patents
Fracture reservoir prediction method based on seismic frequency division coherence attribute Download PDFInfo
- Publication number
- CN109799531B CN109799531B CN201910211333.1A CN201910211333A CN109799531B CN 109799531 B CN109799531 B CN 109799531B CN 201910211333 A CN201910211333 A CN 201910211333A CN 109799531 B CN109799531 B CN 109799531B
- Authority
- CN
- China
- Prior art keywords
- frequency division
- coherence
- frequency
- grids
- seam
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000001427 coherent effect Effects 0.000 claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000005516 engineering process Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000013507 mapping Methods 0.000 abstract description 2
- 206010017076 Fracture Diseases 0.000 description 10
- 208000010392 Bone Fractures Diseases 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
A fracture reservoir prediction method based on seismic frequency division coherence attributes belongs to the field of oil and gas geophysical reservoir prediction and aims to provide a fracture reservoir analysis technology capable of accurately identifying fractures of different scales and predicting distribution rules. The method comprises the following steps: firstly, extracting frequency division coherence attributes of different frequencies of a target layer and performing enhancement processing; calculating the average value and variance of 4 vertexes of each grid according to the frequency division coherent attribute data of the target layer after the gridding processing; thirdly, counting the range of the variance value, taking the median of the range, counting the average value of all grids corresponding to the range, setting the threshold values for dividing the seam units and the non-seam units according to the counting result, and quickly and effectively obtaining the distribution condition of the cracks; fourthly, for the grids which are seam units in the frequency division coherent attribute data of two target positions with different frequencies, the average value of the vertexes of the nearest non-seam units is given to the 4 vertexes of the grid in higher frequency, and the part of the crack prediction result of the higher frequency, which is overlapped with the crack prediction result of the lower frequency, can be effectively removed; fifthly, the frequency division coherence attribute data of different frequencies of the target layer position processed by the steps are utilized to realize the prediction of the fractured reservoir through plane mapping. The method is based on the target horizon frequency division coherence attribute data, divides the seam units and the non-seam units, compares the distribution conditions of the seam units in different frequency data, removes the overlapped parts, and completes the accurate identification of the different-scale seams and the accurate prediction of the distribution rule.
Description
Technical Field
The invention relates to a fractured reservoir prediction problem in the field of oil and gas geophysical, in particular to a fractured reservoir prediction method based on post-stack three-dimensional seismic attribute analysis, which is used for providing a fractured reservoir prediction analysis technology for accurately identifying fractures with different scales and predicting distribution rules.
Background
The post-stack three-dimensional seismic attribute analysis is an important method for identifying and predicting cracks, and the cracks are detected by using the seismic attributes mainly according to the discontinuity of seismic waveforms, such as coherent bodies, curvature, dip angles and the like. Because the cracks have multi-scale property, how to accurately represent different cracks, and how to identify and predict the distribution rules of the cracks with different scales are always the key and difficult points of research. Many scholars at home and abroad propose some technologies and methods, such as Zhang Guangzhi (2011) which takes Curvelet transformation as a basis, gives different reconstruction coefficients in a Curvelet domain to obtain seismic data volumes which highlight different frequency bands and different directions, and then predicts fracture and crack development zones and trends thereof by combining a coherent body edge detection technology; guo Rui et al (2014) propose a fracture analysis method based on high frequency recovery and multi-scale decomposition, widen the frequency band of three-dimensional seismic data, decompose seismic channels into data volumes corresponding to different scales (corresponding to the reconstruction wavelet dominant frequency), analyze fractures by applying coherence and curvature attributes, and depict formation fracture development behaviors and distribution characteristics under different scales; wupeng et al (2011) utilize the coherent attribute technique of wavelet frequency division to interpret the features of high-precision cracks and geologic bodies. The method mainly comprises the steps of carrying out multi-scale decomposition on the post-stack three-dimensional seismic data, and then, combining a post-stack three-dimensional seismic attribute analysis technology related to cracks, identifying the cracks with different scales and predicting the distribution rule. However, the above methods have a problem: the crack prediction is carried out on the basis of the post-stack three-dimensional seismic data with different decomposition scales, and the obtained results have overlapped parts, so that the recognition of the cracks with different scales and the prediction of the distribution rule are inaccurate. Therefore, the method is based on the post-stack three-dimensional seismic attribute analysis technology, and the cracks with different scales and the prediction distribution rule are accurately identified on the basis of the fact that the cracks have the same distribution characteristics in the crack prediction results of the seismic data based on different decomposition scales.
Disclosure of Invention
The invention aims to provide a fracture reservoir prediction analysis technology based on seismic frequency division coherence attributes, which effectively solves the problems of accurately identifying and predicting distribution rules of fractures with different scales. Processing and analyzing the seismic data by adopting a frequency division interpretation technology and extracting frequency division attributes; obtaining frequency division coherence attribute bodies with different frequencies by combining with a seismic coherence attribute analysis technology; extracting frequency division coherence attributes of different frequencies of a target horizon by using the interpreted seismic horizon; gridding the frequency division coherence attribute of the target layer, statistically analyzing the peak average value and variance of all grids, and dividing the grid attribute according to the statistical analysis result; comparing the distribution conditions of the seam unit grids in the data of different frequencies of the target layer, finding out the overlapped part, namely, all seam units in the data of different frequencies of a certain grid, and processing the grid in the data of higher frequency; and predicting the fractured reservoir based on the frequency division coherence property of different frequencies of the processed target layer.
The method comprises the following specific steps:
(1) firstly, time-frequency analysis based on continuous wavelet transformation is carried out on the three-dimensional seismic data after stacking, an effective seismic frequency band range is analyzed and found out, and a series of single-frequency data bodies with different frequencies are generated in the range by utilizing a frequency domain filtering technology.
(2) On the obtained single-frequency data body, performing coherent attribute extraction by adopting a third generation coherent algorithm to obtain frequency division coherent attribute bodies with different frequencies;
(3) enhancing the frequency division coherence attribute bodies with different frequencies;
(4) extracting frequency division coherence attributes of corresponding horizons aiming at specific target horizons on the frequency division coherence attribute body processed in the step (3) by utilizing the interpreted seismic horizons;
(5) gridding the frequency division coherence attribute data of the target horizon, and calculating the average value and variance of 4 vertexes of each grid; counting the variances of the vertexes of all the grids to obtain the distribution range of the vertexes; taking a median value of the variance distribution range, finding out all grids with the variance as the median value, counting the peak average values of the grids, setting a threshold value according to the counting result, and when the peak average value of the grids is larger than the threshold value, determining the grids as seam units, otherwise, determining the grids as non-seam units;
(6) and comparing the distribution of seam units and non-seam units of the frequency division coherent attribute data of the target layer positions with two different frequencies, and assigning values to 4 vertexes of the grid in a higher frequency when the grid is in the seam units in the results of the two frequencies, wherein the assigned value is the average value of the vertexes of the non-seam units nearest to the grid.
(7) And performing crack prediction by using the frequency division coherence property of different frequencies of the processed target layer.
The invention relates to a fracture reservoir prediction method based on seismic frequency division coherence attributes, which has the following characteristics that:
(1) calculating the average value and variance of 4 vertexes of each grid for the gridded seismic frequency division coherence attribute data of the target horizon, counting the range of variance values, taking the median of the range, counting the average values of all grids corresponding to the median, setting the threshold values of a seam dividing unit and a non-seam dividing unit according to the counting result, and quickly and effectively obtaining the distribution condition of the cracks;
(2) the method comprises the steps of assigning values to 4 vertexes of a grid which is a seam unit in frequency division coherent attribute data of two target positions with different frequencies, wherein the 4 vertexes of the grid are assigned to be the average value of the vertexes of non-seam units closest to the grid, so that the part, which is overlapped with a crack prediction result based on a lower frequency seismic frequency division coherent attribute, in a crack prediction result based on the higher frequency seismic frequency division coherent attribute can be effectively removed, and the problems of inaccurate prediction of identification and distribution rules of cracks with different scales are solved
Detailed Description
A fracture reservoir prediction method based on seismic frequency division coherence attributes comprises the following specific implementation steps:
(1) extracting target horizon frequency division coherence attributes
Step 1: wavelet frequency division attribute extraction:
firstly, performing time-frequency analysis on seismic data by adopting continuous wavelet transform with Morlet wavelet as wavelet, and transforming a seismic data body from a time domain to a frequency domain;
calculating the frequency spectrum of all time points of each seismic channel, and rearranging according to the frequency to generate common-frequency seismic data;
step 2: calculating the single-frequency seismic data with different frequencies obtained in the step 1 by adopting a third generation coherence algorithm to obtain coherence attribute bodies of the single-frequency seismic data bodies with different frequencies;
and step 3: carrying out enhancement processing on the coherent attribute bodies with different frequencies obtained by calculation in the step 2;
and 4, step 4: extracting the frequency division coherence attribute of the corresponding horizon aiming at the specific target horizon on the frequency division coherence attribute body processed in the step 4 by utilizing the interpreted seismic horizon;
(2) crack prediction based on target horizon frequency division coherence attribute
Step 1: analyzing the distribution characteristics of cracks in the frequency division coherence attribute data of different frequencies of the target layer:
gridding frequency division coherent attribute data of different frequencies of a target horizon, and calculating the average value and variance of 4 vertexes of each grid;
secondly, counting the variances of the vertexes of all the grids to obtain the distribution range of the vertexes; taking the median of variance distribution range, finding out all the grids with variance as the median, making statistics on the peak average of the grids, and making statistics on the result
Setting a threshold, and when the peak average value of the grid is greater than the threshold, determining the grid as a seam unit, otherwise, determining the grid as a non-seam unit, and obtaining the distribution characteristics of the cracks in the frequency division coherence attribute data of different frequencies of the target layer;
step 2: and comparing the distribution of seam units and non-seam units of the frequency division coherent attribute data of two different frequencies of the target layer, and when the grid is both seam units in the two frequency results, assigning values to 4 vertexes of the grid in higher frequency, wherein the assigned value is the average value of the vertexes of the non-seam units closest to the grid of the unit. The distance between the center points of any two grids is calculated by the following formula:
in the formula: a. b is the vertical line number and the contact line number of the lower right vertex of the 1 st grid, i and j are the vertical line number and the contact line number of the lower right vertex of the 2 nd grid, delta x is the vertical line interval and delta y is the contact line interval.
And step 3: and obtaining the graphic effect of the crack prediction by utilizing the processed frequency division coherence attribute data of different frequencies of the target layer through planar mapping, thereby realizing the prediction of the crack reservoir.
Claims (1)
1. A fracture reservoir prediction method based on seismic frequency division coherence attributes is characterized by adopting a brand-new method, such as the following steps:
(1) taking the three-dimensional seismic data after stacking as a research object, performing time-frequency analysis on the seismic data by adopting continuous wavelet transform with wavelet being Morlet wavelet, transforming the seismic data from a time domain to a frequency domain, generating a series of single-frequency data bodies with different frequencies by using a frequency domain filtering technology, and obtaining seismic frequency division coherence attribute bodies with different frequencies by adopting a third generation coherence algorithm;
(2) enhancing the seismic frequency division coherence attribute bodies with different frequencies, and extracting the frequency division coherence attribute of the corresponding layer by using the interpreted seismic layer;
(3) performing gridding processing on the frequency division coherence attribute data of the target layer, calculating the average value and variance of 4 vertexes of each grid, counting the variances of the vertexes of all the grids to obtain the distribution range of the grids, taking the median of the range, finding out the grids with the variance as the value, counting the vertex average values of the grids, setting a threshold value according to the counting result, and when the vertex average value of the grids is greater than the threshold value, determining the grids as a seam unit, otherwise, determining the grids as a non-seam unit;
(4) comparing the distribution of seam units and non-seam units of target horizon frequency division coherent attribute data with two different frequencies, when the grid is seam units in the results of the two frequencies, assigning values to 4 vertexes of the grid in the higher frequency, wherein the assigned value is the average value of the vertexes of the non-seam units nearest to the grid, and calculating the distance between the central points of any two grids by adopting a formula (1)
In the formula: a. b is the longitudinal survey line number and the contact survey line number of the lower right vertex of the 1 st grid respectively, i and j are the longitudinal survey line number and the contact survey line number of the lower right vertex of the 2 nd grid respectively, Δ x is the longitudinal survey line interval, and Δ y is the contact survey line interval;
(5) and performing crack prediction by using the frequency division coherence property of different frequencies of the processed target layer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910211333.1A CN109799531B (en) | 2019-03-20 | 2019-03-20 | Fracture reservoir prediction method based on seismic frequency division coherence attribute |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910211333.1A CN109799531B (en) | 2019-03-20 | 2019-03-20 | Fracture reservoir prediction method based on seismic frequency division coherence attribute |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109799531A CN109799531A (en) | 2019-05-24 |
CN109799531B true CN109799531B (en) | 2022-04-01 |
Family
ID=66563690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910211333.1A Expired - Fee Related CN109799531B (en) | 2019-03-20 | 2019-03-20 | Fracture reservoir prediction method based on seismic frequency division coherence attribute |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109799531B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112649863A (en) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | Frequency division seismic attribute data optimization method and system |
CN111239817B (en) * | 2020-02-05 | 2022-04-01 | 成都理工大学 | Method for improving fault likelihood attribute analysis precision |
CN111399056B (en) * | 2020-04-29 | 2020-12-08 | 西南石油大学 | Method for predicting crack strength based on divided azimuth filtering |
CN114252914A (en) * | 2020-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for determining distribution of fracture system |
CN112415586B (en) * | 2020-11-16 | 2023-10-03 | 北京孚梅森石油科技有限公司 | Slurry leakage risk assessment method and device |
CN116520419B (en) * | 2023-05-05 | 2024-01-30 | 中国矿业大学(北京) | Hot fluid crack channel identification method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7254091B1 (en) * | 2006-06-08 | 2007-08-07 | Bhp Billiton Innovation Pty Ltd. | Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs |
CN103777244A (en) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | Method for quantitatively analyzing earthquake crack attribute volume |
CN109117851A (en) * | 2018-07-06 | 2019-01-01 | 航天星图科技(北京)有限公司 | A kind of video image matching process based on lattice statistical constraint |
CN109389172A (en) * | 2018-10-11 | 2019-02-26 | 中南大学 | A kind of radio-signal data clustering method based on printenv grid |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0655861B1 (en) * | 1993-11-26 | 2000-08-02 | Koninklijke Philips Electronics N.V. | Image composition method and imaging apparatus for performing said method |
US6407404B1 (en) * | 1999-03-15 | 2002-06-18 | Denso Corporation | Apparatus for the examining defect of monolithic substrate and method for examining the same |
JP4747219B2 (en) * | 2009-04-08 | 2011-08-17 | キヤノン株式会社 | Image processing apparatus and image processing method |
JP5511468B2 (en) * | 2009-04-13 | 2014-06-04 | キヤノン株式会社 | Image processing apparatus, control method therefor, and program |
CN102306376B (en) * | 2009-11-03 | 2014-12-17 | 蒋慧琴 | Method for adaptive medical image enhancement |
CN104749617B (en) * | 2013-12-26 | 2017-05-31 | 中国石油化工股份有限公司 | A kind of multi-scale facture reservoir forward method for establishing model |
US9881413B2 (en) * | 2016-06-09 | 2018-01-30 | Adobe Systems Incorporated | Illumination-guided example-based stylization of 3D renderings |
US9905054B2 (en) * | 2016-06-09 | 2018-02-27 | Adobe Systems Incorporated | Controlling patch usage in image synthesis |
CN108931814B (en) * | 2017-05-26 | 2020-10-20 | 中国石油化工股份有限公司 | Multi-attribute fusion based bedrock fracture prediction method |
CN107909083B (en) * | 2017-11-10 | 2019-03-15 | 中国科学院地理科学与资源研究所 | A kind of hough transform extracting method based on outline optimization |
-
2019
- 2019-03-20 CN CN201910211333.1A patent/CN109799531B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7254091B1 (en) * | 2006-06-08 | 2007-08-07 | Bhp Billiton Innovation Pty Ltd. | Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs |
CN103777244A (en) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | Method for quantitatively analyzing earthquake crack attribute volume |
CN109117851A (en) * | 2018-07-06 | 2019-01-01 | 航天星图科技(北京)有限公司 | A kind of video image matching process based on lattice statistical constraint |
CN109389172A (en) * | 2018-10-11 | 2019-02-26 | 中南大学 | A kind of radio-signal data clustering method based on printenv grid |
Also Published As
Publication number | Publication date |
---|---|
CN109799531A (en) | 2019-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109799531B (en) | Fracture reservoir prediction method based on seismic frequency division coherence attribute | |
CN103728659B (en) | A kind of method improving detection of karst cave precision | |
CN104142519B (en) | Mud rock crack oil deposit predicting method | |
CN105445800A (en) | Thick sand body top differentiation lithologic reservoir identification method | |
CN111158044B (en) | Buried hill fracture reservoir body oil reservoir prediction method | |
NZ603314A (en) | Windowed statistical analysis for anomaly detection in geophysical datasets | |
CN107255831A (en) | A kind of extracting method of prestack frequency dispersion attribute | |
CN103576191B (en) | A kind of method adopting seismic properties identification tomography | |
CN104345341A (en) | Region constraint-based frequency band division energy seismic surface wave processing method | |
CN111679275A (en) | Underground pipeline identification method based on ground penetrating radar | |
CN113887398A (en) | GPR signal denoising method based on variational modal decomposition and singular spectrum analysis | |
CN110954958A (en) | Crack and fault prediction method and system | |
CN104122590B (en) | A kind of gas-oil detecting method based on electromagnetic survey and system | |
CN116520419B (en) | Hot fluid crack channel identification method | |
Wu‐Yang et al. | Research and application of improved high precision matching pursuit method | |
CN112230275A (en) | Seismic waveform identification method and device and electronic equipment | |
CN112444871B (en) | Quantitative analysis method and equipment for crack spacing based on scattered wave seismic response characteristics | |
CN105277981A (en) | Non-consistent time-lapse seismic bin matching method based on wave field continuation compensation | |
CN110989032B (en) | Gravity horizontal total gradient fracture identification method based on inclination angle | |
Bravo et al. | Volume curvature attributes to identify subtle faults and fractures in carbonate reservoirs: Cimarrona Formation, Middle Magdalena Valley Basin, Colombia. | |
CN113075656A (en) | Method for comprehensively predicting coal mine working face fracture through ground penetrating radar and three-dimensional earthquake | |
Tronicke et al. | The redundant wavelet transform to process and interpret GPR data | |
CN105467447A (en) | Phased tendency energy matching seismic amplitude preservation evaluation method | |
CN110907994A (en) | Waveform inversion reservoir prediction method based on frequency division | |
CN113325477B (en) | Method, device and equipment for filtering scattering interference waves of pre-stack seismic data |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220401 |
|
CF01 | Termination of patent right due to non-payment of annual fee |