CN105093312B - Earthquake Relative Wave Impedance Forecasting Methodology and device based on frequency domain pluridifferentiation - Google Patents

Abstract

The embodiment of the present application provides a kind of earthquake Relative Wave Impedance Forecasting Methodology and device based on frequency domain pluridifferentiation.The method includes：Thin interbed model is set up, and synthetic seismogram is obtained according to the thin interbed model.Then, the high frequency lifting of not same order differential is carried out to the synthetic seismogram in frequency domain, and the relation set up between amplitude spectrum differential order and dominant frequency, by earthquake record not same order differential amplitude spectrum addition obtain overlapping amplitude spectrum, and the well logging carried out to earthquake record pluridifferentiation amplitude spectrum by well logging sound wave impedance spectrum is demarcated, overlapping amplitude spectrum after being corrected, three rank differential phases spectrum in conjunction with synthetic seismogram obtains wide frequency-band spectrum, and the high-resolution earthquake relative impedances obtained according to wide frequency-band spectrum, such that it is able to effectively identify thin interbed.The resolution ratio of earthquake Relative Wave Impedance can be effectively improved using technical scheme.

Description

Earthquake Relative Wave Impedance Forecasting Methodology and device based on frequency domain pluridifferentiation

Technical field

The present invention relates to technical field of geophysical exploration, more particularly to a kind of seismic facies based on frequency domain pluridifferentiation To wave impedance Forecasting Methodology and device.

Background technology

Quarter-wave of the thickness much smaller than seismic wavelet in thin interbed, causes the reflective information of different layers in earthquake Interfered on record, severe jamming for thin interbed imaging and identification, therefore, the identification of thin interbed is always that earthquake is surveyed The problem visited.

The technologies such as the existing deconvolution to earthquake record are required to be acted on using the operator with certain time window length Earthquake record, to improve the resolution ratio of earthquake record, and (tens to several by operator length to have the operator of certain time window length Hundred milliseconds, equivalent to tens to hundreds of meters of thickness) in formation information mixed processing, seriously limit deconvolution technology and improve ground The ability of log resolution is shaken, cannot also be finally inversed by can effectively recognize the high-resolution seismic exploration Relative Wave Impedance of thin interbed.

However, differential or derivation can but ask for the rate of change between adjacent 2 points, the limitation of operator length is breached.It is micro- Since sub-argument opinion is founded from newton and Leibnitz, developed rapidly, and obtained in fields such as astronomy, biology, military affairs It is widely applied, but rarely has people to use in the data treatment of field of seismic exploration.Differential can by disclose rate of change come The detailed information of prominent thin layer.Zhao in 1984 it is holy it is bright etc. propose a kind of method that seismic frequency is compensated using differential method, Energy to each layer after image processing of seismic section is kept substantially resolution ratio that is consistent and improve data.The U.S. thickness of cloud is to micro- within 1997 Point can improve seismic resolution comments.The spatial domain differential method is applied to three dimensional seismic data by She De equalitys within 1998 In, achieve preferable effect.Muhammad S in 2014 etc. propose one kind and fast and effectively calculate based on cascade Dipolar Filter Method, several rank differential transforms have been made by earthquake record, are added in the homophasic condition, improve the resolution ratio of seismic data.But It is that above method is essentially confined to differential in single-order or several rank time-domains, high boost effect is limited, at its later stage inverting Manage and can not effectively recognize thin interbed in the wave impedance section for obtaining.

The content of the invention

The purpose of the application is to provide a kind of earthquake Relative Wave Impedance Forecasting Methodology and dress based on frequency domain pluridifferentiation Put, using the larger distinction of signal and noise degree of frequency domain pluridifferentiation and flexibility, high frequency letter is carried out to earthquake record in frequency domain Breath lifting, different frequency range overlapping reaches and improves seismic facies to Resolution of Acoustic Impedance Section, the mesh of effective identification thin interbed to expand frequency range 's.

To achieve these goals, it is pre- this application provides a kind of earthquake Relative Wave Impedance based on frequency domain pluridifferentiation The method of survey, the method includes：

The alternate structure thin interbed model in stratum is preset according to the first default stratum and second, and according to the thin interbed model Determine synthetic seismogram；

The synthetic seismogram is carried out into Fourier transformation and is transformed into frequency domain, and in frequency domain to it is described synthetically Shake record carries out the differential process of varying number rank, determines the earthquake record amplitude spectrum of the varying number rank；

Determine the relation fitting between the dominant frequency of the varying number rank and the varying number rank earthquake record amplitude spectrum Curve；

Equal difference frequency values of predetermined number and corresponding with the equal difference frequency values are determined according to the relation matched curve Differential order, and the differential process of the differential order are done to the earthquake record amplitude spectrum, obtain equal difference dominant frequency earthquake note Record amplitude spectrum；

The equal difference dominant frequency earthquake record amplitude spectrum is normalized, and to the amplitude after the normalized Spectrum be added and obtains overlapping amplitude spectrum；

Determine the envelope curve of the overlapping amplitude spectrum, and obtain on the envelope curve corresponding to the equal difference dominant frequency value Stacked spectrum envelope curve value；

The well logging sound wave impedance amplitude spectrum of the thin interbed model is obtained, the bag of the well logging sound wave impedance amplitude spectrum is determined Network curve, and obtain the well logging sound wave impedance amplitude envelope curve value on the envelope curve corresponding to the equal difference dominant frequency value；

Using the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values to the normalizing Amplitude spectrum after change treatment is corrected treatment and obtains correction overlapping amplitude spectrum；

Three rank differential phases spectrum according to the correction overlapping amplitude spectrum and the earthquake record determines wide frequency-band spectrum, and Earthquake Relative Wave Impedance is obtained according to the wide frequency-band spectrum.

In a preferred embodiment, the method also includes：It is micro- using varying number rank described in default compacting function pair The earthquake record amplitude spectrum for dividing carries out high frequency denoising.

In a preferred embodiment, the default compacting function is as follows：

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

In a preferred embodiment, the method also includes：Thin interbed is recognized according to the earthquake Relative Wave Impedance.

In a preferred embodiment, it is described that the alternate structure thin interbed in stratum is preset according to the first default stratum and second Model includes：

Interval velocity, density and the thickness change on the described first default stratum and the second default stratum are determined respectively, and Described first default stratum and the second default stratum are reduced and alternate structure from shallow-layer and deep layer to intermediate layer thickness equal difference Build the thin interbed model.

In a preferred embodiment, it is described to determine that synthetic seismogram includes according to the thin interbed model：

The wave impedance value of wave impedance value and the second default stratum according to the described first default stratum determines described thin The stratum reflectance factor of mutual layer model；

Convolution treatment is carried out to the stratum reflectance factor and default seismic wavelet and obtains synthetic seismogram.

In a preferred embodiment, the default seismic wavelet is Ricker wavelet, and the Ricker wavelet is zero phase The wavelet that position, perdurabgility are 80ms, dominant frequency is 50Hz.

In a preferred embodiment, it is described that varying number rank is carried out to the synthetic seismogram in frequency domain Differential process, determining the earthquake record amplitude spectrum of the varying number rank includes：

The synthetic seismogram is carried out at the differential of varying number rank using default differential process formula in frequency domain Reason, obtains the earthquake record amplitude spectrum of the varying number rank, and the default differential process formula is as follows：

Wherein, N represents the sampling number of the earthquake record, and S (f) is the frequency spectrum of the synthetic seismogram, and f represents institute State earthquake record frequency, Δ f represent the earthquake record frequency sampling interval (Δ t represents the earthquake note The time sampling interval of record), m is differential order, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record.

In a preferred embodiment, the determination varying number rank and the varying number rank earthquake record are shaken Relation matched curve between the dominant frequency of width spectrum includes：

The varying number rank and the varying number rank earthquake record amplitude spectrum are determined according to default dominant frequency computing formula Dominant frequency between relation discrete point；

Acquisition relation matched curve is fitted to the relation discrete point using least square method.

In a preferred embodiment, the dominant frequency computing formula is as follows：

Wherein, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record, f represents the earthquake record frequency, f_highRepresent the highest effective frequency of the amplitude spectrum of the m rank differential of the earthquake record, f_mRepresent the m of the earthquake record The dominant frequency of the amplitude spectrum of rank differential.

In a preferred embodiment, it is described to be shaken using the stacked spectrum envelope curve value and the well logging sound wave impedance Width spectrum envelope curve values are corrected treatment acquisition correction overlapping amplitude spectrum to the amplitude spectrum after the normalized to be included：

According to the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values, and combine default Iterative formula obtains correction coefficient；

The amplitude spectrum after the normalized is corrected according to the correction coefficient, and combination default updating formula Obtain correction overlapping amplitude spectrum.

In a preferred embodiment, the default iterative formula is as follows：

Wherein,When representing iterations as n, i-th correction coefficient of dominant frequency；When representing iterations as n+1, I-th correction coefficient of dominant frequency；N represents iterations (n >=1), and mf (i) represents i-th equal difference dominant frequency； (b_mf(i)When () represents dominant frequency as mf (i) f, corresponding well logging sound wave impedance amplitude spectrum envelope curve values)；I=1,2 ..., l, l generation Table dominant frequency quantity；cⁿIt is mf (i) that (mf (i)) represents dominant frequency, and iterations is when being n, corresponding stacked spectrum envelope curve value.

In a preferred embodiment, the default updating formula is as follows：

Wherein, FⁿOverlapping amplitude spectrum when representing generation number as n；TheWhen representing iterations as n, i-th dominant frequency Correction coefficient；|F_mf(i)| normalization differential amplitude spectrum when representing dominant frequency as mf (i).

On the other hand the application also provides a kind of earthquake Relative Wave Impedance prediction meanss based on frequency domain pluridifferentiation, should Device includes：

Earthquake record determining unit, for being preset according to the first default stratum and second, stratum is alternate to build thin interbed mould Type, and synthetic seismogram is determined according to the thin interbed model；

Earthquake record amplitude spectrum determining unit, frequency is transformed into for the synthetic seismogram to be carried out into Fourier transformation Domain, and the differential process of varying number rank are carried out to the synthetic seismogram in frequency domain, determine the varying number rank Earthquake record amplitude spectrum；

Preset relation matched curve determining unit, for determining the varying number rank and varying number terrace shake note Record the relation matched curve between the dominant frequency of amplitude spectrum；

Equal difference dominant frequency earthquake record amplitude spectrum obtaining unit, for determining predetermined number according to the relation matched curve Equal difference frequency values and differential order corresponding with the equal difference frequency values, and the differential is done to the earthquake record amplitude spectrum The differential process of exponent number, obtain equal difference dominant frequency earthquake record amplitude spectrum；

Overlapping amplitude spectrum determining unit, for being normalized to the equal difference dominant frequency earthquake record amplitude spectrum, and Amplitude spectrum after the normalized be added and obtains overlapping amplitude spectrum；

First envelope curve value determining unit, for determining the envelope curve of the overlapping amplitude spectrum, and obtains the bag Stacked spectrum envelope curve value on network curve corresponding to the equal difference dominant frequency value；

Second envelope curve value determining unit, the well logging sound wave impedance amplitude for obtaining the thin interbed model is composed, really The envelope curve of the fixed well logging sound wave impedance amplitude spectrum, and obtain on the envelope curve corresponding to the equal difference dominant frequency value Well logging sound wave impedance amplitude envelope curve value；

Correction unit, for utilizing the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve Value is corrected treatment to the amplitude spectrum after the normalized and obtains correction overlapping amplitude spectrum；

Earthquake Relative Wave Impedance obtaining unit, for three ranks according to the correction overlapping amplitude spectrum and the earthquake record Differential phase spectrum determines wide frequency-band spectrum, and obtains earthquake Relative Wave Impedance according to the wide frequency-band spectrum.

In a preferred embodiment, the device also includes：High frequency denoising unit, for using default compacting letter Several earthquake record amplitude spectrums to the varying number rank differential carry out high frequency denoising.

In a preferred embodiment, the default compacting function is as follows：

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

In a preferred embodiment, the device also includes：Recognition unit, for according to the earthquake Relative Wave Impedance Identification thin interbed.

As can be seen here, the technical scheme that the embodiment of the present application is provided by setting up thin interbed model, and according to it is described it is thin mutually Layer model obtains earthquake record.Then, carrying out not same order differential to synthetic seismogram in frequency domain carries out high frequency lifting, and The relation set up between amplitude spectrum differential order and dominant frequency, is overlapped by the addition of earthquake record not same order differential amplitude spectrum Amplitude spectrum, and the well logging demarcation carried out to earthquake record pluridifferentiation amplitude spectrum by well logging sound wave impedance, after being corrected Overlapping amplitude spectrum, the phase spectrum in conjunction with the rank differential of earthquake record three obtains wide frequency-band spectrum, and is obtained according to wide frequency-band spectrum High-resolution seismic exploration Relative Wave Impedance, such that it is able to effectively identify thin interbed.Compared with prior art, it is possible to achieve earthquake Record the high-frequency information lifting of amplitude spectrum, and different frequency range overlapping to expand frequency range, effectively raise earthquake Relative Wave Impedance Resolution ratio, increased the recognition capability of thin interbed.

Brief description of the drawings

In order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor Under, other accompanying drawings can also be obtained according to these accompanying drawings.

Fig. 1 is a kind of stream of the embodiment of the earthquake Relative Wave Impedance Forecasting Methodology based on frequency domain pluridifferentiation of the application Cheng Tu；

Fig. 2 is thin sand-mud interbed model schematic；

Fig. 3 is earthquake record described in the embodiment of the present application and 1 rank of the earthquake record, 2 ranks, 3 ranks and 4 rank differential amplitudes Spectrum schematic diagram；

Fig. 4 is the relation curve schematic diagram of differential order and corresponding rank differential earthquake record dominant frequency of amplitude spectrum；

Fig. 5 is the schematic diagram of the envelope curve that the embodiment of the present application determines overlapping amplitude spectrum；

Fig. 6 is the schematic diagram of the amplitude spectrum with the amplitude spectrum of well logging sound wave impedance of the stratum reflectance factor of thin interbed model；

Fig. 7 is the correction overlapping amplitude spectrum schematic diagram after the amplitude spectrum and 3 iteration of earthquake record in the embodiment of the present application；

Fig. 8 is the correction overlapping after 3 iteration in well logging sound wave impedance amplitude spectrum and Fig. 7 in Fig. 6 shown in b shown in b The fractional amplitude joint and fitting function schematic diagram of amplitude spectrum；

Fig. 9 is the trace integral and the seismic facies pair obtained using technical scheme of earthquake record in the embodiment of the present application Wave impedance records schematic diagram；

Figure 10 is a kind of earthquake Relative Wave Impedance prediction dress based on frequency domain pluridifferentiation that the embodiment of the present application is provided The schematic diagram put.

Specific embodiment

In order that those skilled in the art more fully understand the technical scheme in the application, below in conjunction with the application reality The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation Example is only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection Scope.

Below implementing for the embodiment of the present application is described in detail with several specific examples.

Hereinafter a kind of reality of the earthquake Relative Wave Impedance Forecasting Methodology based on frequency domain pluridifferentiation of the application is introduced first Apply example.With reference to accompanying drawing 1, the embodiment includes：

S110：The alternate structure thin interbed model in stratum is preset according to the first default stratum and second, and according to described thin mutual Layer model determines synthetic seismogram.

In some embodiments it is possible to being preset according to the first default stratum and second, stratum is alternate to build thin interbed model, And synthetic seismogram is determined according to the thin interbed model.

In certain embodiments, presetting the alternate structure thin interbed model in stratum here according to the first default stratum and second can To include：Interval velocity, density and the thickness change on the described first default stratum and the second default stratum are determined respectively, and will Described first default stratum and the second default stratum are reduced and alternate structure from shallow-layer and deep layer to intermediate layer thickness equal difference The thin interbed model.Here the default stratum of first can be mud stone, and the described second default stratum can be sandstone.Specifically , for example, the described first default stratum and the second default stratum are reduced from shallow-layer and deep layer to intermediate layer thickness equal difference And alternate composition sets up the thin interbed model.With reference to accompanying drawing 2, the public affairs of two neighboring layer thickness in thin interbed model here Difference can be 2m, namely thickness become turn to 2m；28m, 26m ..., 2m, 4m ... can be followed successively by from shallow-layer to deep layer thickness, 40m, dark stratum can be mud stone in figure, and the stratum of light color can be sandstone, and mud stone and mud stone interval velocity can be with identical, examples 4000m/s is such as, the density of mud stone can be 3g/cm³, the density of sandstone can be 1.5g/cm³。

Further, the interval velocity on the described first default stratum can be multiplied with the density on the described first default stratum and is obtained Take the wave impedance value on the described first default stratum；And the interval velocity and described second on the described second default stratum is preset into stratum Density is multiplied and obtains the wave impedance value on the described second default stratum.

In certain embodiments, determine that synthetic seismogram can include here according to the thin interbed model：

The wave impedance value of wave impedance value and the second default stratum according to the described first default stratum determines described thin The stratum reflectance factor of mutual layer model.

Specifically, it is possible to use the computing formula of following stratum reflectance factor：

Wherein, r represents the stratum reflectance factor of the thin interbed model；w₁Represent the wave impedance value on the first default stratum；w₂ Represent the wave impedance value on the second default stratum.

Here seismic wavelet can preset, specifically, seismic wavelet can be Ricker wavelet, and the thunder here It zero phase, perdurabgility is that 80ms, dominant frequency are the wavelet of 50Hz that gram wavelet is.

Further, above-mentioned Ricker wavelet and the stratum reflectance factor of above-mentioned thin interbed model are carried out into convolution treatment can be with Obtain synthetic seismogram.

S120：The synthetic seismogram is carried out into Fourier transformation and is transformed into frequency domain, and to described in frequency domain Synthetic seismogram carries out the differential process of varying number rank, determines the earthquake record amplitude spectrum of the varying number rank.

In certain embodiments, after step S110 obtains synthetic seismogram, can be by the synthetic seismogram Carry out Fourier transformation and be transformed into frequency domain, and carry out the differential of varying number rank to the synthetic seismogram in frequency domain Treatment, determines the earthquake record amplitude spectrum of the varying number rank.

In some embodiments it is possible to be carried out at the differential of varying number rank to the synthetic seismogram in frequency domain Reason, determines the earthquake record amplitude spectrum of the varying number rank.Specifically, can be using default differential process formula in frequency domain The differential process of varying number rank are carried out to the earthquake record, the earthquake record amplitude spectrum of the varying number rank, institute is obtained State default differential process formula as follows：

Wherein, N is the sampling number of the earthquake record, and S (f) is the frequency spectrum of the synthetic seismogram, and f is frequency, Δ f be frequency sampling interval (Δ t is the time sampling interval of the earthquake record), m is differential order, F_m It is the amplitude spectrum of the m rank differential of the earthquake record.

It is illustrated in figure 3 the earthquake record and 1 rank of the earthquake record, 2 ranks, 3 ranks and 4 rank differential amplitude spectrums.From Fig. 3 is visible, and with the increase of differential order, its resolution ratio is significantly improved, and middle part thin layer can be highlighted gradually.First differential Phase 90 degree than synthetic seismogram phase difference, and often increase by 2 ranks, phase place change 180 degree.

Further, in certain embodiments, after the earthquake record amplitude spectrum for obtaining different predetermined number ranks, can be with High frequency denoising is carried out using the earthquake record amplitude spectrum of default compacting function pair difference predetermined number rank differential.Specifically, The earthquake record amplitude spectrum of each predetermined number rank can be multiplied with the compacting function.The default compacting function is as follows：

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

S130：Determine the relation between the dominant frequency of the varying number rank and the varying number rank earthquake record amplitude spectrum Matched curve.

In certain embodiments, after step S120, it may be determined that the varying number rank and the varying number rank Relation matched curve between the dominant frequency of earthquake record amplitude spectrum.Specifically, can include：It is true according to default dominant frequency computing formula Fixed relation discrete point between the varying number rank and the dominant frequency of the varying number rank earthquake record amplitude spectrum, and using most Small square law is fitted acquisition relation matched curve to the relation discrete point.

Here varying number rank earthquake record dominant frequency of amplitude spectrum namely the earthquake record amplitude spectrum effective frequency section in The corresponding frequency of centre of area point of frequency axis direction, here effective frequency section can according to collection dynamic range be set It is fixed, therefore dominant frequency computing formula can be with as follows：

Wherein, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record, f represents the earthquake record frequency, f_highRepresent the highest effective frequency of the amplitude spectrum of the m rank differential of the earthquake record, f_mRepresent the m of the earthquake record The dominant frequency of the amplitude spectrum of rank differential, namely the earthquake record m rank differential amplitude spectrum center dot frequency.

Further, the varying number rank and described each predetermined number rank earthquake record are shaken using least square method Relation discrete point between the dominant frequency of width spectrum is fitted acquisition relation matched curve, if Fig. 4 is the relation matched curve Partial schematic diagram, namely differential order and corresponding rank differential earthquake record dominant frequency of amplitude spectrum relation curve schematic diagram.

S140：According to the relation matched curve determine predetermined number equal difference frequency values and with the equal difference frequency values Corresponding differential order, and the differential process of the differential order are done to the earthquake record amplitude spectrum, obtain equal difference dominant frequency ground Shake record amplitude spectrum.

In certain embodiments, differential order and corresponding rank differential earthquake record dominant frequency of amplitude spectrum are obtained in step S130 After relation curve, the frequency values of the equal difference of predetermined number can be determined according to the curve, and determine and the frequency values pair The differential order answered, the frequency values quantity of the equal difference chosen here can as the case may be carry out default setting.Thus may be used Equal difference dominant frequency relation is converted to by differential order, such that it is able to control the scope of dominant frequency value and frequency bandwidth, because simple Integer rank differentiation result is added several times, is likely to result in the exceptional value of a certain frequency range, and in order to avoid such case, we are controlled Amount processed is converted into equal difference dominant frequency by differential order, is so shaken come Earthquake occurrence control record using the differential order corresponding to equal difference dominant frequency The scope and dominant frequency value of width spectrum, can preferably constrain the change of amplitude spectrum.

Further, the differential process that the differential order can be done to the earthquake record amplitude spectrum obtain equal difference dominant frequency Earthquake record amplitude spectrum.

S150：The equal difference dominant frequency earthquake record amplitude spectrum is normalized, and to the normalized after Amplitude spectrum carry out be added obtain overlapping amplitude spectrum.

In some embodiments it is possible to be normalized to the dominant frequency amplitude spectrum of the equal difference dominant frequency earthquake record amplitude spectrum Treatment, so ensures that the area of not same order differential amplitude spectrum is identical；And the amplitude spectrum after the normalized is entered Row is added and obtains overlapping amplitude spectrum, specifically, overlapping amplitude spectrum can be determined with formula is calculated as below：

Wherein mf (i) represents i-th equal difference dominant frequency, | F_mf(i)| the normalization of earthquake record when representing dominant frequency as mf (i) Differential amplitude spectrum, F represents overlapping amplitude spectrum；

S160：Determine the envelope curve of the overlapping amplitude spectrum, and obtain the equal difference dominant frequency value on the envelope curve Corresponding stacked spectrum envelope curve value.

In certain embodiments, with reference to accompanying drawing 5, determining the envelope curve of the overlapping amplitude spectrum can include：

S510：The local maximum on the overlapping amplitude spectrum is obtained successively, and the local maximum shakes for the overlapping Amplitude of the upper amplitude of width spectrum more than two neighboring amplitude.

S520：Local minimum is obtained in the local maximum, the local minimum is the local maximum Amplitude of the middle numerical value less than two neighboring amplitude.

S530：The local minimum is removed, and it is corresponding to be sequentially connected the amplitude adjacent with the local minimum Point obtains the envelope curve of the overlapping amplitude spectrum.

Further, it is determined that after the envelope curve of the overlapping amplitude spectrum, can obtain every on the envelope curve The corresponding stacked spectrum envelope curve value of individual equal difference dominant frequency.

S170：The well logging sound wave impedance amplitude spectrum of the thin interbed model is obtained, the well logging sound wave impedance amplitude is determined The envelope curve of spectrum, and it is bent to obtain the well logging sound wave impedance amplitude envelope on the envelope curve corresponding to the equal difference dominant frequency value Line value.

In certain embodiments, the well logging sound wave impedance amplitude spectrum that the thin interbed model is obtained here can include：Will The well logging sound wave impedance transformation to frequency domain is obtained the well logging by the wave impedance of thin interbed model as well logging sound wave impedance Sound impedance amplitude spectrum.The wave impedance of thin interbed model described here includes the wave impedance on the described first default stratum with described the The wave impedance on two default stratum.

In certain embodiments, determining the envelope curve of the well logging sound wave impedance amplitude spectrum can include：Obtain successively Local maximum in the well logging sound wave impedance amplitude spectrum, the local maximum is in well logging sound wave impedance amplitude spectrum Amplitude of the amplitude more than two neighboring amplitude；Local minimum, the local pole are obtained in the local maximum Small value is amplitude of the numerical value less than two neighboring amplitude in the local maximum；Remove the local minimum, and according to The envelope that the corresponding point of the secondary connection amplitude adjacent with the local minimum obtains the well logging sound wave impedance amplitude spectrum is bent Line.

Further, the well logging sound wave impedance amplitude corresponding to the equal difference dominant frequency value on the envelope curve can be obtained Envelope curve value.

S180：Using the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values to described Amplitude spectrum after normalized is corrected treatment and obtains correction overlapping amplitude spectrum.

In some embodiments it is possible to using the stacked spectrum envelope curve value and well logging sound wave impedance amplitude spectrum bag Network curve values are corrected treatment to the amplitude spectrum after the normalized and obtain correction overlapping amplitude spectrum.As shown in Figure 6 It is the schematic diagram of the amplitude spectrum of the amplitude spectrum and well logging sound wave impedance of the stratum reflectance factor of thin interbed model.In Fig. 6 shown in a It is the schematic diagram of the amplitude spectrum of the stratum reflectance factor of thin interbed model；In Fig. 6 shown in b be well logging sound wave impedance amplitude spectrum Schematic diagram, namely the amplitude spectrum of thin interbed model wave impedance schematic diagram.As seen from Figure 6, the one of the stratum reflectance factor Rank integration is wave impedance, namely well logging sound wave impedance amplitude spectrum is a kind of weighting compacting to stratum reflectance factor amplitude spectrum, two The local frequencies of person are constant.The application is linearly lifted radio-frequency component by carrying out pluridifferentiation to earthquake record in frequency domain Rise, very little is influenceed on local frequency distribution, that is, eliminate the influence of seismic wavelet, and according to the stratum reflectance factor amplitude spectrum Spectrum, can directly using earthquake record amplitude spectrum of the acoustic impedance to varying number rank of logging well with the relation between wave impedance amplitude spectrum It is corrected, so as to obtain earthquake Relative Wave Impedance.Specifically, can include：

According to the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values, and combine default Iterative formula obtains correction coefficient.

In certain embodiments, correction coefficient here can be according to the stacked spectrum envelope curve value and the well logging sound Wave impedance amplitude spectrum envelope curve value, and default iterative formula acquisition is combined, specifically, the iteration coefficient can be included for the first time Iteration coefficient and non-first iteration coefficient, the first iteration coefficient namely well logging sound wave impedance amplitude spectrum envelope curve values, institute Stating non-first iteration coefficient can be obtained by following default iterative formula：

Wherein,When representing iterations as n, i-th correction coefficient of dominant frequency；When representing iterations as n+1, I-th correction coefficient of dominant frequency；N represents iterations (n >=1), and mf (i) represents i-th equal difference dominant frequency； (b_mf(i)When () represents dominant frequency as mf (i) f, corresponding well logging sound wave impedance amplitude spectrum envelope curve values)；I=1,2 ..., l, l generation Table dominant frequency quantity；cⁿIt is mf (i) that (mf (i)) represents dominant frequency, and iterations is when being n, corresponding stacked spectrum envelope curve value.

The amplitude spectrum after the normalized is corrected according to the correction coefficient, and combination default updating formula Obtain correction overlapping amplitude spectrum.

In certain embodiments, here default updating formula is as follows：

Wherein, FⁿOverlapping amplitude spectrum when representing generation number as n；TheWhen representing iterations as n, i-th dominant frequency Correction coefficient；|F_mf(i)| normalization differential amplitude spectrum when representing dominant frequency as mf (i).

Further, after being corrected to the overlapping amplitude spectrum using the correction coefficient after each iteration, can be with root According to the new stacked spectrum envelope curve value of overlapping amplitude spectrum after current correction, and obtain new according to current stacked spectrum envelope curve value Iteration coefficient, recycles current correction factors be corrected to current amplitude spectrum to the good school of acquisition high boost effect Positive overlapping amplitude spectrum.

In certain embodiments, shown in b it is 3 iteration in Fig. 7 by 3 above-mentioned iterative process with reference to accompanying drawing 7 Afterwards correction overlapping amplitude spectrum, in Fig. 7 shown in a be synthetic seismogram amplitude spectrum, it is seen that by being obtained after iteration High boost effect corrects overlapping amplitude spectrum well.

Fig. 8 is the correction overlapping after 3 iteration in well logging sound wave impedance amplitude spectrum and Fig. 7 in Fig. 6 shown in b shown in b The fractional amplitude joint and fitting function schematic diagram of amplitude spectrum, i.e. this amplitude are well logging acoustic impedance values and correction in 0 to 1.2 The joint of amplitude spectrum and the schematic diagram of the fitting function obtained according to the joint are overlapped, and fitting function is y= 2.0105x-0.021, degree of fitting has reached 0.9528, and wherein x represents well logging sound wave impedance amplitude spectrum, after y represents 3 iteration Correction overlapping amplitude spectrum.

S190：Three rank differential phases spectrum according to the correction overlapping amplitude spectrum and the earthquake record determines broadband frequently Spectrum, and earthquake Relative Wave Impedance is obtained according to the wide frequency-band spectrum.

In certain embodiments, the phase spectrum phase shift of the rank differential of earthquake record three is -90 °, its phase and trace integral phase It is identical, so identical with well logging sound wave impedance phase.Therefore, the phase of overlapping amplitude spectrum and the rank differential of earthquake record three will be corrected Spectrum joint obtains wide frequency-band spectrum, and obtains earthquake Relative Wave Impedance according to the wide frequency-band spectrum.

As shown in figure 9, wherein, in Fig. 9 shown in a be the earthquake record trace integral schematic diagram, in Fig. 9 shown in b Be using technical scheme obtain earthquake Relative Wave Impedance record figure.It can be seen that the road of the earthquake record Integration due to the missing of radio-frequency component, cause thin layer fail well identification, and process obtain Relative Wave Impedance record due to Radio-frequency component is improved, is had for the thin layer at model center position and is preferably highlighted effect, in theoretical model, 2m can be recognized Reservoir.

As can be seen here, the technical scheme that the embodiment of the present application is provided by setting up thin interbed model, and according to it is described it is thin mutually Layer model obtains synthetic seismogram.Then, carrying out not same order differential to synthetic seismogram in frequency domain carries out high frequency lift Rise, and the relation set up between amplitude spectrum differential order and dominant frequency, obtained by the addition of earthquake record not same order differential amplitude spectrum Amplitude spectrum spectrum must be overlapped, and the well logging carried out to earthquake record pluridifferentiation amplitude spectrum by well logging sound wave impedance is demarcated, and is obtained Overlapping amplitude spectrum after correction, the phase spectrum in conjunction with the rank differential of earthquake record three obtains wide frequency-band spectrum, and according to broadband The high-resolution seismic exploration Relative Wave Impedance record figure that frequency spectrum is obtained, such that it is able to effectively identify thin interbed.With prior art phase Than, it is possible to achieve the high-frequency information of earthquake record amplitude spectrum is lifted, and different frequency range overlaps to expand frequency range, effectively raises The resolution ratio of earthquake Relative Wave Impedance, increased the discrimination of thin interbed.

On the other hand the application also provides a kind of earthquake Relative Wave Impedance prediction meanss based on frequency domain pluridifferentiation, ties Accompanying drawing 10 is closed, this fills 1000 includes：

Earthquake record determining unit 1010, for according to the first default stratum and the second default alternate structure thin interbed in stratum Model, and synthetic seismogram is determined according to the thin interbed model；

Earthquake record amplitude spectrum determining unit 1020, is transformed into for the synthetic seismogram to be carried out into Fourier transformation Frequency domain, and the differential process of varying number rank are carried out to the synthetic seismogram in frequency domain, determine the different numbers Measure the earthquake record amplitude spectrum of rank；

Preset relation matched curve determining unit 1030, for determining the varying number rank and the varying number terrace Relation matched curve between the dominant frequency of shake record amplitude spectrum；

Equal difference dominant frequency earthquake record amplitude spectrum obtaining unit 1040, for determining present count according to the relation matched curve The equal difference frequency values of amount and differential order corresponding with the equal difference frequency values, and do described to the earthquake record amplitude spectrum The differential process of differential order, obtain equal difference dominant frequency earthquake record amplitude spectrum；

Overlapping amplitude spectrum determining unit 1050, for being normalized place to the equal difference dominant frequency earthquake record amplitude spectrum Reason, and the amplitude spectrum after the normalized is carried out to be added acquisition overlapping amplitude spectrum；

First envelope curve value determining unit 1060, for determining the envelope curve of the overlapping amplitude spectrum, and obtains institute State the stacked spectrum envelope curve value corresponding to the equal difference dominant frequency value on envelope curve；

Second envelope curve value determining unit 1070, the well logging sound wave impedance amplitude for obtaining the thin interbed model Spectrum, determines the envelope curve of the well logging sound wave impedance amplitude spectrum, and obtains the equal difference dominant frequency value institute on the envelope curve Corresponding well logging sound wave impedance amplitude envelope curve value；

Correction unit 1080, for utilizing the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope Curve values are corrected treatment to the amplitude spectrum after the normalized and obtain correction overlapping amplitude spectrum；

Earthquake Relative Wave Impedance record figure obtaining unit 1090, for according to correction overlapping amplitude spectrum and the earthquake The three rank differential phases spectrum of record determines wide frequency-band spectrum, and obtains earthquake Relative Wave Impedance according to the wide frequency-band spectrum.

In a preferred embodiment, the device 1000 also includes：High frequency denoising unit, for using default pressure The earthquake record amplitude spectrum of varying number rank differential carries out high frequency denoising described in function pair processed.

In a preferred embodiment, the default compacting function is as follows：

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

In a preferred embodiment, the device 1000 also includes：Recognition unit, for according to the seismic facies to ripple Impedance recognizes thin interbed.

In a preferred embodiment, it is described that the alternate structure thin interbed in stratum is preset according to the first default stratum and second Model includes：

Determine respectively the interval velocity on the described first default stratum and the second default stratum, density and with thickness change, And reduce and alternate described first default stratum and the second default stratum to intermediate layer thickness equal difference from shallow-layer and deep layer Build the thin interbed model.

In a preferred embodiment, it is described to determine that synthetic seismogram includes according to the thin interbed model：

The wave impedance value of wave impedance value and the second default stratum according to the described first default stratum determines described thin The stratum reflectance factor of mutual layer model；

Convolution treatment is carried out to the stratum reflectance factor and default seismic wavelet and obtains synthetic seismogram.

In a preferred embodiment, the default seismic wavelet is Ricker wavelet, and the Ricker wavelet is zero phase The wavelet that position, perdurabgility are 80ms, dominant frequency is 50Hz.

In a preferred embodiment, it is described that the earthquake record amplitude spectrum is carried out at the differential of varying number rank Reason, determining the earthquake record amplitude spectrum of the varying number rank includes：

The earthquake record amplitude spectrum of the synthetic seismogram is carried out not using default differential process formula in frequency domain With the differential process of quantity rank, the earthquake record amplitude spectrum of the varying number rank is obtained, the default differential process formula is such as Under：

Wherein, N represents the sampling number of the earthquake record, and S (f) is the frequency spectrum of the synthetic seismogram, and f represents institute State earthquake record frequency, Δ f represent the earthquake record frequency sampling interval (Δ t represents the earthquake note The time sampling interval of record), m is differential order, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record.

In a preferred embodiment, the determination varying number rank and the varying number rank earthquake record are shaken Relation matched curve between the dominant frequency of width spectrum includes：

The varying number rank and the varying number rank earthquake record amplitude spectrum are determined according to default dominant frequency computing formula Dominant frequency between relation discrete point；

Acquisition relation matched curve is fitted to the relation discrete point using least square method.

In a preferred embodiment, the dominant frequency computing formula is as follows：

Wherein, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record, f represents the earthquake record frequency, f_highRepresent the highest effective frequency of the amplitude spectrum of the m rank differential of the earthquake record, f_mRepresent the m of the earthquake record The dominant frequency of the amplitude spectrum of rank differential.

In a preferred embodiment, it is described to be shaken using the stacked spectrum envelope curve value and the well logging sound wave impedance Width spectrum envelope curve values are corrected treatment acquisition correction overlapping amplitude spectrum to the amplitude spectrum after the normalized to be included：

According to the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values, and combine default Iterative formula obtains correction coefficient；

The amplitude spectrum after the normalized is corrected according to the correction coefficient, and combination default updating formula Obtain correction overlapping amplitude spectrum.

In a preferred embodiment, the default iterative formula is as follows：

Wherein,When representing iterations as n, i-th correction coefficient of dominant frequency；When representing iterations as n+1, I-th correction coefficient of dominant frequency；N represents iterations (n >=1), and mf (i) represents i-th equal difference dominant frequency； (b_mf(i)When (ω) represents dominant frequency as mf (i), corresponding well logging sound wave impedance amplitude spectrum envelope curve values)；I=1,2 ..., l, l Represent dominant frequency quantity；cⁿIt is mf (i) that (mf (i)) represents dominant frequency, and iterations is when being n, corresponding stacked spectrum envelope curve value.

In a preferred embodiment, the default updating formula is as follows：

Wherein, FⁿOverlapping amplitude spectrum when representing generation number as n；TheWhen representing iterations as n, i-th dominant frequency Correction coefficient；|F_mf(i)| normalization differential amplitude spectrum when representing dominant frequency as mf (i).

By above the application a kind of earthquake Relative Wave Impedance Forecasting Methodology and the reality of device based on frequency domain pluridifferentiation The technical scheme for applying example offer is visible, and the embodiment of the present application is obtained by setting up thin interbed model according to the thin interbed model Take synthetic seismogram.Then, carrying out not same order differential to synthetic seismogram in frequency domain carries out high frequency lifting, and sets up Relation between amplitude spectrum differential order and dominant frequency, by earthquake record not same order differential amplitude spectrum addition obtain overlapping amplitude Spectrum, and the well logging demarcation carried out to earthquake record pluridifferentiation amplitude spectrum by well logging sound wave impedance, the overlapping after being corrected Amplitude spectrum, the phase spectrum in conjunction with the rank differential of earthquake record three obtains wide frequency-band spectrum, and the height obtained according to wide frequency-band spectrum Resolution seismic Relative Wave Impedance, such that it is able to effectively identify thin interbed.Compared with prior art, it is possible to achieve earthquake record The high-frequency information of amplitude spectrum is lifted, and different frequency range overlaps to expand frequency range, effectively raises dividing for earthquake Relative Wave Impedance Resolution, increased the discrimination of thin interbed.

Each embodiment in this specification is described by the way of progressive, what each embodiment was stressed be with The difference of other embodiment, between each embodiment identical similar part mutually referring to.Especially for system reality Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.

Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application's Spirit.

Claims (17)

1. a kind of earthquake Relative Wave Impedance Forecasting Methodology based on frequency domain pluridifferentiation, it is characterised in that the method includes：

The alternate structure thin interbed model in stratum is preset according to the first default stratum and second, and is determined according to the thin interbed model Synthetic seismogram；

The synthetic seismogram is carried out into Fourier transformation and is transformed into frequency domain, and to the synthesis earthquake note in frequency domain Record carries out the differential process of varying number rank, determines the earthquake record amplitude spectrum of the varying number rank；

Determine the relation matched curve between the dominant frequency of the varying number rank and the varying number rank earthquake record amplitude spectrum；

Equal difference frequency values of predetermined number and corresponding micro- with the equal difference frequency values are determined according to the relation matched curve Sublevel number, and the differential process of the differential order are done to the earthquake record amplitude spectrum, obtain equal difference dominant frequency earthquake record and shake Width is composed；

The equal difference dominant frequency earthquake record amplitude spectrum is normalized, and the amplitude spectrum after the normalized is entered Row is added and obtains overlapping amplitude spectrum；

Determine the envelope curve of the overlapping amplitude spectrum, and obtain folded corresponding to the equal difference dominant frequency value on the envelope curve Close spectrum envelope curve values；

The well logging sound wave impedance amplitude spectrum of the thin interbed model is obtained, determines that the envelope of the well logging sound wave impedance amplitude spectrum is bent Line, and obtain the well logging sound wave impedance amplitude envelope curve value on the envelope curve corresponding to the equal difference dominant frequency value；

Using the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values to the normalization at Amplitude spectrum after reason is corrected treatment and obtains correction overlapping amplitude spectrum；

According to the correction overlapping amplitude spectrum and the three rank differential phases spectrum determination wide frequency-band spectrum of the earthquake record, and according to The wide frequency-band spectrum obtains earthquake Relative Wave Impedance.

2. method according to claim 1, it is characterised in that the method also includes：Using described in default compacting function pair The earthquake record amplitude spectrum of varying number rank differential carries out high frequency denoising.

3. method according to claim 2, it is characterised in that the default compacting function is as follows：

$T=\left\{\begin{array}{cc}{e}^{-a(f-f_{high})}& f>f_{high}\\ 1& f\≤f_{high}\end{array}\right.$

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

4. according to any described method of claims 1 to 3, it is characterised in that the method also includes：According to the seismic facies pair Wave impedance recognizes thin interbed.

5. according to any described method of claims 1 to 3, it is characterised in that described pre- according to the first default stratum and second If the alternate thin interbed model that builds in stratum includes：

Determine interval velocity, density and the thickness change on the described first default stratum and the second default stratum respectively, and by institute The first default stratum and the second default stratum is stated to be reduced and alternate structure institute to intermediate layer thickness equal difference from shallow-layer and deep layer State thin interbed model.

6. according to any described method of claims 1 to 3, it is characterised in that described to be determined to close according to the thin interbed model Include into earthquake record：

The wave impedance value of wave impedance value and the second default stratum according to the described first default stratum determines the thin interbed The stratum reflectance factor of model；

Convolution treatment is carried out to the stratum reflectance factor and default seismic wavelet and obtains synthetic seismogram.

7. method according to claim 6, it is characterised in that the default seismic wavelet is Ricker wavelet, and the thunder It zero phase, perdurabgility is that 80ms, dominant frequency are the wavelet of 50Hz that gram wavelet is.

8. according to any described method of claims 1 to 3, it is characterised in that it is described in frequency domain to the synthesis earthquake Record carries out the differential process of varying number rank, and determining the earthquake record amplitude spectrum of the varying number rank includes：

Carry out the differential process of varying number rank to the synthetic seismogram using default differential process formula in frequency domain, obtain The earthquake record amplitude spectrum of the varying number rank is obtained, the default differential process formula is as follows：

$F_m\left(f\right)=\left|S\left(f\right)\right|\·{\[2\sin \left(\frac{\mathrm{\π}}{N}\frac{f}{\mathrm{\Δ}f}\right)\]}^m$

Wherein, N represents the sampling number of the earthquake record, and S (f) is the frequency spectrum of the synthetic seismogram, and f is represented describedly Shake recording frequency, Δ f represents the sampling interval of the earthquake record frequency,Δ t represents the earthquake record Time sampling interval, m is differential order, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record.

9. according to any described method of claims 1 to 3, it is characterised in that the determination varying number rank and described Relation matched curve between the dominant frequency of varying number rank earthquake record amplitude spectrum includes：

The master of the varying number rank and the varying number rank earthquake record amplitude spectrum is determined according to default dominant frequency computing formula Relation discrete point between frequency；

Acquisition relation matched curve is fitted to the relation discrete point using least square method.

10. method according to claim 9, it is characterised in that the dominant frequency computing formula is as follows：

$\underset{f=0}{\overset{f_m}{\mathrm{\Σ}}}{F}_m\left(f\right)=\frac{1}{2}\underset{f=0}{\overset{f_{high}}{\mathrm{\Σ}}}{F}_m\left(f\right)$

Wherein, F_mF () represents the amplitude spectrum of the m rank differential of the earthquake record, f represents the earthquake record frequency, f_highGeneration The highest effective frequency of the amplitude spectrum of the m rank differential of earthquake record described in table, f_mRepresent the m rank differential of the earthquake record Amplitude spectrum dominant frequency.

11. according to any described method of claims 1 to 3, it is characterised in that described to utilize the stacked spectrum envelope curve value Treatment is corrected to the amplitude spectrum after the normalized with the well logging sound wave impedance amplitude spectrum envelope curve values to obtain Correction overlapping amplitude spectrum includes：

According to the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values, and combine default iteration Formula obtains correction coefficient；

According to the correction coefficient, and the default updating formula of combination is corrected acquisition to the amplitude spectrum after the normalized Correction overlapping amplitude spectrum.

12. methods according to claim 11, it is characterised in that the default iterative formula is as follows：

$a_i^{n+1}=a_i^n\·b\left(mf\right(i\left)\right)/c^n\left(mf\right(i\left)\right)$

Wherein,When representing iterations as n, i-th correction coefficient of dominant frequency；When representing iterations as n+1, i-th The correction coefficient of dominant frequency；N represents iterations, and n >=1, and mf (i) represents i-th equal difference dominant frequency；b_mf(i) When () represents dominant frequency as mf (i) f, corresponding well logging sound wave impedance amplitude spectrum envelope curve values；I=1,2 ..., l, l represents dominant frequency Quantity；cⁿIt is mf (i) that (mf (i)) represents dominant frequency, and iterations is when being n, corresponding stacked spectrum envelope curve value.

13. methods according to claim 11, it is characterised in that the default updating formula is as follows：

$F^n=\underset{i=1}{\overset{l}{\Σ}}{a}_i^n\left|F_{mf\left(i\right)}\right|$

Wherein, FⁿOverlapping amplitude spectrum when representing generation number as n；TheWhen representing iterations as n, i-th correction system of dominant frequency Number；|F_mf(i)| normalization differential amplitude spectrum when representing dominant frequency as mf (i).

14. a kind of earthquake Relative Wave Impedance prediction meanss based on frequency domain pluridifferentiation, it is characterised in that the device includes：

Earthquake record determining unit, for presetting the alternate structure thin interbed model in stratum according to the first default stratum and second, and Synthetic seismogram is determined according to the thin interbed model；

Earthquake record amplitude spectrum determining unit, frequency domain is transformed into for the synthetic seismogram to be carried out into Fourier transformation, And the differential process of varying number rank are carried out to the synthetic seismogram in frequency domain, determine the ground of the varying number rank Shake record amplitude spectrum；

Preset relation matched curve determining unit, for determining that the varying number rank and the varying number rank earthquake record are shaken Relation matched curve between the dominant frequency of width spectrum；

Equal difference dominant frequency earthquake record amplitude spectrum obtaining unit, the equal difference for determining predetermined number according to the relation matched curve Frequency values and differential order corresponding with the equal difference frequency values, and the differential order is done to the earthquake record amplitude spectrum Differential process, obtain equal difference dominant frequency earthquake record amplitude spectrum；

Overlapping amplitude spectrum determining unit, for being normalized to the equal difference dominant frequency earthquake record amplitude spectrum, and to institute Stating the amplitude spectrum after normalized carries out being added acquisition overlapping amplitude spectrum；

First envelope curve value determining unit, the envelope curve for determining the overlapping amplitude spectrum, and it is bent to obtain the envelope Stacked spectrum envelope curve value on line corresponding to the equal difference dominant frequency value；

Second envelope curve value determining unit, the well logging sound wave impedance amplitude for obtaining the thin interbed model is composed, and determines institute The envelope curve of well logging sound wave impedance amplitude spectrum is stated, and obtains the well logging on the envelope curve corresponding to the equal difference dominant frequency value Sound impedance amplitude envelope curve values；

Correction unit, for utilizing the stacked spectrum envelope curve value and the well logging sound wave impedance amplitude spectrum envelope curve values pair Amplitude spectrum after the normalized is corrected treatment and obtains correction overlapping amplitude spectrum；

Earthquake Relative Wave Impedance obtaining unit, for the three rank differential according to the correction overlapping amplitude spectrum and the earthquake record Phase spectrum determines wide frequency-band spectrum, and obtains earthquake Relative Wave Impedance according to the wide frequency-band spectrum.

15. devices according to claim 14, it is characterised in that the device also includes：High frequency denoising unit, is used for High frequency denoising is carried out using the earthquake record amplitude spectrum of varying number rank differential described in default compacting function pair.

16. devices according to claim 15, it is characterised in that the default compacting function is as follows：

$T=\left\{\begin{array}{cc}{e}^{-a(f-f_{high})}& f>f_{high}\\ 1& f\≤f_{high}\end{array}\right.$

Wherein, f represents frequency variable, f_highHighest cut-off frequency is represented, a represents blanket factor.

17. according to any described device of claim 14 to 16, it is characterised in that the device also includes：Recognition unit, is used for Thin interbed is recognized according to the earthquake Relative Wave Impedance.