CN107765329A - Shale reservoir seismic identification - Google Patents
Shale reservoir seismic identification Download PDFInfo
- Publication number
- CN107765329A CN107765329A CN201610673190.2A CN201610673190A CN107765329A CN 107765329 A CN107765329 A CN 107765329A CN 201610673190 A CN201610673190 A CN 201610673190A CN 107765329 A CN107765329 A CN 107765329A
- Authority
- CN
- China
- Prior art keywords
- shale
- interval
- seismic
- favourable
- reservoir
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000011161 development Methods 0.000 claims abstract description 32
- 230000002349 favourable effect Effects 0.000 claims abstract description 29
- 208000010392 Bone Fractures Diseases 0.000 claims abstract description 15
- 206010017076 Fracture Diseases 0.000 claims abstract description 15
- 238000011156 evaluation Methods 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 12
- 238000010276 construction Methods 0.000 claims abstract description 10
- 230000001133 acceleration Effects 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims abstract description 6
- 230000005251 gamma ray Effects 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 102000000584 Calmodulin Human genes 0.000 claims description 3
- 108010041952 Calmodulin Proteins 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 2
- 230000008021 deposition Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000000034 method Methods 0.000 description 10
- 238000011160 research Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007306 turnover Effects 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
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses shale reservoir seismic identification, utilize structural configuration, depth, thickness and its distribution of seismic prospecting identification shale formation, complete region, which is developed, in the shale formation identified carries out prospect pit construction, logging operation is carried out to prospect pit, shale interval sedimentary micro is divided according to well log interpretation, favourable rammell sedimentary micro development interval is filtered out, determines its areal extent；Pre-stack seismic inversion is implemented to the favourable shale interval of the deposition filtered out again, filter out favourable oily interval, analyze the structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and variance attribute of such interval, filter out that curvature is big while geological data is irrelevant or discontinuous interval, as favourable fracture development interval, with reference to shale brittleness evaluation, high-quality shale reservoir is determined, realizes the purpose of the shale gas reservoir RESERVOIR RECOGNITION based on seismic technology.
Description
Technical field
The present invention relates to geological exploration and development field, relates in particular to shale reservoir seismic identification.
Background technology
Shale gas reservoir is to be present in organic matter and clay mineral surface in source bed with ADSORPTION STATE, or is present in free state
In shale matrix pores and crack, typical " being born from storage certainly " is shown as into the feature of Tibetan.Shale oil gas is as new life's in recent years
Energy extraction source, the exploitation of petroleum resources is carried out from the source bed for thinking to preserve oil gas in the past, in the world
With extensive prospect.And it is different from the pore channel in conventional sandstone reservoir, oil-gas migration is carried out using crack in rammell,
Crack can both provide aggregation space for shale gas, or the production of shale gas provides migration pathway.Shale gas has been me
The focus of state's unconventionaloil pool area research, Ministry of Land and Resources point out in national shale gas resource potential research and appraisal achievement,
China's shale gas exploration-development potential is huge, and industrial prospect is wide.However, at present exploration and development of the China to shale gas also in
Preferably with guide area experimental stage, technological means is also immature for resource assessment, Favorable Areas, and particularly seismic exploration technique lacks system
The research method of system, also rest on the research mode of conventional oil gas reservoir.
The content of the invention
It is an object of the invention to propose shale reservoir seismic identification, to solve the exploration of shale gas in the prior art
Development tool lacks the problem of systems approach, realizes the purpose of the shale gas reservoir RESERVOIR RECOGNITION based on seismic technology.
The present invention is achieved through the following technical solutions：
Shale reservoir seismic identification, including following technical step：（a）In a certain exploration and development block, carry out earthquake and survey
Visit, using seismic prospecting, surface geology data, identify structural configuration, depth, thickness and its distribution of shale formation；（b）Knowing
The shale formation not gone out develops complete region and carries out prospect pit construction, carries out logging operation to prospect pit, well logging project includes：Resistance
Rate well logging, acoustic travel time logging, nutural potential logging, gamma ray log, compensated neutron log, density log, VSP well loggings,
Well log interpretation is made, calmodulin binding domain CaM construction, divides shale interval sedimentary micro, filters out favourable rammell sedimentary micro development
Interval, determine its areal extent；（c）To step（b）In the favourable shale interval that filters out implement pre-stack seismic inversion, further
Favourable oily interval is filtered out, while shale brittleness evaluation is carried out to such interval；（d）To step（c）It is middle to be filtered out
Favourable oily interval, further Study of Seismic data, analyze the structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and variance of such interval
Attribute, filters out that curvature is big while geological data is irrelevant or discontinuous interval, as favourable fracture development interval；（e）
With reference to step（c）In shale brittleness evaluation and step（d）In fracture development interval, filter out stratum in fracture development interval
The high interval of brittleness index, is defined as high-quality shale reservoir.
For shale gas in the prior art exploration and development means lack systems approach the problem of, the present invention proposes shale
Seismic reservoir recognition methods.This method first in a certain exploration and development block, carry out seismic prospecting, using seismic prospecting,
Surface geology data, structural configuration, depth, thickness and its distribution of shale formation are identified, be to identify and chase after on seismic profile
Track rammell, demarcation rammell top bottom interface, analyzes its seismic wave individual features, determine the depth of rammell, thickness and point
Cloth scope, obtain rammell structural configuration, faults distribution and deposit thickness feature.Afterwards according to seismic profile, shale is chosen
Complete region is developed in rock stratum, carries out prospect pit construction, carries out the logging operation of full well section to the prospect pit, and well logging project is included but not
Be limited to resistivity logging, acoustic travel time logging, nutural potential logging, gamma ray log, compensated neutron log, density log,
VSP logs well, and makes well log interpretation.Wherein, VSP well loggings are vertical seismic profile log, and explosive source is located at earth's surface, in well
Different depth receives to seismic wave, the vertical change of geologic section near research well.VSP compared with surface seismic signal to noise ratio, point
Resolution is higher, and the kinematics and dynamic characteristic of ripple become apparent from.According to above-mentioned all well log interpretations, the structural region construction deposition back of the body
Scape is studied, individual well sedimentary section well log interpretation and seismic attributes analysis, and division shale interval sedimentary micro type and deposition are micro-
The changing rule of phase, it is determined that favourable rammell sedimentary micro development interval and its areal extent.
Afterwards, the shale interval favourable to the deposition determined implements pre-stack seismic inversion.Pre-stack seismic inversion is that stratum contains
One kind in gas detection means, the size of rammell air content determine whether shale gas exploitation is commercially valuable, so
Carry out Gas potential detection to be very important, in addition to pre-stack seismic inversion technology, use poststack wave impedance inversion, prestack bullet
The STUDIES ON SEISMIC WAVES methods such as property Impedance Inversion, spectrum analysis technique can also complete the detection of gas-bearing property.Wherein, made in the present invention
Prestack seismic inversion method, it is theoretical according to rock physicses theory and amplitude variation with Offset, by Zuo Yipuli
Hereby equation, the change to reflected amplitude are turned into least square fitting, until theoretical value is fitted well with observation, using amplitude with
Geophone offset variation relation curve, two parameters of intercept P and gradient G are calculated, then it is anti-using mathematical modeling by the two parameters
Elastic parameter needed for performance, and then lithology, fluid identification are carried out, it is standby to make stratum brittleness evaluation.Simultaneously as shale gas
The rammell of rich region, its bulk density can reduce, elastic wave velocity can reduce, therefore have to Gas potential detection parameter bright
Development rings, therefore can accurately be judged in the prestack inversion of seismic wave.
Hereafter, to the favourable oily interval filtered out using pre-stack seismic inversion, further Study of Seismic data, divide
Analyse the structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and variance attribute of such interval, filter out that curvature is big while geological data is irrelevant or
Discontinuous interval, as favourable fracture development interval.Curvature is a certain curve of reflection, the parameter of curved surface degree of crook.It is bent
The structural main curvature in face is bigger, easier generation crack, therefore structural main curvature can generally reflect the distribution in crack.
When curvature increases to stratum elastic limit, crack will be produced bending larger area, therefore in the both sides of fold axis, structure
Some higher curvature positions such as turnover position and the both sides of the plane of disruption are made, are all the development regions in crack.Acceleration Algorithm in Seismic Coherence Cube and side
The research object of poor attribute is similitude and discontinuity between adjacent seismic trace signal in three dimensional seismic data, and they can dash forward
Go out those irrelevant or discontinuous geological datas, the density in crack, orientation, intensity etc. in shale gas-bearing formation are identified with this, so as to
Select favourable fracture development interval.During the exploration and development of shale gas, intrinsic fracture is vital, and it can
Increase shale gas reservoir space, improve the extremely low matrix permeability of shale, the effect of artificial pressure is improved, therefore, in this step
Using structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and variance property analysis method, it is essential to select the shale interval of fracture development
's.
Finally, with reference to above-mentioned steps（c）In shale brittleness evaluation and step（d）In fracture development interval, filter out
Brittleness index of a relatively high interval in stratum in fracture development interval, is defined as high-quality shale reservoir.Shale reservoir fragility is bigger,
Sand fracturing is easier to be pressed off, and is easy to reservoir reconstruction, is advantageous to production development.Shale fragility earthquake prediction can pass through step
（c）In prestack inversion, the elastic parameter related to rock-mechanics property is obtained, so as to carry out shale brittleness evaluation.
Further, the step（b）In prospect pit work progress in, carry out rammell coring operation, draw taken rock
The modulus of elasticity of the heart, Poisson's ratio, anti-tensile, shearing resistance, compression strength, step is made with this（c）Described in shale brittleness evaluation.I.e.
Using the actual institute's coring in rammell, elastic parameter is obtained in a manner of experiment, it is more smart so as to be drawn to rammell brittleness index
True conclusion.The coring operation can carry out the shaft bottom in drilling process according to the rammell position that seismic prospecting is predicted and take
It the heart, can also be layered according to the geological logging in drilling process, sidewall coring is carried out after confirming to bore chance shale.
Further, the step（e）The high interval of middle stratum brittleness index, refer to the interval that brittleness index is more than 40.
The present invention compared with prior art, has the following advantages and advantages：
Shale reservoir seismic identification of the present invention, using seismic prospecting identify the structural configuration of shale formation, depth, thickness and
It is distributed, and developing complete region in the shale formation identified carries out prospect pit construction, logging operation is carried out to prospect pit, according to survey
Well explains division shale interval sedimentary micro, filters out favourable rammell sedimentary micro development interval, determines its areal extent；
Pre-stack seismic inversion is implemented to the favourable shale interval of the deposition filtered out again, filters out favourable oily interval, analyzes such
Structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and the variance attribute of interval, filter out curvature greatly while geological data is irrelevant or discontinuous
Interval, as favourable fracture development interval, with reference to shale brittleness evaluation, determine high-quality shale reservoir, realize with earthquake skill
Art is the purpose of the shale gas reservoir RESERVOIR RECOGNITION of main body.
Brief description of the drawings
Accompanying drawing described herein is used for providing further understanding the embodiment of the present invention, forms one of the application
Point, do not form the restriction to the embodiment of the present invention.In the accompanying drawings：
Fig. 1 is the schematic process flow diagram of a specific embodiment of the invention.
Embodiment
For the object, technical solutions and advantages of the present invention are more clearly understood, with reference to embodiment and accompanying drawing, to this
Invention is described in further detail, and exemplary embodiment of the invention and its explanation are only used for explaining the present invention, do not make
For limitation of the invention.
Embodiment 1：
Shale reservoir seismic identification as shown in Figure 1, including following technical step：（a）In a certain exploration and development block,
Carry out seismic prospecting, using seismic prospecting, surface geology data, identify the structural configuration of shale formation, depth, thickness and its point
Cloth；（b）Complete region is developed in the shale formation identified and carries out prospect pit construction, and logging operation is carried out to prospect pit, item of logging well
Mesh includes：Resistivity logging, acoustic travel time logging, nutural potential logging, gamma ray log, compensated neutron log, density are surveyed
Well, VSP well loggings, well log interpretation is made, calmodulin binding domain CaM construction, shale interval sedimentary micro is divided, filters out favourable rammell
Sedimentary micro develops interval, determines its areal extent；（c）To step（b）In the favourable shale interval that filters out with implementing prestack
Inverting is shaken, further filters out favourable oily interval, while shale brittleness evaluation is carried out to such interval；（d）To step
（c）Middle filtered out favourable oily interval, further Study of Seismic data, analyze such interval structural main curvature,
Seismic phase stem body and variance attribute, filter out that curvature is big while geological data is irrelevant or discontinuous interval, are split as favourable
Seam development interval；（e）With reference to step（c）In shale brittleness evaluation and step（d）In fracture development interval, filter out crack
The interval that stratum brittleness index in interval is more than 40 is developed, is defined as high-quality shale reservoir.Realized with this based on seismic technology
The purpose of the shale gas reservoir RESERVOIR RECOGNITION of body.
Above-described embodiment, the purpose of the present invention, technical scheme and beneficial effect are carried out further
Describe in detail, should be understood that the embodiment that the foregoing is only the present invention, be not intended to limit the present invention
Protection domain, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., all should include
Within protection scope of the present invention.
Claims (3)
1. shale reservoir seismic identification, it is characterised in that including following technical step：
（a）In a certain exploration and development block, seismic prospecting is carried out, using seismic prospecting, surface geology data, identifies shale rock
Structural configuration, depth, thickness and its distribution of layer；
（b）Complete region is developed in the shale formation identified and carries out prospect pit construction, and logging operation is carried out to prospect pit, item of logging well
Mesh includes：Resistivity logging, acoustic travel time logging, nutural potential logging, gamma ray log, compensated neutron log, density are surveyed
Well, VSP well loggings, well log interpretation is made, calmodulin binding domain CaM construction, shale interval sedimentary micro is divided, filters out favourable rammell
Sedimentary micro develops interval, determines its areal extent；
（c）To step（b）In the favourable shale interval that filters out implement pre-stack seismic inversion, further filter out favourable oil-containing
Gas-bearing formation section, while shale brittleness evaluation is carried out to such interval；
（d）To step（c）Middle filtered out favourable oily interval, further Study of Seismic data, analyzes such interval
Structural main curvature, Acceleration Algorithm in Seismic Coherence Cube and variance attribute, filter out that curvature is big while geological data is irrelevant or discontinuous layer
Section, as favourable fracture development interval；
（e）With reference to step（c）In shale brittleness evaluation and step（d）In fracture development interval, filter out fracture development layer
The high interval of stratum brittleness index, is defined as high-quality shale reservoir in section.
2. shale reservoir seismic identification according to claim 1, it is characterised in that：The step（b）In prospect pit
In work progress, rammell coring operation is carried out, draws modulus of elasticity, Poisson's ratio, anti-tensile, shearing resistance, the pressure resistance of institute's coring
Degree, step is made with this（c）Described in shale brittleness evaluation.
3. shale reservoir seismic identification according to claim 1, it is characterised in that：The step（e）Middle stratum is crisp
The high interval of sex index, refer to the interval that brittleness index is more than 40.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610673190.2A CN107765329A (en) | 2016-08-17 | 2016-08-17 | Shale reservoir seismic identification |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610673190.2A CN107765329A (en) | 2016-08-17 | 2016-08-17 | Shale reservoir seismic identification |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107765329A true CN107765329A (en) | 2018-03-06 |
Family
ID=61259568
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610673190.2A Withdrawn CN107765329A (en) | 2016-08-17 | 2016-08-17 | Shale reservoir seismic identification |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107765329A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112034514A (en) * | 2020-07-15 | 2020-12-04 | 成都理工大学 | Shale oil-gas system based on seismic sedimentology theory and identification method thereof |
-
2016
- 2016-08-17 CN CN201610673190.2A patent/CN107765329A/en not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112034514A (en) * | 2020-07-15 | 2020-12-04 | 成都理工大学 | Shale oil-gas system based on seismic sedimentology theory and identification method thereof |
| CN112034514B (en) * | 2020-07-15 | 2021-09-14 | 成都理工大学 | Shale oil-gas system based on seismic sedimentology theory and identification method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Grammer et al. | Integration of outcrop and modern analogs in reservoir modeling: Overview with examples from the Bahamas | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| CN101158724B (en) | Reservoir thickness prediction method based on dipolar wavelet | |
| WO2019062655A1 (en) | Method and device for determining thin interlayer | |
| CN102466815B (en) | Triassic clastic rock petroleum reservoir identification method | |
| CN102121995A (en) | Seismic inversion reservoir stratum forecasting method for overthrust-containing complex structure | |
| CN105116449B (en) | A kind of recognition methods of weak reflection reservoir | |
| CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN105938203A (en) | Detection method and apparatus of reservoir features | |
| CN105425299A (en) | Method and device for determining formation fracture distribution | |
| Laudon et al. | Machine learning applied to 3D seismic data from the Denver-Julesburg basin improves stratigraphic resolution in the Niobrara | |
| Liu et al. | Seismic characterization of fault and fractures in deep buried carbonate reservoirs using CNN-LSTM based deep neural networks | |
| Chashkov et al. | Use of the cluster analysis and artificial neural network technology for log data interpretation | |
| CN111983683A (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| Irofti et al. | Characterization of a Tight Gas Reservoir Using the Integration of Electrofacies and Fracture Aperture, Ahnet, Algeria | |
| CN107765329A (en) | Shale reservoir seismic identification | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| Patel et al. | Seismic illumination of small-throw seismogenic faults, Anadarko Basin, Oklahoma | |
| 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 | |
| Ningkai et al. | Stepped and detailed seismic prediction of shallow-thin reservoirs in Chunfeng oilfield of Junggar Basin, NW China | |
| Ross et al. | Tomographic fracture imaging: Examples of induced fracture and reservoir-scale observations during wellbore stimulations, Niobrara and Bakken plays, USA | |
| Dowdell et al. | An Integrated Study of a Mississippian Tripolitic Chert Reservoir-Osage County, Oklahoma, USA | |
| CN107765304A (en) | Shale reservoir seismic identification | |
| Jafri et al. | Improved reservoir characterisation using fuzzy logic platform: an integrated petrophysical, seismic structural and poststack inversion study |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180306 |
|
| WW01 | Invention patent application withdrawn after publication |