CN108761535A - A kind of recognition methods of intrusive rock distribution - Google Patents
A kind of recognition methods of intrusive rock distribution Download PDFInfo
- Publication number
- CN108761535A CN108761535A CN201810895651.XA CN201810895651A CN108761535A CN 108761535 A CN108761535 A CN 108761535A CN 201810895651 A CN201810895651 A CN 201810895651A CN 108761535 A CN108761535 A CN 108761535A
- Authority
- CN
- China
- Prior art keywords
- intrusive rock
- rock
- intrusive
- wave field
- thickness
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 238000000926 separation method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 claims description 2
- 230000009545 invasion Effects 0.000 abstract description 6
- 239000010410 layer Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 230000009466 transformation 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of recognition methods of intrusive rock distribution, steps:Wave field forward modeling analysis determines intrusive rock reflectance signature;Intrusive rock wave field is detached using wavelet reconfiguration technique;Sensitive Attributes analysis determines the plane distribution of intrusive rock;Intrusive rock thickness quantitative forecast.The present invention can be complicated in igneous invasion rock spatial distribution, in the case of seismic data frequency is lower, accurately identifies and portrays igneous invasion rock spatial distribution;By establishing the quantitative relationship of Sensitive Attributes and intrusive rock thickness, the purpose of identification thin layer intrusive rock thickness is realized.
Description
Technical field
The present invention relates to a kind of intrusion magmatic rock spatial distribution recognition methods of petroleum exploration field, especially with regard to one kind
The recognition methods of intrusive rock distribution in the case of and igneous invasion rock spatial distribution complexity relatively low in seismic data frequency.
Background technology
Conventional igneous invasion rock depicting method can successfully portray the plane distributing scope of intrusive rock, and may be implemented to invade
Enter the thickness prediction of rock.But when intrusive rock distribution is complex, physical properties of rock is more close with country rock, meanwhile, seismic data
In the case of frequency is lower, existing method identifies that the multi-solution of intrusive rock increases, and the thickness prediction of intrusive rock is also difficult to realize.Separately
On the one hand, can be realized by the simple combination of multi-method reduces the purpose of prediction multi-solution, but often improvement effect is limited,
And effective identification of thin layer intrusive rock space distribution rule can not be completed.These problems make conventional method carry out under complex situations
Intrusive rock space identity difficulty is larger.
Invention content
In view of the above-mentioned problems, the object of the present invention is to provide a kind of recognition methods of intrusive rock distribution, it can be in magma
Intrusive rock spatial distribution is complicated, in the case of seismic data frequency is lower, accurately identifies and portrays igneous invasion rock spatial distribution.
To achieve the above object, the present invention takes following technical scheme:A kind of recognition methods of intrusive rock distribution comprising
Following steps:1) wave field forward modeling analysis determines intrusive rock reflectance signature;2) wavelet reconfiguration technique is applied to carry out intrusive rock wave field
Separation;3) Sensitive Attributes analysis determines the plane distribution of intrusive rock;4) intrusive rock thickness quantitative forecast.
Further, it in the step 1), is fluctuated by the forward modeling analysis and one-dimensional individual well that carry out two-dimentional typical earthquake section
The forward modeling of equation is as a result, determine the reflected wave field feature of known location intrusive rock.
Further, in the step 2), according to wave field forward modeling as a result, right in particular dimensions using wavelet reconfiguration technique
The seismic wave field of intrusive rock is detached, and ground coincideing with wave field forward modeling result, can reflecting intrusive rock feature is finally obtained
Shake data volume.
Further, the separation process is:In conjunction with intrusive rock reflected wave field feature, different scale seismic data wave is searched for
, it will reflect that one or more scale seismic wavefield datas of intrusive rock wave field information preferably come out, obtain to reflect
The seismic data cube of thin layer intrusive rock wave field characteristics in reservoir, and then these seismic data cubes are reconstructed.
Further, acquired in step 2) according to the intrusive rock reflectance signature determined in step 1) in the step 3)
Data volume on carry out Sensitive Attributes analysis, so that it is determined that going out the plane distribution of intrusive rock.
Further, the Sensitive Attributes analytic process is:By the forward modeling analysis result of step 1), intrusive rock is determined
The difference of wave field characteristics and country rock reservoir, and tentative calculation is carried out on the obtained seismic data cube of step 2), knot is met with practical bore
Fruit is hard constraint, and the geological knowledge of intrusive rock distribution is foundation, preferably goes out the seismic wave field category that can most meet above-mentioned two conditions
Property, the plane distribution of final preferably negative amplitude attribute prediction intrusive rock.
Further, in the step 4), according to the one-dimensional individual well Wave equation forward modeling of forward modeling in step 1) as a result, based on step
It is rapid 2) based on the Sensitive Attributes selected by the obtained seismic data of wavelet reconfiguration technique, it is sensitive with this to establish intrusive rock thickness
The quantitative relationship of attribute, so that it is determined that the thickness distribution of intrusive rock.
Further, determine that the process of the intrusive rock thickness distribution is:The forward modeling provided according to well data is as a result, determination is invaded
The presence for entering rock affects greatly reflected wave field amplitude;With this, the wavelet based on step 2) reconstructs data decimation sensitivity category
Property, by the method for linear regression, the intrusive rock thickness met and the well point position amplitude attribute prediction knot are bored by part drilling well
Fruit, quantitative analysis different-thickness intrusive rock amplitude attribute response characteristic, and determine in peak swing attribute and intrusive rock thickness
In quantitative relationship, in the intrusive rock thickness range of 0-100m, intrusive rock thickness has linear well close with peak swing value
System, related coefficient can be reached for 0.91;On this basis, the thickness for calculating the intrusive rock that step 3) determines, to obtain
The spatial distribution of intrusive rock.
The invention adopts the above technical scheme, which has the following advantages:1, the present invention can be in igneous invasion rock point
In the case of cloth complexity, pass through the spatial distribution of multi-method comprehensive constraint and effective information extraction identification intrusive rock.2, of the invention
By establishing the quantitative relationship of Sensitive Attributes and intrusive rock thickness, the purpose of identification thin layer intrusive rock thickness is realized.
Description of the drawings
Fig. 1 is the overall flow schematic diagram of the present invention;
Fig. 2 a are intrusive rock Seismic reflection character forward simulation model schematics;
Fig. 2 b are intrusive rock Seismic reflection character forward modeling diagrammatic cross-sections;
Fig. 3 is individual well Wave equation forward modeling schematic diagram;
Fig. 4 a are seismic profile wavelet reconstruct wave field separation result schematic diagrams;
Fig. 4 b are reconstruct section wavelet reconstruct wave field separation result schematic diagrams;
Fig. 5 is intrusive rock thickness and peak swing cross plot.
Specific implementation mode
The present invention is described in detail below with reference to the accompanying drawings and embodiments.
As shown in Figure 1, the present invention provides a kind of recognition methods of intrusive rock distribution comprising following steps:
1) wave field forward modeling analysis determines intrusive rock reflectance signature;
Forward modeling by carrying out two-dimentional typical earthquake section is analyzed and the forward modeling of one-dimensional individual well wave equation is as a result, can be true
Determine the reflected wave field feature of known location intrusive rock.
In the present embodiment, it is illustrated by taking the complex situations at thin layer intrusive rock surface layer top interface as an example.Pass through two dimension
Forward modelling result (as shown in Fig. 2 a, Fig. 2 b) is it can be found that intrusive rock reflectance signature is integrally buried in the strong amplitude in reservoir top
In secondary lobe, without apparent interface feature.But simultaneously it has also been found that intrusive rock causes amplitude to enhance about 40%, and stratum has thickening to become
Gesture.In addition, by all-wave earthquake analogue technique, the sunykatuib analysis of one-dimensional wave equation has been carried out to fixed well.Analog result
Show (as shown in Figure 3), the presence of intrusive rock makes seismic amplitude be remarkably reinforced.But percentage is enhanced for positive amplitude and negative amplitude
Than and differ.The presence of intrusive rock makes positive amplitude enhancing about 30%, and negative amplitude then enhances 70%.
2) intrusive rock wave field is detached using wavelet reconfiguration technique.
According to wave field forward modeling as a result, being carried out to the seismic wave field of intrusive rock in particular dimensions using wavelet reconfiguration technique
Separation, finally obtain with wave field forward modeling result coincide, can reflect the seismic data cube of intrusive rock feature.
In the present embodiment, separation process is said by taking the complex situations at thin layer intrusive rock nearly reservoir top interface as an example
It is bright:
Using small wave converting method, different scale small echo may separate out more set seismic data cubes.In conjunction with invading in step 1)
Enter rock reflected wave field feature, searches for different scale seismic data wave field, will reflect one or more of intrusive rock wave field information
A scale seismic wavefield data preferably comes out, and obtains to reflect the seismic data cube of thin layer intrusive rock wave field characteristics in reservoir,
And then these seismic data cubes are reconstructed.The reflected wave field of intrusive rock and stratum are successfully pushed up interface by wavelet reconfiguration technique
Wave field detached, the spatial position of plutone meets preferable (as shown in Fig. 4 a, Fig. 4 b) with fixed well.
3) Sensitive Attributes analysis determines the plane distribution of intrusive rock
According to the intrusive rock reflectance signature determined in step 1), sensitive category is carried out on obtained data volume in step 2)
Property analysis, so that it is determined that going out the plane distribution of intrusive rock.
In the present embodiment, Sensitive Attributes are analyzed by taking the complex situations at thin layer intrusive rock nearly reservoir top interface as an example and is carried out
Explanation:
By the forward modeling analysis result of step 1), the wave field characteristics for determining intrusive rock and the difference of country rock reservoir are (general
Carry out but be not limited to the analysis of amplitude, frequency, phase and structure generic attribute), and on the obtained seismic data cube of step 2)
Carry out tentative calculation, result met as hard constraint using practical bore, the geological knowledge of intrusive rock distribution is foundation, preferably go out can most meet it is above-mentioned
The seismic wave field attribute of two conditions, the plane distribution of final preferably negative amplitude attribute prediction intrusive rock.
4) intrusive rock thickness quantitative forecast;
According to the one-dimensional individual well Wave equation forward modeling of forward modeling in step 1) as a result, reconstructing skill based on wavelet based on step 2)
Sensitive Attributes selected by the obtained seismic data of art, establish the quantitative relationship of intrusive rock thickness and this Sensitive Attributes, to
Determine the thickness distribution of intrusive rock.
In the present embodiment, it is illustrated by taking the complex situations at thin layer intrusive rock nearly reservoir top interface as an example.According to well number
According to the forward modeling provided as a result, determining that the presence of intrusive rock affects greatly reflected wave field amplitude.With this, based on step 2)
Wavelet reconstructs data decimation Sensitive Attributes, by the method for linear regression, bores the intrusive rock thickness met by part drilling well and is somebody's turn to do
Well point position amplitude attribute prediction result, quantitative analysis different-thickness intrusive rock amplitude attribute response characteristic, and determine in maximum
In the quantitative relationship of amplitude attribute and intrusive rock thickness, in the intrusive rock thickness range of 0-100m, intrusive rock thickness and maximum
There is amplitude good linear relationship, related coefficient can be reached for 0.91 (as shown in Figure 5).It on this basis, can be with
The thickness for calculating the intrusive rock that step 3) determines, to obtain the spatial distribution of intrusive rock.
The various embodiments described above are merely to illustrate the present invention, and structure and size, installation position and the shape of each component are all can be with
It is varied from, based on the technical solution of the present invention, all improvement that individual part is carried out according to the principle of the invention and waits
With transformation, should not exclude except protection scope of the present invention.
Claims (8)
1. a kind of recognition methods of intrusive rock distribution, which is characterized in that include the following steps:
1) wave field forward modeling analysis determines intrusive rock reflectance signature;
2) intrusive rock wave field is detached using wavelet reconfiguration technique;
3) Sensitive Attributes analysis determines the plane distribution of intrusive rock;
4) intrusive rock thickness quantitative forecast.
2. method as described in claim 1, it is characterised in that:In the step 1), by carrying out two-dimentional typical earthquake section
Forward modeling is analyzed and the forward modeling of one-dimensional individual well wave equation is as a result, determine the reflected wave field feature of known location intrusive rock.
3. method as described in claim 1, it is characterised in that:In the step 2), according to wave field forward modeling as a result, using wavelet
Reconfiguration technique detaches the seismic wave field of intrusive rock in particular dimensions, finally obtain with wave field forward modeling result coincide,
It can reflect the seismic data cube of intrusive rock feature.
4. method as claimed in claim 3, it is characterised in that:The separation process is:In conjunction with intrusive rock reflected wave field feature, search
Rope different scale seismic data wave field will reflect that one or more scale seismic wavefield datas of intrusive rock wave field information are excellent
It elects, obtains to reflect the seismic data cube of thin layer intrusive rock wave field characteristics in reservoir, and then by these seismic data cubes
It is reconstructed.
5. method as described in claim 1, it is characterised in that:It is anti-according to the intrusive rock determined in step 1) in the step 3)
Feature is penetrated, Sensitive Attributes analysis is carried out on obtained data volume in step 2), so that it is determined that going out the plane distribution of intrusive rock.
6. method as claimed in claim 5, it is characterised in that:The Sensitive Attributes analytic process is:Pass through the forward modeling of step 1)
Analysis result determines the difference of the wave field characteristics and country rock reservoir of intrusive rock, and in the obtained seismic data cube of step 2)
Upper carry out tentative calculation meets result as hard constraint using practical bore, and the geological knowledge of intrusive rock distribution is foundation, and preferably going out can most meet
The seismic wave field attribute of two conditions is stated, the plane distribution of final preferably negative amplitude attribute prediction intrusive rock.
7. method as described in claim 1, it is characterised in that:In the step 4), according to the one-dimensional individual well of forward modeling in step 1)
Wave equation forward modeling based on the sensitivity selected by the obtained seismic data of wavelet reconfiguration technique as a result, being belonged to based on step 2)
Property, the quantitative relationship of intrusive rock thickness and this Sensitive Attributes is established, so that it is determined that the thickness distribution of intrusive rock.
8. method as claimed in claim 7, it is characterised in that:Determine that the process of the intrusive rock thickness distribution is:According to well number
According to the forward modeling provided as a result, determining that the presence of intrusive rock affects greatly reflected wave field amplitude;With this, based on step 2)
Wavelet reconstructs data decimation Sensitive Attributes, by the method for linear regression, bores the intrusive rock thickness met by part drilling well and is somebody's turn to do
Well point position amplitude attribute prediction result, quantitative analysis different-thickness intrusive rock amplitude attribute response characteristic, and determine in maximum
In the quantitative relationship of amplitude attribute and intrusive rock thickness, in the intrusive rock thickness range of 0-100m, intrusive rock thickness and maximum
There is amplitude good linear relationship, related coefficient can be reached for 0.91;On this basis, calculate what step 3) determined
The thickness of intrusive rock, to obtain the spatial distribution of intrusive rock.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810895651.XA CN108761535B (en) | 2018-08-08 | 2018-08-08 | Method for identifying distribution of invaded rocks |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810895651.XA CN108761535B (en) | 2018-08-08 | 2018-08-08 | Method for identifying distribution of invaded rocks |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761535A true CN108761535A (en) | 2018-11-06 |
CN108761535B CN108761535B (en) | 2020-04-07 |
Family
ID=63969302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810895651.XA Active CN108761535B (en) | 2018-08-08 | 2018-08-08 | Method for identifying distribution of invaded rocks |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761535B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002369933A (en) * | 2001-06-13 | 2002-12-24 | Takao:Kk | Pachinko game machine |
CN1595107A (en) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | Multiple optical spectrum alteration zoning method based on wavelet packet variation |
US20080232193A1 (en) * | 2007-03-20 | 2008-09-25 | Geocyber Solutions, Inc. | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection |
CN104280770A (en) * | 2014-09-28 | 2015-01-14 | 中国石油大港油田勘探开发研究院 | Prediction method of compact transition rock reservoir stratum |
CN108226999A (en) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | The processing method of the small scale fracture hole body information of carbonate rock |
-
2018
- 2018-08-08 CN CN201810895651.XA patent/CN108761535B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002369933A (en) * | 2001-06-13 | 2002-12-24 | Takao:Kk | Pachinko game machine |
CN1595107A (en) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | Multiple optical spectrum alteration zoning method based on wavelet packet variation |
US20080232193A1 (en) * | 2007-03-20 | 2008-09-25 | Geocyber Solutions, Inc. | Methods for Noise Removal and/or Attenuation from Seismic Data by Wavelet Selection |
CN104280770A (en) * | 2014-09-28 | 2015-01-14 | 中国石油大港油田勘探开发研究院 | Prediction method of compact transition rock reservoir stratum |
CN108226999A (en) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | The processing method of the small scale fracture hole body information of carbonate rock |
Also Published As
Publication number | Publication date |
---|---|
CN108761535B (en) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8892388B2 (en) | Box counting enhanced modeling | |
CN105445800B (en) | A kind of recognition methods of the different lithological pool in thick sand bodies top | |
CN102736107B (en) | Energy constraint heterogeneous reservoir thickness identification system | |
CN105093313B (en) | A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device | |
CN104237945B (en) | A kind of seismic data self adaptation high resolution processing method | |
CN111580181B (en) | Water guide collapse column identification method based on multi-field multi-feature information fusion | |
CN106597543B (en) | A kind of Sedimentary Facies division methods | |
CN109541685B (en) | River channel sand body identification method | |
CN109001801B (en) | Fault variable-scale identification method based on multiple iteration ant colony algorithm | |
CN106291682B (en) | A kind of poststack sound impedance inversion method based on base method for tracing | |
CN107678072B (en) | Based on magnetic force, earthquake, the united igneous reservoirs prediction technique of drilling well | |
CN105093306A (en) | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration | |
CN109763814B (en) | Stratum matching visual analysis method based on multi-dimensional logging data | |
CN111650643B (en) | Method and system for extracting training set of seismic fault machine learning model | |
Pepper et al. | PSAdvances in Seismic Fault Interpretation Automation | |
CN110261905A (en) | Complex value based on pitch angle control is concerned with microfault recognition methods | |
CN111257933B (en) | Novel method for predicting oil and gas reservoir based on low-frequency shadow phenomenon | |
CN117075221A (en) | Reservoir fluid identification method based on combination of seismic elastic properties, resistivity and polarization rate | |
CN108680950B (en) | A kind of desert seismic signal method for detecting position based on Self-adaptive Block Matching | |
CN112711068A (en) | Method and device for predicting effective oil and gas reservoir in sandstone | |
CN108761535A (en) | A kind of recognition methods of intrusive rock distribution | |
CN116027454A (en) | Method for imaging contour map of sandstone thickness and intelligent drawing device | |
Cheng et al. | Multiscale fracture prediction technique via deep learning, seismic gradient disorder, and aberrance: Applied to tight sandstone reservoirs in the Hutubi block, southern Junggar Basin | |
CN109581489B (en) | Nested seismic facies extraction method and system | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked 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 |