CN110297270A - High-Resolution Seismic Data method based on structure constraint - Google Patents
High-Resolution Seismic Data method based on structure constraint Download PDFInfo
- Publication number
- CN110297270A CN110297270A CN201910496190.3A CN201910496190A CN110297270A CN 110297270 A CN110297270 A CN 110297270A CN 201910496190 A CN201910496190 A CN 201910496190A CN 110297270 A CN110297270 A CN 110297270A
- Authority
- CN
- China
- Prior art keywords
- resolution
- data
- signal
- seismic data
- structure constraint
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000010010 raising Methods 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 6
- 238000003325 tomography Methods 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000003908 quality control method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012790 confirmation Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000003672 processing method Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 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
- 238000002360 preparation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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
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 High-Resolution Seismic Data method based on structure constraint that the invention discloses a kind of, it is related to Seismic Data Processing Technique field, itself the following steps are included: S1, from VSP well logging in obtain Q value, and then the Q value of individual well and multiple wells is shifted by layer position structure constraint by full work area, space-variant Q body when formation, then longitudinal resolution is improved by the method for Q compensation, S2, carries out random Beta deconvolution on the basis of previous step, can further improves longitudinal resolution.The High-Resolution Seismic Data method based on structure constraint, it makes full use of the wells data such as existing acoustic logging, VSP to carry out well control and improves resolution processes, accomplish the raising resolution ratio for having foundation, the variation of signal-to-noise ratio is considered while improving resolution ratio, the potentiality of data are sufficiently excavated using Tau-P domain construction constraint filtering technique, accomplish resolution ratio and the double experiences for improving, recognizing using the personnel of explanation data of signal-to-noise ratio, accomplishes sufficient integrated processing and interpretation during processing.
Description
Technical field
The present invention relates to Seismic Data Processing Technique field, specially a kind of high-resolution seismic exploration money based on structure constraint
Expect processing method.
Background technique
The low problem of conventional seism processing longitudinal resolution, and unlimited raising longitudinal resolution will lead to following three
Kind result: (1) data signal-to-noise ratio decreases after longitudinal resolution improves, quality of achievement decline;(2) longitudinal resolution improves
Too much, false construction is generated;(3) cause tomographic imaging to be deteriorated after improving resolution ratio, and then the lateral resolution of data is caused to reduce.
Achievement data is influenced by firsthand information quality and pre-stack processing method, some are asked there are resolution ratio is relatively low
Topic, is mainly reflected in that longitudinal resolution is lower, for example can not identify the sand body of thinner thickness, to be difficult to portray the spread of sand body
Rule influences the difficulty of direction in space identification.But the unconfined point longitudinal resolution that improves can generate false leveling signal, simultaneously
Since the quality of data high frequency signal is limited, the reduction of whole signal-to-noise ratio is easy to cause when improving the specific gravity of high-frequency signal.
It is badly in need of the new high-resolution technology that proposes of one kind now and solves problem above.
Summary of the invention
(1) the technical issues of solving
In view of the deficiencies of the prior art, the present invention provides a kind of High-Resolution Seismic Datas based on structure constraint
It is low to solve conventional seism processing longitudinal resolution for method, but the unconfined point longitudinal resolution that improves can generate falseness
Leveling signal be easy to cause when improving the specific gravity of high-frequency signal whole simultaneously because the quality of data high frequency signal is limited
The problem of reduction of body signal-to-noise ratio.
(2) technical solution
To achieve the above objectives, the technical solution adopted by the present invention is that: a kind of high-resolution seismic exploration based on structure constraint
Data processing method, comprising the following steps:
S1, Q value is obtained from VSP well logging, and then the Q value of individual well and multiple wells is shifted onto entirely by layer position structure constraint
Work area, space-variant Q body when formation, the method for then passing through Q compensation improve longitudinal resolution.
S2, random Beta deconvolution is carried out on the basis of previous step, can further improve longitudinal resolution, still
It can not confirm the reasonability of processing result, in order to there is the raising resolution ratio of foundation, need to obtain mesh from the information of acoustic logging
Layer reflection coefficient, to carry out Quality Control to the result of deconvolution, it is ensured that the high-resolution data acquired are genuine and believable.
S3, signal-to-noise ratio is improved after upper two step improves resolution ratio using newest homologous Fourier transformation antialiasing technology,
Obtain the higher seismic data of quality.
S4, after completing above-mentioned 3 step, after improving signal-to-noise ratio and having the raising longitudinal resolution of foundation, high-frequency signal is more preferably
It is abundant, and low frequency signal is no longer dominant relatively, therefore not can guarantee tomography then and be still apparent, so needing using new Tau-
The domain P filtering technique can guarantee that tomography is clear, and lateral resolution improves simultaneously.
In addition the processing step of data and layer data is explained performance data, layer position are as follows:
T1, the performance data of input, layer position are explained that data and layer data carry out encryption interval sampling, then is sought
Gain factor is simultaneously applied.
T2, by data obtained in step T1 carry out the domain Tau-P filtering remove dryness with the domain Tau-P filtering signal enhance, thus
So that amplitude enhancing and frequency equilibrium.
T3, the data for obtaining step T2 carry out it is counter fall gain factor, and according to crude sampling rate progress resampling.
Preferably, the step S1 to S4 can make data explain that personnel obtain signal-to-noise ratio, longitudinal resolution and transverse direction
The data volume for the Trinity high quality that resolution ratio all improves.
Preferably, rapid S1 Q into S4 is formation absorption coefficient, and Tau-P is a kind of point domain mode of seismic data.
(3) beneficial effect
The beneficial effects of the present invention are:
1, it is somebody's turn to do the High-Resolution Seismic Data method based on structure constraint, makes full use of existing acoustic logging, VSP
Equal wells data carries out well control and improves resolution processes, accomplishes the raising resolution ratio for having foundation, considers while improving resolution ratio
The variation of signal-to-noise ratio is sufficiently excavated the potentiality of data using Tau-P domain construction constraint filtering technique, accomplishes resolution ratio and signal-to-noise ratio
Double raisings, the experience recognized using the personnel of explanation data, accomplish sufficient integrated processing and interpretation during processing.
2, it is somebody's turn to do the High-Resolution Seismic Data method based on structure constraint, obtains true subsurface structure, avoids beating
Well positional fault, so that the economic loss for the well that avoids fanning the air, beats the cost of a well at least 2000 in central/western China
It is more than ten thousand RMB, if with new technology can to avoid the above economic loss, new achievement signal-to-noise ratio can be better than currently at
Fruit can find new trap and oil reservoir on the basis of new results, make more contributions to the energy security of country.
Detailed description of the invention
Fig. 1 is data amplitudes frequency processing method flow schematic diagram of the present invention;
Fig. 2 is tomography site location cross-reference schematic diagram of the present invention;
Fig. 3 is data amplitudes frequency variation curve schematic diagram of the present invention;
Fig. 4 is the first Tau-P domain construction constraint filtering front and back of the invention contrast schematic diagram;
Fig. 5 is the 2nd Tau-P domain construction constraint filtering front and back of the invention contrast schematic diagram;
Fig. 6 is that structure constraint of the present invention filters front and back correlation curve schematic diagram;
Fig. 7 is that log data of the present invention generates composite traces and seismic trace near well and matches schematic diagram.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
As shown in figs. 1-7, the present invention provides a kind of technical solution: a kind of high-resolution seismic data based on structure constraint
Processing method, comprising the following steps:
S1, Q value is obtained from VSP well logging, and then the Q value of individual well and multiple wells is shifted onto entirely by layer position structure constraint
Work area, space-variant Q body when formation, the method for then passing through Q compensation improve longitudinal resolution.
S2, random Beta deconvolution is carried out on the basis of previous step, can further improve longitudinal resolution, still
It can not confirm the reasonability of processing result, in order to there is the raising resolution ratio of foundation, need to obtain mesh from the information of acoustic logging
Layer reflection coefficient, to carry out Quality Control to the result of deconvolution, it is ensured that the high-resolution data acquired are genuine and believable.
S3, signal-to-noise ratio is improved after upper two step improves resolution ratio using newest homologous Fourier transformation antialiasing technology,
Obtain the higher seismic data of quality.
S4, after completing above-mentioned 3 step, after improving signal-to-noise ratio and having the raising longitudinal resolution of foundation, high-frequency signal is more preferably
It is abundant, and low frequency signal is no longer dominant relatively, therefore not can guarantee tomography then and be still apparent, so needing using new Tau-
The domain P filtering technique can guarantee that tomography is clear, and lateral resolution improves simultaneously, and step S1 to S4 can make data explain personnel
The data volume of signal-to-noise ratio, the Trinity high quality that longitudinal resolution and lateral resolution all improve is obtained, wherein Q is stratum
Absorption coefficient, Tau-P are a kind of point domain modes of seismic data.
In addition the processing step of data and layer data is explained performance data, layer position are as follows:
T1, the performance data of input, layer position are explained that data and layer data carry out encryption interval sampling, then is sought
Gain factor is simultaneously applied.
T2, by data obtained in step T1 carry out the domain Tau-P filtering remove dryness with the domain Tau-P filtering signal enhance, thus
So that amplitude enhancing and frequency equilibrium.
T3, the data for obtaining step T2 carry out it is counter fall gain factor, and according to crude sampling rate progress resampling.
Wherein the first lines are that Q is equal to 50 from the bottom up when TWT, that is, 250ms in Fig. 3, and the second lines Q is equal to
100, third lines Q are equal to 150, and the 4th lines Q is equal to 200, and front end attenuation amplitude is greater than low frequency end, and target zone is deeper, more
It is the longitudinal resolution for influencing data:
A (ω, τ)=A (ω, 0) exp (- ω τ/2Q)
Scheme left to be the application of the prior art as a result, can not confirm the accurate location at tomography position in Fig. 2, at the same lineups it
Between can not distinguish, longitudinal resolution is lower, right figure be invention new technology application effect, interrupting layer position is high-visible,
What layer interdigit can also be apparent is distinguished.
It is carried out with performance data using VSP corridor stack corresponding in Fig. 7, is generated by composite traces and well using log data
Seismic channel matches.
Theoretically longitudinal resolution can be improved infinitely, but how accomplish the raising resolution ratio for having foundation, mainly be used
The mode of well control:
Q value is extracted from VSP, it is anti-to generate the random Beta of Beta factor progress using sound wave measuring well curve for Lai Jinhang Q compensation
Convolution, is carried out using VSP corridor stack with performance data corresponding, and log data is utilized to generate composite traces and seismic trace near well phase
Matching.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects
It is described in detail, it should be understood that being not limited to this hair the foregoing is merely a specific embodiment of the invention
Bright, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in the present invention
Protection scope within.
Claims (3)
1. a kind of High-Resolution Seismic Data method based on structure constraint, it is characterised in that: the following steps are included:
S1, Q value is obtained from VSP well logging, and then the Q value of individual well and multiple wells is shifted by layer position structure constraint by full work
Area, space-variant Q body when formation, the method for then passing through Q compensation improve longitudinal resolution;
S2, random Beta deconvolution is carried out on the basis of previous step, can further improve longitudinal resolution, but can not
The reasonability of confirmation processing result needs to obtain target zone from the information of acoustic logging to there is the raising resolution ratio of foundation
Reflection coefficient, to carry out Quality Control to the result of deconvolution, it is ensured that the high-resolution data acquired are genuine and believable;
S3, signal-to-noise ratio is improved after upper two step improves resolution ratio using newest homologous Fourier transformation antialiasing technology, obtained
The higher seismic data of quality;
S4, after completing above-mentioned 3 step, after improving signal-to-noise ratio and having the raising longitudinal resolution of foundation, high-frequency signal is more preferably abundant,
And low frequency signal is no longer dominant relatively, therefore not can guarantee tomography then and be still apparent, so needing to filter using the new domain Tau-P
Wave technology can guarantee that tomography is clear, and lateral resolution improves simultaneously;
In addition the processing step of data and layer data is explained performance data, layer position are as follows:
T1, the performance data of input, layer position are explained that data and layer data carry out encryption interval sampling, then carries out seeking gain
The factor is simultaneously applied;
T2, by data obtained in step T1 carry out the domain Tau-P filtering remove dryness with the domain Tau-P filtering signal enhance so that
Amplitude enhancing and frequency equilibrium;
T3, the data for obtaining step T2 carry out it is counter fall gain factor, and according to crude sampling rate progress resampling.
2. the High-Resolution Seismic Data method according to claim 1 based on structure constraint, it is characterised in that: institute
Stating step S1 to S4 can make data explain that personnel obtain signal-to-noise ratio, three that longitudinal resolution and lateral resolution all improve
The data volume of the integrated high quality in position.
3. the High-Resolution Seismic Data method according to claim 1 based on structure constraint, it is characterised in that: institute
Stating rapid S1 Q into S4 is formation absorption coefficient, and Tau-P is a kind of point domain mode of seismic data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910496190.3A CN110297270A (en) | 2019-06-10 | 2019-06-10 | High-Resolution Seismic Data method based on structure constraint |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910496190.3A CN110297270A (en) | 2019-06-10 | 2019-06-10 | High-Resolution Seismic Data method based on structure constraint |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110297270A true CN110297270A (en) | 2019-10-01 |
Family
ID=68027713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910496190.3A Pending CN110297270A (en) | 2019-06-10 | 2019-06-10 | High-Resolution Seismic Data method based on structure constraint |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110297270A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普(北京)能源科技研究院有限公司 | Well control Q compensation method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5426618A (en) * | 1993-05-03 | 1995-06-20 | Chen; Hong-Bin | Method of high resolution and high SNR data acquisition for probing using pulse-compression |
US5970023A (en) * | 1999-02-19 | 1999-10-19 | Baker Hughes Incorporated | Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms |
WO2007075295A2 (en) * | 2005-12-10 | 2007-07-05 | Landray Technology, Inc. | Linear fm radar |
US20120215453A1 (en) * | 2011-02-22 | 2012-08-23 | Cggveritas Services Sa | Device and method for multi-dimensional coherency driven denoising data |
CN105652319A (en) * | 2016-01-19 | 2016-06-08 | 中国海洋石油总公司 | Estimation method of near-surface stratum Q value of complex mediums |
CN105954801A (en) * | 2016-06-12 | 2016-09-21 | 成都理工大学 | Frequency-conversion-attribute-based reservoir permeability evaluation method |
CN106154330A (en) * | 2016-06-20 | 2016-11-23 | 吉林大学 | Electromagnetic type controlled source parabolic Radon transform suppresses parallel focus crosstalk noise method |
-
2019
- 2019-06-10 CN CN201910496190.3A patent/CN110297270A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5426618A (en) * | 1993-05-03 | 1995-06-20 | Chen; Hong-Bin | Method of high resolution and high SNR data acquisition for probing using pulse-compression |
US5970023A (en) * | 1999-02-19 | 1999-10-19 | Baker Hughes Incorporated | Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms |
WO2007075295A2 (en) * | 2005-12-10 | 2007-07-05 | Landray Technology, Inc. | Linear fm radar |
US20120215453A1 (en) * | 2011-02-22 | 2012-08-23 | Cggveritas Services Sa | Device and method for multi-dimensional coherency driven denoising data |
CN105652319A (en) * | 2016-01-19 | 2016-06-08 | 中国海洋石油总公司 | Estimation method of near-surface stratum Q value of complex mediums |
CN105954801A (en) * | 2016-06-12 | 2016-09-21 | 成都理工大学 | Frequency-conversion-attribute-based reservoir permeability evaluation method |
CN106154330A (en) * | 2016-06-20 | 2016-11-23 | 吉林大学 | Electromagnetic type controlled source parabolic Radon transform suppresses parallel focus crosstalk noise method |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普(北京)能源科技研究院有限公司 | Well control Q compensation method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE60212575T2 (en) | METHOD FOR SPECTRAL EQUALIZATION OF SEISMIC DATA WITH OFFSET | |
CN109738945B (en) | Method for directly generating construction diagram by using prestack depth migration result | |
CN106526678B (en) | A kind of wave field separation method and device of reflected acoustic wave well logging | |
CN103336303A (en) | Method of conducting seismic frequency extension by utilizing acoustic logging | |
CN107356964A (en) | Q value estimation and compensation method of the S-transformation domain based on variation principle | |
CN106033125B (en) | The trace gather interference of compacting prestack wide-angle carries frequency method | |
CN103675904B (en) | Explanatory target processing method and processing device is mated in a kind of well shake | |
EP1068545B1 (en) | Ground roll attenuation method | |
CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
CN107703546A (en) | A kind of new threshold function table seismic data denoising method based on wavelet transformation | |
CN104635264B (en) | The processing method of earthquake data before superposition and equipment | |
CN110297270A (en) | High-Resolution Seismic Data method based on structure constraint | |
CN109164492B (en) | Method for extracting sound wave velocity of cased well stratum | |
CN100412569C (en) | Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering | |
CN105093291B (en) | A method of restoring oil and gas reservoir Seismic reflection character | |
CN109471162A (en) | Interbed multiple processing method, system, electronic equipment and readable medium | |
CN106371140B (en) | A kind of method for improving mid-deep strata seismic data resolution | |
CN101285892B (en) | Wide band seismic exploration method | |
CN107728213A (en) | A kind of new threshold function table seismic data denoising method of small echo | |
CN106597549A (en) | Reflection wave extracting method and device, and acoustic reflection image logging instrument | |
Chen et al. | Improving the Precision of Surface Seismic Data Processing by Walkaway VSP | |
CN109975873A (en) | A kind of method and system of reverse-time migration imaging removal low frequency noise | |
CN112213775B (en) | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data | |
Yaskevich et al. | A case study on receiver-clamping quality assessment from the seismic-interferometry processing of downhole seismic noise recordings | |
Zhang et al. | Study and practice of wide-angle seismic data processing |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191001 |