CN109143326A - The method and apparatus for identifying coal seam top bottom interface - Google Patents
The method and apparatus for identifying coal seam top bottom interface Download PDFInfo
- Publication number
- CN109143326A CN109143326A CN201810792633.9A CN201810792633A CN109143326A CN 109143326 A CN109143326 A CN 109143326A CN 201810792633 A CN201810792633 A CN 201810792633A CN 109143326 A CN109143326 A CN 109143326A
- Authority
- CN
- China
- Prior art keywords
- curve
- coal seam
- acoustic
- weighted
- identification
- 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.)
- Pending
Links
- 239000003245 coal Substances 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000035945 sensitivity Effects 0.000 claims abstract description 46
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 15
- 238000005516 engineering process Methods 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 14
- 239000002131 composite material Substances 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000005251 gamma ray Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000427 antigen Substances 0.000 description 1
- 102000036639 antigens Human genes 0.000 description 1
- 108091007433 antigens Proteins 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000003317 industrial substance Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/24—Recording seismic data
-
- 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
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 application is the method and apparatus about a kind of identification coal seam top bottom interface.The described method includes: obtaining multiple for reflecting the sensitivity curve of different lithology；Weighting pseudo-acoustic curve is generated according to the sensitivity curve and sound wave curve, the weighting pseudo-acoustic curve is used to reflect the nuance of formation velocity, the variation of wave impedance and lithology；Global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify the top bottom interface in coal seam.The application improves the identification accuracy of the top bottom interface in coal seam, reduces the interference of limitation and multi-solution brought by monotechnics.
Description
Technical field
The application is about fossil fuel technologies field, especially with respect to the method and dress of a kind of identification coal seam top bottom interface
It sets.
Background technique
Coal seam refers to changed as plant remains made of the flammable mineral products of layered solid.The distribution in coal seam, thickness and bury depth
Degree is the important evidence of determining coalfield economic value and construction of coal mine.
Coal bed gas belongs to Unconventional gas, hides oneself in coal seam, being high-grade energy and industrial chemicals, to predict coal seam
Gas needs first accurately to depict by the distribution in coal seam；In addition Seam Roof And Floor rich water, corresponding coal mining is extremely important safely
, to predict its roof and floor situation, coal seam distribution, which first has to accurately depict, to be come.
Research in terms of having carried out geology, geophysics for coal seam resource at present, achieves a collection of technological achievement.This
A little achievements mainly include: Seismic reflection character research, post-stack inversion etc..These means can predict point in coal seam to a certain extent
Cloth, but every kind of method has multi-solution, and inaccurate.
Summary of the invention
In order to overcome the problems, such as present in the relevant technologies, this application provides it is a kind of identification coal seam top bottom interface method and
Device.The technical solution is as follows:
According to the embodiment of the present application in a first aspect, providing a kind of method of identification coal seam top bottom interface, the method packet
It includes:
It obtains multiple for reflecting the sensitivity curve of different lithology；
Weighting pseudo-acoustic curve is generated according to the sensitivity curve and sound wave curve, the weighting pseudo-acoustic curve is for anti-
Reflect the nuance of formation velocity, the variation of wave impedance and lithology；
Global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify the top bottom interface in coal seam.
It is further, described that weighting pseudo-acoustic curve is generated according to the sensitivity curve and sound wave curve, comprising:
It is weighted based on the sensitivity curve and coalbed coring relationship, the weighting for obtaining coal seam for identification is bent
Line；
The low-frequency information of the high-frequency information of the weighted curve and sound wave curve is merged, the weighting onomatopoeia is obtained
Wave profile.
Further, the method also includes:
Using Wavelet Decomposition Technology, the high-frequency information of the weighted curve in coal seam for identification is extracted；And/or
Using Wavelet Decomposition Technology, the low-frequency information of the sound wave curve is extracted.
Further, described to be weighted based on the sensitivity curve and coalbed coring relationship, it obtains for knowing
The weighted curve in other coal seam, comprising:
Nonlinear weight coefficient is calculated according to sensitivity curve and the coal geology feature goodness of fit；
Using the nonlinear weight coefficient, weighted calculation obtains the weighted curve in the coal seam for identification.
Further, described that global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify coal
The top bottom interface of layer, comprising:
Obtain basic data, the basic data includes at least one of: poststack data, modeling data, layer position data,
Wavelet data；
According to the basic data, objective function is established, the objective function is used at least one of: actual seismic note
Record, composite traces, prior-constrained, transverse gradients constraint；
Using the fast inversion algorithm of global optimizing, amendment is iterated to initial geological model, is obtained high-resolution
Sound impedance/quasi-acoustic inversion model；
According to the sound impedance/quasi-acoustic inversion model and the objective function, inversion result is determined；
Based on the inversion result, the top bottom interface in the coal seam is identified.
According to the second aspect of the embodiment of the present application, a kind of device of identification coal seam top bottom interface, described device packet are provided
It includes:
Acquiring unit, it is multiple for reflecting the sensitivity curve of different lithology for obtaining；
Generation unit, for generating weighting pseudo-acoustic curve according to the sensitivity curve and sound wave curve, the weighting is quasi-
Sound wave curve is used to reflect the nuance of formation velocity, the variation of wave impedance and lithology；
Inverting unit, for carrying out global optimizing broad-band constrained inversion based on the weighting pseudo-acoustic curve, to identify coal
The top bottom interface of layer.
Further, the generation unit includes:
Computing module is obtained for being weighted based on the sensitivity curve and coalbed coring relationship for knowing
The weighted curve in other coal seam；
Fusion Module is obtained for merging the low-frequency information of the high-frequency information of the weighted curve and sound wave curve
To the weighting pseudo-acoustic curve.
Further, described device further include:
First extraction module extracts the height of the weighted curve in coal seam for identification for utilizing Wavelet Decomposition Technology
Frequency information；And/or
Second extraction module extracts the low-frequency information of the sound wave curve for utilizing Wavelet Decomposition Technology.
Further, the computing module includes:
First computational submodule, for calculating nonlinear weight system according to sensitivity curve and the coal geology feature goodness of fit
Number；
Second computational submodule, for utilizing the nonlinear weight coefficient, weighted calculation obtains the coal for identification
The weighted curve of layer.
Further, the inverting unit includes:
Basic data obtains module, and for obtaining basic data, the basic data includes at least one of: poststack number
According to, modeling data, layer position data, wavelet data；
Objective function establishes module, for according to the basic data, establishing objective function, the objective function be used for
It is at least one lower: real seismic record, composite traces, prior-constrained, transverse gradients constraint；
Iterated revision module is iterated initial geological model and repairs for the fast inversion algorithm using global optimizing
Just, high-resolution sound impedance/quasi-acoustic inversion model is obtained；
Inversion result determining module is used for according to the sound impedance/quasi-acoustic inversion model and the objective function, really
Determine inversion result；
Bottom interface identification module is pushed up, for being based on the inversion result, identifies the top bottom interface in the coal seam.
The technical solution that embodiments herein provides can include the following benefits:
It is multiple for reflecting the sensitivity curve of different lithology by obtaining in the embodiment of the present application, according to described sensitive bent
Line and sound wave curve generate weighting pseudo-acoustic curve, and the weighting pseudo-acoustic curve is used to reflect the change of formation velocity, wave impedance
Change the nuance with lithology, and global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify coal
The top bottom interface of layer improves the identification accuracy of the top bottom interface in coal seam, reduces limitation and more brought by monotechnics
The interference of solution property.
It should be understood that above general description and following detailed description be only it is exemplary and explanatory, not
The application can be limited.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows the implementation for meeting the application
Example, and together with specification it is used to explain the principle of the application.In the accompanying drawings:
Fig. 1 is a kind of schematic diagram of the method for identification coal seam top shown according to an exemplary embodiment bottom interface；
Fig. 2 is the schematic diagram of relationship between a kind of sound impedance shown according to an exemplary embodiment and natural gamma；
Fig. 3 is a kind of schematic diagram of pseudo-acoustic curve synthesis shown according to an exemplary embodiment；
Fig. 4 is the comparison signal of a kind of original acoustic wave, natural gamma and onomatopoeia wave shown according to an exemplary embodiment
Figure；
Fig. 5 is the schematic diagram of the method for identification coal seam shown according to an exemplary embodiment top bottom interface；
Fig. 6 is a kind of block diagram of the device of identification coal seam top shown according to an exemplary embodiment bottom interface.
Through the above attached drawings, it has been shown that the specific embodiment of the application will be hereinafter described in more detail.These attached drawings
It is not intended to limit the range of the application design in any manner with verbal description, but is by referring to specific embodiments
Those skilled in the art illustrate the concept of the application.
Specific embodiment
Example embodiments are described in detail here, and the example is illustrated in the accompanying drawings.Following description is related to
When attached drawing, unless otherwise indicated, the same numbers in different drawings indicate the same or similar elements.Following exemplary embodiment
Described in embodiment do not represent all embodiments consistent with the application.On the contrary, they be only with it is such as appended
The example of the consistent device and method of some aspects be described in detail in claims, the application.
The embodiment of the present application provides a kind of method of identification coal seam top bottom interface, as shown in Figure 1, which comprises
Step S101 is obtained multiple for reflecting the sensitivity curve of different lithology.
In this step, it is studied by the methods of Seismic reflection character research, Rock physical analysis, a variety of energy can be obtained and know
The sensitivity curve of other different lithology.Since coal seam speed, density are all smaller, coal seam top surface seismic wave is caused to be entered by high speed low
Speed, seismic reflection are high-amplitude wave paddy.On the contrary, coal seam bottom surface is that seismic wave is entered at a high speed by low speed, seismic reflection is strong wave peak.Simultaneously
Coal seam is low speed, low-density, low gamma, high resistance.When curve intersection is analyzed in this way, it is low that coal seam would generally fall in sound wave, density
It is worth end, distinguishes obvious.Very strong sharp pulse is all shown as on sound wave curve, density curve for girdle.Sharp pulse is wide
It spends related with coal seam thickness.But sound wave curve resolution ratio is without density curve height.And density curve does not have coal seam response amplitude
Sound wave curve is big, in the application, is comprehensively considered to the two.Fig. 2 is provided between a kind of sound impedance and natural gamma
The schematic diagram of relationship.
Step S102 generates weighting pseudo-acoustic curve, the weighting onomatopoeia wave according to the sensitivity curve and sound wave curve
Curve is used to reflect the nuance of formation velocity, the variation of wave impedance and lithology.
In this step, by analysis sensitivity curve and coalbed coring relationship, the weight coefficient of calculated curve is melted by information
It closes analysis and obtains the curve in coal seam for identification.Since coal seam shows as this feature of apparent low value on sensitivity curve, together
When consider to divide example, therefore the information of sound wave difference sensitivity curve is fused together by the application proposition.Become using reflection lithology
Change than more sensitive log, construct the new curve with sound wave dimension, in conjunction with the low frequency model of sound wave, synthesizes onomatopoeia curve
Line makes it that can reflect the variation of formation velocity and wave impedance and reflects the nuance of lithology.
Step S103 carries out global optimizing broad-band constrained inversion based on the weighting pseudo-acoustic curve, to identify coal seam
Push up bottom interface.
It is described that weighting onomatopoeia curve is generated according to the sensitivity curve and sound wave curve in a kind of optional embodiment
Line, comprising:
It is weighted based on the sensitivity curve and coalbed coring relationship, the weighting for obtaining coal seam for identification is bent
Line；
The low-frequency information of the high-frequency information of the weighted curve and sound wave curve is merged, the weighting onomatopoeia is obtained
Wave profile.
In this embodiment it is possible to be calculated by the following formula weighting pseudo-acoustic curve:
DTgr (h)=W1*AC (h)+W2*GR (h)
Wherein: DTgr is indicated: pseudo-acoustic curve；
W1, W2 are indicated: weighting coefficient, controlling curve low frequency, high frequency content；
AC is indicated: sound wave curve；
GR is indicated: gamma ray curve.
Fig. 3 is a kind of schematic diagram of pseudo-acoustic curve synthesis shown according to an exemplary embodiment, as shown in figure 3, sound
The high-frequency information fusion of the high-frequency information and reflection lithologic character curve of wave profile, can be obtained pseudo-acoustic curve.The application passes through
The high-frequency information and sound wave (time difference) curve low-frequency information of above-mentioned weighted curve are studied, the onomatopoeia of reflection coalbed coring is established in synthesis
Wave profile.As shown in figure 4, showing the contrast schematic diagram of a kind of original acoustic wave, natural gamma and onomatopoeia wave.Natural gamma in figure
Curve (in) lithology is reflected clearly, and sound wave curve (left side) only has reflection to coal seam, has no to reflect to sandstone, at this moment
We combine the low frequency of sound wave curve using the high frequency (i.e. curve rises and falls, and reflects lithologic character) of gamma ray curve, and information is melted
It is combined, and the dimension of curve new after synthesis is demarcated to sound wave curve dimension range, just obtain pseudo-acoustic curve
(right side).From fig. 4, it can be seen that pseudo-acoustic curve not only reflects sensitivity to lithology, but also there is sound wave dimension, can use acoustic resistance
Antigen reason carries out inverting.
In a kind of optional embodiment, the method also includes:
Using Wavelet Decomposition Technology, the high-frequency information of the weighted curve in coal seam for identification is extracted；And/or
Using Wavelet Decomposition Technology, the low-frequency information of the sound wave curve is extracted.
In this embodiment, the weighted curve of above-mentioned identification coalbed coring is believed its high frequency by Wavelet Decomposition Technology
Breath extracts；By sound wave curve by Wavelet Decomposition Technology, its low-frequency information is extracted；In turn, it is merged using information
Technology gets up weighted curve high-frequency information and the fusion of sound wave curve low-frequency information, synthesis weighting pseudo-acoustic curve.
In a kind of optional embodiment, it is described be weighted by the sensitivity curve and coalbed coring relationship based on
It calculates, obtains the weighted curve in coal seam for identification, comprising:
Nonlinear weight coefficient is calculated according to sensitivity curve and the coal geology feature goodness of fit；
Using the nonlinear weight coefficient, weighted calculation obtains the weighted curve in the coal seam for identification.
In this embodiment, sensitivity curve and coalbed coring relationship are nonlinear function, using sensitivity curve and coal seam
The matter feature goodness of fit determines its contribution.And then nonlinear weight is calculated by sensitivity curve and the coal geology feature goodness of fit
Coefficient, weighted information merge to obtain the curve of coalbed coring for identification.
It is described that the constraint of global optimizing broadband is carried out based on the weighting pseudo-acoustic curve in a kind of optional embodiment
Inverting, to identify the top bottom interface in coal seam, comprising:
Obtain basic data, the basic data includes at least one of: poststack data, modeling data, layer position data,
Wavelet data；
According to the basic data, objective function is established, the objective function is used at least one of: actual seismic note
Record, composite traces, prior-constrained, transverse gradients constraint；
Using the fast inversion algorithm of global optimizing, amendment is iterated to initial geological model, is obtained high-resolution
Sound impedance/quasi-acoustic inversion model；
According to the sound impedance/quasi-acoustic inversion model and the objective function, inversion result is determined；
Based on the inversion result, the top bottom interface in the coal seam is identified.
In this embodiment, the basic data of inverting is inputted, wherein the basic data includes poststack number
According to, modeling data, layer position data, wavelet data etc..According to basic data, objective function is established, including real seismic record,
Composite traces, prior-constrained (it is too far that inversion solution cannot deviate priori value), transverse gradients constraint (have centainly laterally consecutive
Property, keep solution more reasonable).Inverting is carried out using mature general borehole restraint non-linear inversion technology, high-frequency information is made full use of, mentions
The precision and resolution ratio of high inverting carry out iteration repeatedly to initial geological model using the fast inversion algorithm of global optimizing
Amendment, obtains high-resolution sound impedance/quasi-acoustic inversion model.Above-mentioned high-resolution sound impedance/quasi-acoustic inversion
Model corresponds to final inversion result.
Fig. 5 is the schematic diagram of the method for identification coal seam shown according to an exemplary embodiment top bottom interface.Such as Fig. 5 institute
Show, the method for identification coal seam top bottom interface includes:
Step A, the sensitivity curve for distinguishing different lithology is obtained.
It is studied by the methods of Seismic reflection character research, Rock physical analysis, the sensitivity for obtaining distinguishing different lithology is bent
Line.
Step B, it by the above sensitivity curve and coalbed coring relationship analysis, calculates with the weight coefficient of upper curve, information is melted
It closes analysis and obtains the curve in coal seam for identification.
Step C, the high-frequency information and interval transit time curve low-frequency information of above-mentioned weighted curve are studied, synthesis, which is established, to be reflected
The pseudo-acoustic curve of coalbed coring.
Step D, on the basis of above-mentioned weighting pseudo-acoustic curve, global optimizing broad-band constrained inversion, quantitative judge coal are carried out
The top bottom interface of layer.
The above embodiments of the present application can reduce the interference of limitation and multi-solution brought by monotechnics, more smart
Girdle really is portrayed, the top bottom interface of girdle can be identified.The application by coalbed coring sensitivity curve, weight coefficient,
The key links such as information fusion, onomatopoeia wave, inverting are together in series, and a variety of valuable informations are fused together, single skill is reduced
The interference of limitation and multi-solution brought by art.Girdle top bottom interface is predicted by this method, it is objective clear, with coal seam
Geology cognition is consistent.The advantage that the inversion method is high, with a high credibility with precision, has a wide range of application.
The application another exemplary embodiment provides a kind of device of identification coal seam top bottom interface, as shown in fig. 6, the dress
It sets and includes:
Acquiring unit 60, it is multiple for reflecting the sensitivity curve of different lithology for obtaining；
Generation unit 61, for generating weighting pseudo-acoustic curve, the weighting according to the sensitivity curve and sound wave curve
Pseudo-acoustic curve is used to reflect the nuance of formation velocity, the variation of wave impedance and lithology；
Inverting unit 62, for carrying out global optimizing broad-band constrained inversion based on the weighting pseudo-acoustic curve, with identification
The top bottom interface in coal seam.
In a kind of optional embodiment, the generation unit includes:
Computing module is obtained for being weighted based on the sensitivity curve and coalbed coring relationship for knowing
The weighted curve in other coal seam；
Fusion Module is obtained for merging the low-frequency information of the high-frequency information of the weighted curve and sound wave curve
To the weighting pseudo-acoustic curve.
In a kind of optional embodiment, described device further include:
First extraction module extracts the height of the weighted curve in coal seam for identification for utilizing Wavelet Decomposition Technology
Frequency information；And/or
Second extraction module extracts the low-frequency information of the sound wave curve for utilizing Wavelet Decomposition Technology.
In a kind of optional embodiment, the computing module includes:
First computational submodule, for calculating nonlinear weight system according to sensitivity curve and the coal geology feature goodness of fit
Number；
Second computational submodule, for utilizing the nonlinear weight coefficient, weighted calculation obtains the coal for identification
The weighted curve of layer.
In a kind of optional embodiment, the inverting unit includes:
Basic data obtains module, and for obtaining basic data, the basic data includes at least one of: poststack number
According to, modeling data, layer position data, wavelet data；
Objective function establishes module, for according to the basic data, establishing objective function, the objective function be used for
It is at least one lower: real seismic record, composite traces, prior-constrained, transverse gradients constraint；
Iterated revision module is iterated initial geological model and repairs for the fast inversion algorithm using global optimizing
Just, high-resolution sound impedance/quasi-acoustic inversion model is obtained；
Inversion result determining module is used for according to the sound impedance/quasi-acoustic inversion model and the objective function, really
Determine inversion result；
Bottom interface identification module is pushed up, for being based on the inversion result, identifies the top bottom interface in the coal seam.
It is multiple for reflecting the sensitivity curve of different lithology by obtaining in the embodiment of the present application, according to described sensitive bent
Line and sound wave curve generate weighting pseudo-acoustic curve, and the weighting pseudo-acoustic curve is used to reflect the change of formation velocity, wave impedance
Change the nuance with lithology, and global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify coal
The top bottom interface of layer improves the identification accuracy of the top bottom interface in coal seam, reduces limitation and more brought by monotechnics
The interference of solution property.
Those skilled in the art will readily occur to its of the application after considering specification and practicing application disclosed herein
Its embodiment.This application is intended to cover any variations, uses, or adaptations of the application, these modifications, purposes or
The common knowledge in the art that person's adaptive change follows the general principle of the application and do not apply including the application
Or conventional techniques.The description and examples are only to be considered as illustrative, and the true scope and spirit of the application are by following
Claim is pointed out.
It should be understood that the application is not limited to the precise structure that has been described above and shown in the drawings, and
And various modifications and changes may be made without departing from the scope thereof.Scope of the present application is only limited by the accompanying claims.
Claims (10)
1. a kind of method of identification coal seam top bottom interface, which is characterized in that the described method includes:
It obtains multiple for reflecting the sensitivity curve of different lithology；
Weighting pseudo-acoustic curve is generated according to the sensitivity curve and sound wave curve, the weighting pseudo-acoustic curve is for reflecting ground
Interval velocity, the variation of wave impedance and lithology nuance；
Global optimizing broad-band constrained inversion is carried out based on the weighting pseudo-acoustic curve, to identify the top bottom interface in coal seam.
2. the method according to claim 1, wherein described generated according to the sensitivity curve and sound wave curve adds
Weigh pseudo-acoustic curve, comprising:
It is weighted based on the sensitivity curve and coalbed coring relationship, obtains the weighted curve in coal seam for identification；
The low-frequency information of the high-frequency information of the weighted curve and sound wave curve is merged, the weighting onomatopoeia curve is obtained
Line.
3. according to the method described in claim 2, it is characterized in that, the method also includes:
Using Wavelet Decomposition Technology, the high-frequency information of the weighted curve in coal seam for identification is extracted；And/or
Using Wavelet Decomposition Technology, the low-frequency information of the sound wave curve is extracted.
4. according to the method described in claim 2, it is characterized in that, described be based on the sensitivity curve and coalbed coring relationship
It is weighted, obtains the weighted curve in coal seam for identification, comprising:
Nonlinear weight coefficient is calculated according to sensitivity curve and the coal geology feature goodness of fit；
Using the nonlinear weight coefficient, weighted calculation obtains the weighted curve in the coal seam for identification.
5. the method according to claim 1, wherein described sought based on the weighting pseudo-acoustic curve progress overall situation
Excellent broad-band constrained inversion, to identify the top bottom interface in coal seam, comprising:
Basic data is obtained, the basic data includes at least one of: poststack data, modeling data, layer position data, wavelet
Data；
According to the basic data, objective function is established, the objective function is at least one of: real seismic record,
Composite traces, prior-constrained, transverse gradients constraint；
Using the fast inversion algorithm of global optimizing, amendment is iterated to initial geological model, obtains high-resolution sound wave
Impedance/quasi-acoustic inversion model；
According to the sound impedance/quasi-acoustic inversion model and the objective function, inversion result is determined；
Based on the inversion result, the top bottom interface in the coal seam is identified.
6. a kind of device of identification coal seam top bottom interface, which is characterized in that described device includes:
Acquiring unit, it is multiple for reflecting the sensitivity curve of different lithology for obtaining；
Generation unit, for generating weighting pseudo-acoustic curve, the weighting onomatopoeia wave according to the sensitivity curve and sound wave curve
Curve is used to reflect the nuance of formation velocity, the variation of wave impedance and lithology；
Inverting unit, for carrying out global optimizing broad-band constrained inversion based on the weighting pseudo-acoustic curve, to identify coal seam
Push up bottom interface.
7. device according to claim 6, which is characterized in that the generation unit includes:
Computing module obtains coal for identification for being weighted based on the sensitivity curve and coalbed coring relationship
The weighted curve of layer；
Fusion Module obtains institute for merging the low-frequency information of the high-frequency information of the weighted curve and sound wave curve
State weighting pseudo-acoustic curve.
8. device according to claim 7, which is characterized in that described device further include:
First extraction module extracts the high frequency letter of the weighted curve in coal seam for identification for utilizing Wavelet Decomposition Technology
Breath；And/or
Second extraction module extracts the low-frequency information of the sound wave curve for utilizing Wavelet Decomposition Technology.
9. device according to claim 7, which is characterized in that the computing module includes:
First computational submodule, for calculating nonlinear weight coefficient according to sensitivity curve and the coal geology feature goodness of fit；
Second computational submodule, for utilizing the nonlinear weight coefficient, weighted calculation obtains the coal seam for identification
Weighted curve.
10. device according to claim 6, which is characterized in that the inverting unit includes:
Basic data obtains module, and for obtaining basic data, the basic data includes at least one of: poststack data,
Modeling data, layer position data, wavelet data；
Objective function establishes module, for according to the basic data, establishing objective function, the objective function be used for down toward
It is one of few: real seismic record, composite traces, prior-constrained, transverse gradients constraint；
Iterated revision module is iterated amendment to initial geological model, obtains for the fast inversion algorithm using global optimizing
To high-resolution sound impedance/quasi-acoustic inversion model；
Inversion result determining module, for determining anti-according to the sound impedance/quasi-acoustic inversion model and the objective function
Drill result；
Bottom interface identification module is pushed up, for being based on the inversion result, identifies the top bottom interface in the coal seam.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810792633.9A CN109143326A (en) | 2018-07-18 | 2018-07-18 | The method and apparatus for identifying coal seam top bottom interface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810792633.9A CN109143326A (en) | 2018-07-18 | 2018-07-18 | The method and apparatus for identifying coal seam top bottom interface |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109143326A true CN109143326A (en) | 2019-01-04 |
Family
ID=64800840
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810792633.9A Pending CN109143326A (en) | 2018-07-18 | 2018-07-18 | The method and apparatus for identifying coal seam top bottom interface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143326A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112394392A (en) * | 2019-08-13 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and device for evaluating distribution condition of hydrocarbon source rock |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150260034A1 (en) * | 2014-03-14 | 2015-09-17 | Schlumberger Technology Corporation | Determining mineralogy of an earth formation using linear regressions |
| CN105089652A (en) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | Pseudo-acoustic curve rebuilding and sparse pulse joint inversion method |
| CN105223608A (en) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | A kind of earthquake prediction and describing method shielding by force fractured-vuggy reservoir containing coal |
| CN106526696A (en) * | 2016-11-03 | 2017-03-22 | 中国矿业大学(北京) | Coal seam roof and floor lithology prediction method |
| CN107817522A (en) * | 2017-10-31 | 2018-03-20 | 西安恒泰艾普能源发展有限公司 | A kind of Deterministic Methods of marine qualitative recognition hydrate top bottom interface |
-
2018
- 2018-07-18 CN CN201810792633.9A patent/CN109143326A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150260034A1 (en) * | 2014-03-14 | 2015-09-17 | Schlumberger Technology Corporation | Determining mineralogy of an earth formation using linear regressions |
| CN105089652A (en) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | Pseudo-acoustic curve rebuilding and sparse pulse joint inversion method |
| CN105223608A (en) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | A kind of earthquake prediction and describing method shielding by force fractured-vuggy reservoir containing coal |
| CN106526696A (en) * | 2016-11-03 | 2017-03-22 | 中国矿业大学(北京) | Coal seam roof and floor lithology prediction method |
| CN107817522A (en) * | 2017-10-31 | 2018-03-20 | 西安恒泰艾普能源发展有限公司 | A kind of Deterministic Methods of marine qualitative recognition hydrate top bottom interface |
Non-Patent Citations (12)
| Title |
|---|
| 付志方等: "拟声波构建技术在砂泥岩薄互层储层预测中的应用", 《石油物探》 * |
| 刘代志: "《地球物理环境与国家安全》", 31 December 2014, 西安：西安地图出版社 * |
| 周连敏等: "拟声波地震反演在三角洲前缘河道砂体预测中的应用", 《勘探开发》 * |
| 安晶等: "多曲线融合技术在储层随机反演中的应用", 《物探与化探》 * |
| 尹继尧等: "多测井曲线拟声波重构技术的应用", 《新疆石油地质》 * |
| 张静等: "储层特征曲线重构技术在储层预测中的应用研究", 《天然气地球科学》 * |
| 彭小东等: "拟声波地震反演技术在合川气田须二段储层预测中的应用", 《内江科技》 * |
| 王允诚: "《油气储层地质学》", 31 January 2008, 北京：地质出版社 * |
| 王香增: "《鄂尔多斯盆地南部天然气勘探》", 31 March 2016, 北京：地质出版社 * |
| 赵继龙等: "伽马拟声波分频重构反演在储层预测中的应用", 《地球物理学进展》 * |
| 陈强等: "《沁水盆地高瓦斯矿井三维地震技术应用实践》", 30 September 2012, 北京：煤炭工业出版社 * |
| 魏文希等: "拟声波反演技术在识别煤层顶底板砂泥岩中的应用", 《物探与化探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112394392A (en) * | 2019-08-13 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and device for evaluating distribution condition of hydrocarbon source rock |
| CN112394392B (en) * | 2019-08-13 | 2023-09-26 | 中国石油天然气股份有限公司 | Method and device for evaluating distribution condition of hydrocarbon source rock |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8930170B2 (en) | Attribute importance measure for parametric multivariate modeling | |
| US20200132869A1 (en) | Seismic rock physics inversion method based on large area tight reservoir | |
| CN101680960B (en) | Optimizing amplitude inversion utilizing statistical comparisons of seismic to well control data | |
| CN101999086A (en) | Method for determining seismic data quality | |
| Verma et al. | The early Paleozoic structures and its influence on the Permian strata, Midland Basin: Insights from multi-attribute seismic analysis | |
| Lehocki et al. | Seismic methods for fluid discrimination in areas with complex geologic history—A case example from the Barents Sea | |
| Pancha et al. | Determination of 3D basin shear‐wave velocity structure using ambient noise in an urban environment: A case study from Reno, Nevada | |
| Kulikowski et al. | 3D seismic analysis investigating the relationship between stratigraphic architecture and structural activity in the intra-cratonic Cooper and Eromanga basins, Australia | |
| Burgin et al. | Determining paleo-structural environments through natural fracture and calcite twin analyses: a case study in the Otway Basin, Australia | |
| CN109143326A (en) | The method and apparatus for identifying coal seam top bottom interface | |
| Grit et al. | Integrated seismic survey for detecting landslide effects on high speed rail line at Istanbul–Turkey | |
| AlKawai et al. | Integrating statistical rock physics and pressure and thermal history modeling to map reservoir lithofacies in the deepwater Gulf of Mexico | |
| Zellou et al. | Integrated fractured reservoir characterization using neural networks and fuzzy logic: three case studies | |
| Fu et al. | Rock property-and seismic-attribute analysis of a chert reservoir in the Devonian Thirty-one Formation, west Texas, USA | |
| Chellingsworth et al. | Human factors in seismic uncertainty—Restoring a realistic uncertainty range | |
| Wotherspoon et al. | Site period characteristics across the Canterbury region of New Zealand | |
| RU22830U1 (en) | AUTOMATED SYSTEM FOR SEARCHING, EXPLORING AND EVALUATING THE OPERATIONAL PROPERTIES OF DEPOSITS AND DEPOSITS OF MINERAL RESOURCES AND FORECAST OF TECTONIC AND LITHOPHYSICAL PROPERTIES OF GEOLOGICAL MEDIA | |
| Sarzalejo et al. | Stratigraphy and lithologic heterogeneity in the Mannville Group (southeast Saskatchewan) defined by integrating 3-D seismic and log data | |
| Clemons et al. | Seismic attributes: Exploiting seismic data to understand heterogeneous reservoir performance in the Eagle Ford Shale, south Texas, USA | |
| Foley et al. | Fast-Tracking Gold Exploration Below 300m around a mature mine complex–3D Seismic Case History of the Darlot–Centenary Gold Mine | |
| Rathore et al. | Prediction of petrophysical properties through comparative post-stack inversion techniques using advance Neural Networking | |
| Langlais et al. | Comparison of several methods of lithofacies simulation on the fluvial gypsy sandstone of Oklahoma | |
| Gutierrez | Rock physics-based interpretation framework: Key to de-risking amplitude-supported exploration plays & prospects, offshore Mexico | |
| Lupo et al. | A case study in the Pennsylvanian Cleveland sandstone on the Nemaha Ridge: Leveraging high resolution 3D seismic & stratigraphic analysis to create the conditions for repeatable drilling opportunities in complexly distributed reservoirs | |
| Stephens et al. | 3D Reservoir Characterization and Integrated Completion Optimization for Understanding Horizontal Well Spacing, Frac Staging, and Well Performance in the Niobrara Formation, DJ Basin |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhang Yan Inventor after: Zhang Zhirang Inventor after: Wang Zhaofeng Inventor after: Tai Lixun Inventor after: Zhao Shifeng Inventor after: Cui Gang Inventor before: Zhang Yan Inventor before: Zhang Zhirang Inventor before: Sun Gengwen Inventor before: Wang Zhaofeng Inventor before: Cui Gang Inventor before: Tang Chengfeng |
|
| CB03 | Change of inventor or designer information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190104 |
|
| RJ01 | Rejection of invention patent application after publication |