CN110297270A - High-Resolution Seismic Data method based on structure constraint - Google Patents

High-Resolution Seismic Data method based on structure constraint Download PDF

Info

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
Application number
CN201910496190.3A
Other languages
Chinese (zh)
Inventor
孙振刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Youlong Technology Service Co Ltd
Original Assignee
Beijing Youlong Technology Service Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Youlong Technology Service Co Ltd filed Critical Beijing Youlong Technology Service Co Ltd
Priority to CN201910496190.3A priority Critical patent/CN110297270A/en
Publication of CN110297270A publication Critical patent/CN110297270A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing 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

High-Resolution Seismic Data method based on structure constraint
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.
CN201910496190.3A 2019-06-10 2019-06-10 High-Resolution Seismic Data method based on structure constraint Pending CN110297270A (en)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111060969A (en) * 2019-12-25 2020-04-24 恒泰艾普(北京)能源科技研究院有限公司 Well control Q compensation method

Citations (7)

* Cited by examiner, † Cited by third party
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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