CN111897008B - Fracture grading prediction method based on seismic frequency division technology - Google Patents
Fracture grading prediction method based on seismic frequency division technology Download PDFInfo
- Publication number
- CN111897008B CN111897008B CN202010790656.3A CN202010790656A CN111897008B CN 111897008 B CN111897008 B CN 111897008B CN 202010790656 A CN202010790656 A CN 202010790656A CN 111897008 B CN111897008 B CN 111897008B
- Authority
- CN
- China
- Prior art keywords
- fracture
- data
- frequency division
- seismic
- prediction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005516 engineering process Methods 0.000 title claims abstract description 17
- 206010017076 Fracture Diseases 0.000 claims abstract description 127
- 208000010392 Bone Fractures Diseases 0.000 claims abstract description 107
- 238000011161 development Methods 0.000 claims abstract description 25
- 238000005553 drilling Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000011160 research Methods 0.000 claims abstract description 9
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 9
- 230000005856 abnormality Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 5
- 230000005012 migration Effects 0.000 claims description 5
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 238000013480 data collection Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 4
- 230000002159 abnormal effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 208000013201 Stress fracture Diseases 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011295 pitch Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
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/282—Application of seismic models, synthetic seismograms
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
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 discloses a fracture grading prediction method based on a seismic frequency division technology, which comprises the following steps: collecting basic data; making a single-well synthetic record, and calibrating a research target layer; analyzing the original post-stack amplitude data frequency spectrum of the target layer to determine the effective signal range of the target layer data; defining the approximate scale of the target horizon fracture development; establishing a forward modeling for simulation; performing frequency division processing on a simulation result, and preferably selecting fracture response frequencies of different orders; performing targeted frequency division processing on the original amplitude data; analyzing each level of fracture properties and optimizing an algorithm; analyzing fracture horizontal splitting and three-dimensional space prediction effects at all levels; fusing the fractures of all levels into a graph, and performing well drilling data matching analysis; and eliminating the fracture detection abnormal area, and completing the fracture grading prediction. The invention has low requirement on data processing conditions, small operand and high processing speed; the problem of multiple solution of a large number of fracture interpretations is solved, the complex problem is simplified, the fracture release efficiency is improved, the workload of seismic interpretation is reduced, and the accuracy of fracture prediction is increased.
Description
Technical Field
The invention relates to the technical field of petroleum exploration and development, in particular to a fracture grading prediction method based on a seismic frequency division technology.
Background
At present, fracture (crack) prediction methods popular at home and abroad mainly include three types of recognition prediction methods of geology, well logging and earthquake, and from the recognition scale, the traditional geology and well logging recognition method mainly recognizes and predicts cracks with a certain distribution range and a smaller scale by means of data such as actual drilling, outcrop, rock core, well logging and the like, and the prediction range is limited.
Due to the influence of factors such as regional tectonic events, fracture development stages, strength and the like, the fractures in the research region are generally divided into main-stem (large-scale) fractures, secondary (medium-scale) fractures and micro (small-scale) fractures, and the prominent manifestations of the fractures in seismic data are that the fracture pitches are different in size and form. With the rapid development of three-dimensional seismic technology, crack identification methods based on earthquake are increasing and can be divided into prediction methods based on pre-stack data and post-stack data, and the former has relatively harsh processing conditions on data and relatively high cost, so that various crack earthquake detection attributes based on post-stack data are the main means for crack prediction. At present, the mainstream post-stack fracture prediction methods are all used for carrying out prediction work on full-frequency-band seismic data. The basic method comprises the following steps: firstly, seismic data filtering (structure-oriented filtering), then carrying out optimization work of methods such as waveform difference, coherence, edge detection, dip angle, curvature property, ant body, maximum likelihood property and the like, and finally completing fracture prediction. That is to say, for fractures with different geological scales, a certain preferred prediction method is uniformly used, so that the differences of the response of seismic data to fractures with different grades under different frequency conditions and the differences of the response of different algorithms to fractures with different grades are ignored, the research on fracture grades is weakened, and the deviation is not avoided.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a fracture grading prediction method based on the seismic frequency division technology, which perfects the problems in the technology.
In order to realize the purpose, the technical scheme adopted by the invention is as follows:
a fracture grading prediction method based on a seismic frequency division technology comprises the following steps:
s1, basic data collection: collecting data from the area of interest, the data comprising: data such as zone adjustment data, post-stack migration seismic body, conventional logging curve, drilling and logging data, well deviation data, coring data, layering data, seismic interpretation horizon, single well testing productivity and the like;
s2, making a single-well synthetic record, and calibrating a research target layer and an upper layer and a lower layer;
s3, analyzing the frequency spectrum of the original post-stack amplitude data of the target layer to determine the effective signal range of the data of the target layer;
s4, directly obtaining the approximate scale of the fracture development of the target horizon according to relevant geological data such as actual drilling, well logging, exposure and zone adjustment;
s5, designing a theoretical forward model according to the fault development scale, and setting forward model parameters according to original amplitude data volume parameters;
s6, frequency division processing is carried out on the forward modeling earthquake result, section identifiability of different-grade fracture is observed under different frequency data, and optimal response frequency of different fracture grades is selected;
s7, carrying out frequency division processing of the preferred frequency in S6 on the data body within the effective signal range of the original amplitude data body to obtain data bodies of a plurality of frequency bands;
s8, analyzing each level of fracture property and optimizing an algorithm;
s9, analyzing the fracture bisection and three-dimensional space prediction effect of each level;
and S10, fusing all levels of fractures into a graph, and matching and performing coincidence analysis on data related to fractures (fractures) such as actual drilling fracture and drilling conditions, fracture development conditions, reservoir development conditions, single well testing productivity and the like.
S11, eliminating the detection abnormality caused by the data, and completing the classification prediction of the fracture.
Further, the slice acquisition parameter in S4 is preceded by a fault distance.
Further, in S5, the forward modeling is to ensure that the parameters of the obtained seismic forward data and the original amplitude data are consistent as much as possible, and mainly includes: dominant frequency, bandwidth, etc.
Further, in said S7, the inventor proposes to perform frequency division processing on the original amplitude data volume by a frequency division technique (frequency magnifying glass) based on fast matching pursuit.
Further, in the S8, in the crack property and algorithm optimization, the inventor proposes that the low frequency volume is used for the identification and prediction calculation of the large-scale fracture, the medium frequency volume is used for the identification and prediction calculation of the medium-scale fracture, and the high frequency volume is used for the identification and prediction calculation of the small-scale or micro fracture (crack). Waveform difference attribute computation may be used for large scale fractures, coherence, similarity, etc. for mesoscale fractures, curvature analysis, maximum likelihood, ant tracking, etc. for small scale or microfractures (fractures)
Furthermore, in the S9 plano-section effect and the three-dimensional spatial prediction analysis, it is noted that a background is formed by combining the basic structural features and fractures of the actual research area, and the comprehensive analysis of the original amplitude data volume is combined, so that the prediction misalignment caused by the abnormality generated in the frequency division processing is avoided.
Further, the S10 matches and fits the data related to fracture (fracture) such as actual drilling fracture situation, fracture development situation, reservoir development situation, single well testing productivity, etc., and needs to specifically analyze the influence of the fracture cause, reservoir type and fracture on oil and gas migration and accumulation in the actual research area.
Further, the S11 eliminates detection abnormality caused by the data, and fracture detection is often abnormal for a non-full coverage area of the three-dimensional data.
Compared with the prior art, the invention has the advantages that:
(1) the method is a comprehensive application and innovation of the prior art, the process can be realized under the existing seismic exploration technical condition, and the time period from development to application of the crack prediction attribute or algorithm is saved.
(2) The method adopts post-stack migration data as basic seismic data, and has smaller calculated amount and higher calculating speed compared with a pre-stack fracture (crack) prediction means.
(3) The overall thought of fracture-fracture system frequency division hierarchical prediction is adopted, fractures are hierarchically depicted, the primary and secondary relations of the fractures are clarified, for seismic interpreters, the problem of multiple resolutions of a large number of fracture interpretations is eliminated, the complex problem is simplified, the fracture splitting release efficiency is improved, and fracture flat-cutting and space depiction are clearer.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a diagram illustrating a spectrum analysis of a target horizon for original post-stack data according to an embodiment of the present invention;
FIG. 3 is a simulation section of the size of the fracture development of the target zone and the earthquake forward modeling of the embodiment of the invention;
FIG. 4 is a forward analog frequency division profile according to an embodiment of the present invention;
FIG. 5a is a schematic diagram of the spatial distribution of stem fracture in accordance with an embodiment of the present invention;
FIG. 5b is a graphical representation of a stem fracture profile of an embodiment of the present invention;
FIG. 5c is a secondary fracture space spread characteristic diagram according to an embodiment of the present invention;
FIG. 5d is a graphical representation of the secondary fracture profile spread characteristics of an embodiment of the present invention;
FIG. 5e is a graph of the spatial distribution of cracks in accordance with an embodiment of the present invention;
FIG. 5f is a graphical representation of the fracture profile spread characteristics of an embodiment of the present invention;
FIG. 6 is a graph of fracture plane spread characteristics and single well productivity matching for an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be further described in detail below with reference to the accompanying drawings by way of examples.
As shown in fig. 1, a fracture stage prediction method based on seismic frequency division technology includes the following steps:
s1, basic data collection: collecting data from the area of interest, the data comprising: zone adjustment data, a post-stack migration seismic body, a conventional logging curve, drilling and logging data, well deviation data, coring data, layering data, seismic interpretation horizon and single well testing productivity data;
s2, making a single-well synthetic record by using the acoustic wave and the density logging curve, and calibrating a top-bottom interface of a target layer and interfaces of upper and lower layers of the target layer;
s3, analyzing the original superposed amplitude data spectrum of the target layer, and determining the effective signal range of the target layer data, the dominant frequency of the target layer data is 37.5Hz, the effective bandwidth is 17-60 Hz, and the absolute bandwidth is 43 Hz;
as shown in fig. 2, it is clear that the effective signal range of the original post-stack seismic data is the basic condition for determining the frequency division range.
S4, directly obtaining the approximate scale of the fracture development of the target horizon according to the regional adjustment data and the related geological data such as actual drilling, well logging and the like;
s5, respectively designing a theoretical forward modeling with the fault distance from small to large according to the development scale of the zonal adjustment fault, acquiring the speed value and the thickness value of the theoretical stratigraphic unit from the actual drilling data, and selecting an excitation receiving mode of self excitation and self collection in the forward modeling mode. Setting the forward modeling parameters according to the original amplitude data volume parameters;
as shown in FIG. 3, the size of the fracture development in the area is determined and forward simulation is performed, so that the seismic response characteristics of the fractures of different levels can be effectively determined.
S6, frequency division processing is carried out on the forward modeling earthquake result, section identifiability of different-grade fracture is observed under different frequency data, and optimal response frequency of different fracture grades is selected;
s7, carrying out frequency division processing of the preferred frequency in S6 on the data body within the effective signal range of the original amplitude data body to obtain data bodies of a plurality of frequency bands; distributing and acquiring narrow-band frequency dividers with main frequencies of 25Hz, 40Hz and 60Hz by using a frequency division technology of a rapid matching pursuit algorithm (MPD);
as shown in fig. 4, the frequency division processing is performed by using the theoretical forward modeling data, so that the fracture response frequencies of different orders can be better screened, the frequency division processing frequency of the original post-stack data is determined, and the data processing time is reduced.
S8, fracture attribute analysis and algorithm optimization; the large-scale fracture is identified by optimizing waveform difference attributes to the low and narrow frequency with the main frequency of 25Hz, the medium and narrow frequency with the main frequency of 40Hz adopts high eigenvalue coherence to identify the medium-scale fracture, and the high and narrow frequency with the main frequency of 60Hz adopts maximum likelihood attributes to identify the small-scale fracture.
S9, fracture bisection and three-dimensional space prediction effect analysis; the main stem fracture distribution area is small, the wave group dislocation section is obvious, the secondary fracture distance is small, the wave group dislocation section is often derived beside the main stem fracture, cracks are developed at local twisting, flexing and micro dislocation positions of the same phase axis, and the cracks are relatively widely developed in the whole area.
And S10, matching and matching the data related to fracture (fracture) such as actual drilling fracture and drilling encounter condition, fracture development condition, reservoir development condition, single well testing productivity and the like. The regional fracture development degree of the developing crack-hole type reservoir body is in a relative positive correlation with the reservoir and the single-well test yield.
As shown in fig. 5a, 5b, 5c, 5d, 5e, and 5f, the horizontal section, spatial distribution characteristics, and rationality of each stage of fracture were analyzed in conjunction with the background formed by the region fracture.
S11, eliminating the detection abnormality caused by the data, and completing the classification prediction of the fracture.
As shown in fig. 6, the multi-stage fracture is fused into a graph, and is matched with relevant data of actual drilling well in an analysis manner, so that abnormal fracture detection caused by data reasons is eliminated, and the graded prediction of fracture is completed.
It will be appreciated by those of ordinary skill in the art that the examples described herein are intended to assist the reader in understanding the manner in which the invention is practiced, and it is to be understood that the scope of the invention is not limited to such specifically recited statements and examples. Those skilled in the art can make various other specific changes and combinations based on the teachings of the present invention without departing from the spirit of the invention, and these changes and combinations are within the scope of the invention.
Claims (7)
1. A fracture grading prediction method based on a seismic frequency division technology is characterized by comprising the following steps:
s1, basic data collection: collecting data from the area of interest, the data comprising: zone adjustment data, a post-stack migration seismic body, a conventional logging curve, drilling and logging data, well deviation data, coring data, layering data, seismic interpretation horizon and single well testing productivity data;
s2, making a single-well synthetic record, and calibrating a research target layer and an upper layer and a lower layer;
s3, analyzing the frequency spectrum of the original post-stack amplitude data of the target layer to determine the effective signal range of the data of the target layer;
s4, directly obtaining the approximate scale of the target horizon fracture development according to the relevant geological data of actual drilling, well logging, head exposure and regional phase modulation;
s5, designing a theoretical forward model according to the fault development scale, and setting forward model parameters according to original amplitude data volume parameters;
s6, frequency division processing is carried out on the forward modeling earthquake result, section identifiability of different-grade fracture is observed under different frequency data, and optimal response frequency of different fracture grades is selected;
s7, carrying out frequency division processing of the preferred frequency in S6 on the data body within the effective signal range of the original amplitude data body to obtain data bodies of a plurality of frequency bands;
s8, analyzing each level of fracture property and optimizing an algorithm;
in S8, the low frequency body is used for the identification, prediction and calculation of large-scale fracture, the medium frequency body is used for the identification, prediction and calculation of medium-scale fracture, and the high frequency body is used for the identification, prediction and calculation of small-scale fracture or crack; calculating using waveform difference attributes for large-scale fracture, calculating using coherence and similarity attributes for medium-scale fracture, and calculating using curvature analysis, maximum likelihood attributes and ant tracking attributes for small-scale fracture or crack;
s9, analyzing the fracture bisection and three-dimensional space prediction effect of each level;
s10, fusing all levels of fractures into a graph, and matching and performing coincidence analysis on data related to fracture or fracture of actual drilling fracture and drilling conditions, fracture development conditions, reservoir development conditions and single well testing productivity;
s11, eliminating the detection abnormality caused by the data, and completing the classification prediction of the fracture.
2. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: and in the step S4, acquiring the approximate scale of the fracture development of the target horizon, wherein the fault acquisition parameters are with the fault distance as the first factor.
3. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: in S5, the forward modeling of the forward model is performed to obtain seismic forward data that is consistent with the original amplitude data parameters, and the method mainly includes: main frequency, bandwidth information.
4. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: in S7, the original amplitude data volume is divided by a frequency division technique based on fast matching pursuit.
5. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: in the S9 plano-section effect and three-dimensional space prediction analysis, a background is formed by combining basic structural characteristics and fracture of an actual research area, and the prediction misalignment caused by abnormity generated by frequency division processing is avoided by combining the comprehensive analysis of an original amplitude data volume.
6. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: and S10, matching and matching the data of actual drilling fracture and drilling encounter condition, fracture development condition, reservoir development condition and single well test productivity about fracture or fracture, and specifically analyzing the influence of the actual research area fracture cause, reservoir type and fracture on oil and gas transportation and gathering.
7. The fracture grading prediction method based on the seismic frequency division technology as claimed in claim 1, wherein: and S11, eliminating detection abnormality caused by data reasons, wherein the fracture detection abnormality often occurs aiming at the non-full coverage area of the three-dimensional data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010790656.3A CN111897008B (en) | 2020-08-07 | 2020-08-07 | Fracture grading prediction method based on seismic frequency division technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010790656.3A CN111897008B (en) | 2020-08-07 | 2020-08-07 | Fracture grading prediction method based on seismic frequency division technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111897008A CN111897008A (en) | 2020-11-06 |
CN111897008B true CN111897008B (en) | 2022-03-01 |
Family
ID=73246296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010790656.3A Active CN111897008B (en) | 2020-08-07 | 2020-08-07 | Fracture grading prediction method based on seismic frequency division technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111897008B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113296166A (en) * | 2021-05-27 | 2021-08-24 | 北京艾迪博科油气技术有限公司 | Method for constructing crack model |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1228385A1 (en) * | 1999-11-02 | 2002-08-07 | Geco AS | Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot |
CN102230974A (en) * | 2011-04-08 | 2011-11-02 | 中国石油大学(华东) | Three-dimensional high-precision bin fractionation processing and evaluation technology for seismic data |
CN103728659A (en) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | Method for improving underground karst detecting precision |
CN104316958A (en) * | 2014-10-20 | 2015-01-28 | 中国石油天然气集团公司 | Coherent processing method for identifying different scales of formation fractures |
GB201502027D0 (en) * | 2015-02-06 | 2015-03-25 | Foster Findlay Ass Ltd | A method for determining sedimentary facies using 3D seismic data |
CN105182424A (en) * | 2015-08-03 | 2015-12-23 | 中国石油天然气股份有限公司 | Method and device of reservoir porosity quantitative forecast based on patchy saturation model |
WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
CN105954797A (en) * | 2016-04-19 | 2016-09-21 | 中国石油天然气集团公司 | Fracture identification method and fracture identification device of seismic data |
CN108254783A (en) * | 2016-12-29 | 2018-07-06 | 中国石油化工股份有限公司 | A kind of poststack earthquake fluid recognition methods based on time frequency analysis |
CN109581485A (en) * | 2018-12-04 | 2019-04-05 | 成都捷科思石油天然气技术发展有限公司 | A method of carrying out automatic slit detection directly on pre-stack depth migration seismic data |
CN110568493A (en) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | Identification method of complex fault block basin hidden fault |
CN111290020A (en) * | 2020-03-25 | 2020-06-16 | 北京奥能恒业能源技术有限公司 | Fracture detection method and device based on structural filtering processing and frequency division attribute fusion |
CN111399056A (en) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | Method for predicting crack strength based on divided azimuth filtering |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2525586A (en) * | 2014-03-31 | 2015-11-04 | Foster Findlay Ass Ltd | Improved interpretation of seismic survey data using synthetic modelling |
CN105182415B (en) * | 2015-09-02 | 2018-02-13 | 中国石油化工股份有限公司 | A kind of fault recognizing method based on scaling down processing |
CN106291715B (en) * | 2016-09-24 | 2018-04-03 | 中国地质大学(北京) | A kind of low-grade fault law of development Forecasting Methodology based on fracture Self-similarity Theory |
GB2565526A (en) * | 2017-06-12 | 2019-02-20 | Foster Findlay Ass Ltd | A method for validating geological model data over corresponding original seismic data |
-
2020
- 2020-08-07 CN CN202010790656.3A patent/CN111897008B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1228385A1 (en) * | 1999-11-02 | 2002-08-07 | Geco AS | Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot |
CN102230974A (en) * | 2011-04-08 | 2011-11-02 | 中国石油大学(华东) | Three-dimensional high-precision bin fractionation processing and evaluation technology for seismic data |
CN103728659A (en) * | 2012-10-12 | 2014-04-16 | 中国石油化工股份有限公司 | Method for improving underground karst detecting precision |
WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
CN104316958A (en) * | 2014-10-20 | 2015-01-28 | 中国石油天然气集团公司 | Coherent processing method for identifying different scales of formation fractures |
GB201502027D0 (en) * | 2015-02-06 | 2015-03-25 | Foster Findlay Ass Ltd | A method for determining sedimentary facies using 3D seismic data |
CN105182424A (en) * | 2015-08-03 | 2015-12-23 | 中国石油天然气股份有限公司 | Method and device of reservoir porosity quantitative forecast based on patchy saturation model |
CN105954797A (en) * | 2016-04-19 | 2016-09-21 | 中国石油天然气集团公司 | Fracture identification method and fracture identification device of seismic data |
CN108254783A (en) * | 2016-12-29 | 2018-07-06 | 中国石油化工股份有限公司 | A kind of poststack earthquake fluid recognition methods based on time frequency analysis |
CN109581485A (en) * | 2018-12-04 | 2019-04-05 | 成都捷科思石油天然气技术发展有限公司 | A method of carrying out automatic slit detection directly on pre-stack depth migration seismic data |
CN110568493A (en) * | 2019-08-21 | 2019-12-13 | 中国石油化工股份有限公司 | Identification method of complex fault block basin hidden fault |
CN111290020A (en) * | 2020-03-25 | 2020-06-16 | 北京奥能恒业能源技术有限公司 | Fracture detection method and device based on structural filtering processing and frequency division attribute fusion |
CN111399056A (en) * | 2020-04-29 | 2020-07-10 | 西南石油大学 | Method for predicting crack strength based on divided azimuth filtering |
Non-Patent Citations (2)
Title |
---|
Porosity prediction from seismic inversion of a similarity attribute based on a pseudo-forward equation (PFE): a case study from the North Sea Basin, Netherlands;Saeed Mojeddifar 等;《Petroleum Science》;20150722;428–442 * |
文南油田低级序复杂断块精细刻画研究;马龙;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》;20180615;B019-9 * |
Also Published As
Publication number | Publication date |
---|---|
CN111897008A (en) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108931814B (en) | Multi-attribute fusion based bedrock fracture prediction method | |
US7113869B2 (en) | Method for prospect identification in asset evaluation | |
Hunt et al. | Quantitative estimate of fracture density variations in the Nordegg with azimuthal AVO and curvature: A case study | |
CN111596978A (en) | Web page display method, module and system for lithofacies classification by artificial intelligence | |
CN110703354B (en) | Lithology boundary identification method based on well-seismic combination | |
CN110658556A (en) | Seismic technology combination method for identifying and evaluating broken zone of sliding fracture of carbonate rock | |
CN106405644B (en) | Crack determines method and apparatus | |
RU2563323C1 (en) | Method of reconstruction of fine structure of geological object and forecast of its fluid saturation | |
CN114529110A (en) | Lithofacies inversion method and system based on deep neural network model | |
CN111897008B (en) | Fracture grading prediction method based on seismic frequency division technology | |
CN108957526B (en) | The method for obtaining crack fused data body | |
CN113484910A (en) | Tunnel advanced geological prediction method and system based on seismic interference method | |
CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
CN105242312A (en) | Method used for characterizing small and medium sized extension rift basin growth reverse fault | |
US20220120933A1 (en) | Method of detection of hydrocarbon horizontal slippage passages | |
CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
CN112946751A (en) | Method and system for predicting plane distribution of different reservoirs | |
CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
Lorentzen et al. | Mapping Cretaceous faults using a convolutional neural network-A field example from the Danish North Sea. | |
CN116027454A (en) | Method for imaging contour map of sandstone thickness and intelligent drawing device | |
CN113671566A (en) | Method for calculating crack parameters based on depth domain seismic data | |
Warren et al. | Final Phase 1 Report DE-EE0007698: A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography | |
CN112394392A (en) | Method and device for evaluating distribution condition of hydrocarbon source rock | |
CN117930380A (en) | Carbonate reservoir prediction method and system based on seismology | |
RU2758416C1 (en) | Method for reconstruction of fine structure of geological features and differentiation thereof into fractured and cavernous |
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 |