• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Inversion-based attenuation compensation with dip constraint

    2022-06-02 04:59:54XiongLiLiHuoGuoLiHoLiQingLongMeng
    Petroleum Science 2022年2期

    Xiong M ,Li-Li Huo ,Guo-F Li ,*,Ho Li ,Qing-Long Meng

    a School of Geophysics,China University of Petroleum-Beijing,State Key Lab of Petroleum Resources and Prospecting,Key Lab of Geophysical Exploration of CNPC,Changping,102249,Beijing,China

    b Exploration and Development Research Institute,Jidong Oilfield,China National Petroleum Corporation,Tangshan,Hebei 063000,China

    c Exploration and Development Research Institute,Dagang Oilfield,China National Petroleum Corporation,Tianjin,300280,China

    Keywords:Attenuation compensation Instability Inverse scheme Seismic dip Seismic resolution Spatial continuity

    ABSTRACT Instability is an inherent problem with the attenuation compensation methods and has been partially relieved by using the inverse scheme.However,the conventional inversion-based attenuation compensation approaches ignore the important prior information of the seismic dip.Thus,the compensated result appears to be distorted spatial continuity and has a low signal-to-noise ratio (S/N).To alleviate this issue,we have incorporated the seismic dip information into the inversion framework and have developed a dip-constrained attenuation compensation (DCAC) algorithm.The seismic dip information,calculated from the poststack seismic data,is the key to construct a dip constraint term.Benefiting from the introduction of the seismic dip constraint,the DCAC approach maintains the numerical stability and preserves the spatial continuity of the compensated result.Synthetic and field data examples demonstrate that the proposed method can not only improve seismic resolution,but also protect the continuity of seismic data.

    1.Introduction

    The attenuation phenomenon of seismic wave propagation in attenuating media has been acknowledged in theory for many years.Many attenuation models,including the Kolsky-Futterman model (Kolsky,1956;Futterman,1962),the power-law model(Strick,1967),Kjartansson's constant-Q model(Kjartansson,1979),Muller's model(Muller,1983),and the standard linear solid model(Zener,1948),have been established to describe these attenuationdispersion effects.The existence of these effects results in the distortion of a propagating wavelet and the decrease of the seismic resolution(Bickel and Natarajan 1985;Li and Wang,2007;Xie and Yang,2018).To compensate the seismic attenuation and enhance the seismogram resolution,a variety of inverse-Q filtering algorithms have been developed by many researchers(Hargreaves and Calvert,1991;Wang,2006;Oliveira and Lupinacci,2013;Chai et al.,2014;Wang and Chen,2014;Li et al.,2015;Yuan et al.,2016;Morozov et al.,2018).

    Amplitude compensation performed in inverse Q filtering is inherently unstable and can cause undesirable noise amplification as well as bandwidth enhancement in the seismic data.Therefore,the current various inverse Q filtering methods are aimed to enhance the bandwidth while suppressing the noise.For example,Robinson (1979) proposes a phase-only inverse Q filtering algorithm which only corrects the phase distortion but neglects the amplitude compensation.Since the amplitude compensation is not considered,this method is unconditionally stable.Wang (2002)deals with the instability of inverse Q filtering by adding a regularization factor to the amplitude compensation operator.Furthermore,Zhang and Ulrych (2007) use the least-squares inversion with Cauchy-Gauss type regularization to stabilize the attenuation compensation solution.Braga and Moraes (2013) accomplish the inverse Q filtering in the wavelet domain and use the L2norm constraint to stabilize the compensation result.Wang et al.(2018)exploit a nonconvex L1-2norm constraint for the stabilized absorption compensation.Nevertheless,the above inversion-based compensation methods ignore an important prior information of the structural orientation or the seismic dip,thus,the invited results may appear an incorrect or distorted spatial continuity(Li and Oldenburg,2000;Lelievre and Oldenburg,2009;Zhang et al.,2013;Gholami,2015;Yuan et al.,2015;Karimi,2015;Wang and Wang,2017;Ma et al.,2020a).

    The seismic dip,which reflects the orientations of the seismic events,is an essential attribute in a 2D seismic profile (Marfurt,2006).This information has been widely used in geophysical inversions to produce a result with improved S/N and better continuity (Yao et al.,2020).For example,Li and Oldenburg (2000)incorporate the geological dip information into resistivity inversion and induced-polarization inversion.Clapp et al.(2004)exploit the structural orientation information for the reflection tomography.Hamid and Pidlisecky (2016) use the structural angle information for the seismic impedance inversion.In addition,Zu et al.(2019) introduce the dip field information to the denoising algorithm and further propose a dip-oriented dictionary-learning approach for the random noise attenuation.

    In this paper,we incorporate the seismic dip information into the attenuation compensation algorithm,and further propose a dip-constrained attenuation compensation (DCAC) algorithm.Benefiting from the introduction of the dip constraint,the proposed DCAC method will maintain the numerical stability and improve the spatial continuity of the compensation results.The compensation performance of this approach is verified by applying both synthetic and field data examples.

    In the following of this paper,firstly,we briefly review the principle of seismic wave propagation in the attenuation medium and introduce the dip-unconstrained attenuation compensation(DUAC) algorithm.Then,we construct a dip constraint term by applying the seismic dip information.Next,we incorporate the dip constraint term into the inversion system and develop a DCAC algorithm.And then we use synthetic experiments and field data to testify the stability and superiority of the proposed DCAC method.Finally,we discuss the influence of the seismic dip estimation,the determination of the regularization parameters and Q analysis before drawing some conclusions.

    2.Methodology

    2.1.Seismic wave propagation in the attenuation medium

    In the elastic medium,the seismogram can be synthesized by the well-known convolution model which can be expressed as,

    where t and τ represent the record time,d(t)is the stationary seismic trace,w(t)is the source wavelet,and r(τ)is the reflectivity series.The schematic diagram of the convolution model is shown in Fig.1.As shown in Fig.1,the wavelet w(t)is constant when propagating in the non-attenuated medium.

    Fig.1.The schematic diagram of the convolution model:(a) the stationary seismic trace,(b) the source wavelet matrix,(c) the reflectivity series.

    When propagating in the attenuation medium,the propagation wavelet is time-varying and then the stationary convolution model(Equation (1)) can be modified as (Wang,2011):

    where s(t)is the attenuated seismic trace,and(t,τ)is the timevarying wavelet which can be expressed as (Wang,2002):

    According to the convolution theorems of Fourier transform,the frequency-domain expressions of Equation(1)and Equation(2)are,

    where D(ω)and S(ω)are the frequency-domain non-attenuated trace and the frequency-domain attenuated trace respectively,W(ω)and W(ω,τ)are the frequency-domain source wavelet and the frequency-domain time-varying wavelet respectively,and R(ω)is the frequency-domain reflectivity sequences.

    Substituting Equations (4) and (5) into Equation (6),we have,

    Transforming Equation(7)into the time domain,we obtain the following formula,

    Fig.2.The schematic diagram of the non-stationary convolution model:(a) the attenuated seismic trace,(b) the time-varying wavelet matrix,(c) the reflectivity series.

    where N is the sampling points,s and d represent the discrete attenuated record and the discrete non-attenuated record respectively,and the matrix A stands for the time-domain attenuation matrix due to the earth's Q filtering effects.An illustration of this matrix-vector product is shown in Fig.3.

    When considering a 2D seismic data with M traces,Equation(9)can be generalized to a multichannel system(Ma et al.,2020b):

    where y and m are the concatenated attenuated data vector and the concatenated non-attenuated data vector respectively,and G is the block diagonal matrix standing for the multichannel attenuation effects.

    Fig.3.An illustration of the matrix-vector product (Equation (6)):(a) the attenuated seismic trace,(b) the time-domain absorption matrix,(c) the non-attenuated seismic record.

    The main objective of the attenuation compensation is to eliminate the seismic attenuation effects from the recorded seismic data,that is,we desire to obtain the non-attenuated data m by solving Equation (10) with the known (or estimated) Q value and the recorded attenuation data y.However,the attenuation function contains an exponential decay term,which means the block attenuation matrix G is ill-conditioned and the numerical solution of Equation (10) is often unstable.

    2.2.Dip-unconstrained absorption compensation

    To obtain a numerical stable result,the Tikhonov regularization is often used to establish the following objective function(Ma et al.,2020c):

    where λ is a trade-off parameter which controls the relative strength of the constraint term (second term) to the data misfit(first term).Setting the derivative of Equation (11) with respect to m equals to zero,we obtain its solution which can be expressed as:

    where the superscript T denotes the transpose.Since this attenuation compensation algorithm lacks of the dip constraint,we refer it to as DUAC method.

    2.3.Constructing a dip constraint term

    To improve the spatial continuity of the attenuation compensation algorithm,we incorporate the additional dip information into the inversion system.The dip orientation,estimated from the post-stack seismic data,plays an important role in constructing the dip constraint term.Therefore,in this section,we focus on estimating the dip orientation and constructing a dip constraint term.

    According to Hamid and Pidlisecky (2016),the seismic dip is defined as the angle between the x-axis and the vector describing the minimum gradient direction (Fig.4).Based on the above definition,we calculate the seismic dip by the following expression:

    where θ(t,x)is the seismic dip,dxand dtare the horizontal and vertical derivatives of the seismic data respectively,and Dx(·)and Dt(·)are the horizontal and vertical derivative operators respectively.

    After estimating the seismic dip,we need to rotate the horizontal and vertical derivatives in the original coordinate system.And then we obtain two rotation operators,that is,the one parallel to the local dip and the other perpendicular to the local dip.According to Fig.4,the rotated and original derivative operators satisfy the following coordinate transformation formula:

    where Dparlis the derivative operator parallel to the local dip direction and Dperpis the derivative operator perpendicular to the local dip direction.The matrix-vector form of Equation(15)can be expressed as,

    The derivative parallel to the local dip can be used to enforce the smoothness and continuity of the seismic data in the dip direction,thus,we can construct the following dip constraint term:

    Fig.4.An illustration for the seismic dip estimation.The blue curve represents a seismic event.dx and dt are the horizontal and vertical derivatives of the seismic event respectively.θ is the seismic dip.dparl is the derivative which is oriented parallel to the local dip direction and dperp is the derivative perpendicular to the local dip direction.

    To better understand this constraint,we first consider a special case of the horizontally layered medium.In this case,the seismic dip is 0 and the operator Dparlwill degenerate to the horizontal derivative matrix Dx.As discussed by Hamid and Pidlisecky(2015),the effect of the horizontal derivative matrix Dxis to minimize the differences between the seismic data in the horizontal direction.This means the lateral constraintcan enforce the lateral continuity of the inverted result m (because the inverted result with the horizontal events can minimize the lateral constraint).Compared with the horizontal derivative matrix Dx,the dip-related derivative matrix Dparlcontains the extra dip information at each sampling point.Therefore,we can find that the effect of the diprelated derivative matrix Dparlis to minimize the differences between the seismic data in the local dip direction.Similarly,the dip constraintcan enforce the spatial continuity of the inverted result m along the dip direction (because the inverted result m with the similar seismic dip field to the dip-related derivative matrix Dparlcan minimize the dip constraint).It can be seen that the core issue of this dip constraint is the reliability of the seismic dip information.Theoretically,the most reliable seismic dip should be calculated from the true inverted result (or the nonattenuated data).In the field data,we do not have the true inverted result (or the non-attenuated data),but only have the attenuated data.Generally,the difference of the dip information between the attenuated and the compensated profiles will not be great.Thus,we can estimate the seismic dip information from the attenuated data.In summary,the role of the dip constraint is to ensure that the inverted results are structurally similar to the attenuated seismic data.

    2.4.Dip-constrained attenuation compensation

    The DUAC algorithm ignores the seismic dip information,thus,the inverted results may appear an incorrect or distorted spatial continuity when seismic data are contaminated by the random noise.To eliminate this issue,we incorporate the dip constraint terms (Equation (17)) into the inversion system and establish the following objective function:

    where μ is the regularization parameter controlling the relative strength of the dip constraint (third term) to the data misfit term(first term).The least-squares solution of this problem is:

    Since the proposed attenuation compensation algorithm includes a dip constraint term,we refer it to as DCAC method.It needs to be explained that although we use the matrix inversion to express the numerical results of the DCAC method,we do not directly use the inversion operation in programming(because the inversion operation is computationally expensive).Instead,we use the conjugate gradient method to solve the equivalent equation=GTy.Therefore,the algorithm has high computational efficiency.

    The objective function of the DCAC method contains three terms,that is,the data misfit termthe vertical smoothness constraintand the dip constraintThe data misfit term contorls the inversion accuracy.The second term requires the inversion results to be smoothness (for suppressing the amplification of the high frequencies) and the parameter λ contorls the smoothing strength of the compensation results.The parameter λ has a similar action to the stabilization factor in Wang's method(Wang,2002).If the parameter λ is small,more high frequencies (both the effective signal and noise) can be recovered or amplified.In contrast,if the parameter λ is large,the high-frequency recovery or amplification(both the effective signals and noises) will be suppressed.For any frequency component,the magnification of the effective signals and noises is the same.In the proposed method,we have an additional dip constraint.The role of the dip constraint is to ensure that the compensation results are structurally similar to the seismic section before compensation.In other word,the dip constraint term is mainly to enhance the signals consistent with the seismic dip (or suppress the signals inconsistent with the seismic dip).This means the dip constraint can selectively enhance the effective signals while suppress the seismic noises(at the same bandwidth).In summary,the data misfit term is to ensure the accuracy of the inversion results.The second term is mainly to improve the seismic resolution and partially suppress the amplification of the high frequencies (similar to a low-pass filtering).The third term is to further selectively suppress the signals inconsistent with the seismic dip (similar to a dip-filtering).

    3.Example

    3.1.Synthetic data experiments

    To testify the validity and superiority of the DCAC algorithm,we use Marmousi model shown in Fig.5a to conduct the attenuation compensation tests.Assuming that the model has a constant density,we obtain the reflectivity model displayed in Fig.5b.Convolving a 30 Hz Ricker wavelet with the above reflectivity model provides the non-attenuated seismic data(Fig.5c).The nonattenuated data can be served as the reference data to evaluate the compensation performance of Wang's method (Wang,2002),the DUAC method and the DCAC approach.Fig.5d shows the noisy attenuated seismogram with the quality factor Q=40 and contaminated by 20%random noise.This noisy data is exploited to conduct the compensation experiments and used to analyze the performance of Wang's,the DUAC and the DCAC approaches.We find from Fig.5d that the resolution of the attenuated seismic data is decreased due to the seismic attenuation,and the reflection events in the deep are almost drowned in the random noise.To compensate the seismic attenuation effects and enhance seismic resolution,we apply Wang's method,the DUAC approach and the DCAC algorithm to process the attenuated seismic data.

    In order to implement the proposed DCAC algorithm,we first calculate the seismic dip by the following steps:1)Computing the vertical and horizontal derivatives of the attenuated seismic data.Note that the derivative calculation will amplify the seismic noise,thus,we estimate the derivatives from the smoothed seismic data.The obtained vertical and horizontal derivatives are shown in Fig.6a and b respectively.2) Estimating the seismic dip at every location using Equation(13)and depicting the result in Fig.6c.3)As we can see,the estimated dips are also effected by the seismic noise,thus,we use spatial Gaussian filter to smooth them and show the smoothed seismic dip section in Fig.6d.Based on the seismic dip information,we then construct the dip constraint and further implement the DCAC algorithm to compensate the seismic attenuation.For quantitatively assessing the compensation results,we use the average correlation coefficient (ACC) as an indicator:

    where djandare the reference (non-attenuated) and the compensated traces respectively.

    In Wang's method,we choose σ2=0.01 and the inverse-Q filtering result is showed in Fig.7a.We can find that the seismic amplitude is enhanced especially in the mid-and deep-layers,and the spatial continuity of the compensated profile is slightly low.The ACC between it and the non-attenuated seismic data (Fig.5c) is 0.5975.In the DUAC method,we select the regularization parameter λ=0.007 and display the corresponding compensation result in Fig.7b.We see that the compensated result improves the seismic energy and seismic resolution,but the amplification of the high frequencies (both the effective signals and noises) is also evident.The ACC between it and the non-attenuated seismic data(Fig.5c)is 0.6124.To prove the superiority of the DCAC algorithm,we apply the proposed DCAC algorithm to process the attenuated data.In the DCAC algorithm,we fix the parameter λ=0.007 and determine the dip constraint parameter μ=0.1 by trial and error.The DCAC compensation result is shown in Fig.7c.Compared with Wang's and the DUAC compensation results,the DCAC compensation result shows better spatial continuity (see red arrows) and maintains a relatively higher S/N.Its ACC with the ideal data reaches 0.8602.

    To clearly view the differences between the reference and the compensated data,we display the seismic traces extracted from the non-attenuated data(Fig.5c)and the compensated data(Fig.7) at CDP=101 and their spectra in Fig.8.Comparing the compensated results with the reference trace,we find that Wang's and the DUAC algorithms overcompensate the attenuated seismic trace due to the amplification of the high frequencies(both the effective signals and noises).The comparison of their spectra further confirms that the proposed method enhances the seismic bandwidth and maintains a relatively higher S/N.

    We also use the attenuated data shown in Fig.5d to investigate the influence of the dip constraint parameter on the compensation results.In this experiment,we fix the parameters λ=0.005 and respectively choose μ as 1,0.1,0.01,and 0.001.The corresponding compensated data are displayed in Fig.9a-d.As shown in Fig.9a,when the dip constraint parameter μ is too large,the compensated profile appears to be over-smoothed and the resolution is slightly low.The ACC between it and the non-attenuated seismic data(Fig.5c)is 0.7874.However,when the parameter μ is too small,the amplification of the high frequencies is evident and the spatial continuity of the compensated result is poor (Fig.9d).In this case,the ACC between it and the non-attenuated seismic data (Fig.5c)becomes 0.8103.When the parameter μ is moderate(Fig.9b and c),the compensated results achieve a good balance between the noise suppression and the resolution enhancement.And the ACC between them and the non-attenuated seismic data (Fig.5c) are 0.8477 and 0.8298 respectively.This means the selection of the good regularization parameters is important to the compensation results and we will discuss the parameter selection strategy in the DISCUSSION.

    Fig.5.Forward modeling for generating the synthetic data.(a)The velocity model,(b)the reflectivity model,(c)the synthetic noise-free non-attenuated seismic data serving as the reference,(d) the attenuated data with 20% Gaussian noise.

    Fig.6.Estimation of the seismic dip.(a)The second derivative of the attenuated seismic data in the vertical direction,(b)the second derivative of the attenuated seismic data in the horizontal direction,(c) the estimated dip,(d) the smoothed seismic dip.

    Fig.7.Attenuation compensation for the noisy attenuated seismic data.(a)Wang's method with σ2=0.01,(b)the DUAC algorithm with the regularization parameter λ=0.007,(c)the DCAC algorithm with the regularization parameters λ=0.007 and μ=0.1.

    The accuracy and stability of the attenuation compensation results are affected by the noise level.Therefore,we compare the compensation results of Wang's method,the DUAC algorithm and the DCAC approach at different noise levels (5%,10%,15%,20% and 25% random noises).Different methods are used to process the noisy records and the parameters in these methods are fixed at different cases for an objective comparison.In Wang's method,we choose the stabilization factor σ2=0.01.In the DUAC method,the regularization parameter is λ=0.007.In the DCAC method,we fix the parameter λ=0.007 and select the dip constraint parameter μ=0.1.The ACCs of the compensation results are reported in Table 1.We can see that when the noise level is low(5%and 10%),all three methods provide a good result with high ACCs.With the noise level increases,the ACCs of the attenuation compensation results are decrease.Furthermore,for the relatively strong nosie cases(20%and 25%),the ACCs are decrease faster than that in the low noise levels.The main reason is that the high frequencies are partially over-compensated with above parameters for strong nosie cases.And in strong nosie cases,the accuracy of the dip estimation is decreased,which will further influence the accuracy of the proposed method.For the influence of the seismic dip estimation,we will discuss in the DISCUSSION.In addition,we find that the proposed method outperformes other methods in all cases,which further proves the superiority of the proposed method.

    Table 1The ACC of the attenuation compensation results using different methods at different noise levels.

    3.2.Field data tests

    The field data shown in Fig.10a are used to further demonstrate the practicability of the DCAC algorithm.To compensate the energy attenuation and enhance seismic resolution,we exploit Wang's,the DUAC and the DCAC methods to process the original seismic profile.In Wang's method,we select σ2=0.02 and the inverse-Q filtering result is showed in Fig.10b.For the DUAC algorithm,we choose the regularization parameter λ=0.008.While in the DCAC algorithm,we fix the parameter λ=0.008,and then determine μ=0.2 by trial and error.The compensation results using the DUAC and DCAC methods are shown in Figs.10c and d respectively.The averaged amplitude spectra of the three compensation results are displayed in Fig.10e.From the five Figs.,we can draw the following conclusions:Firstly,we find that all three methods compensate the seismic attenuation and enhance the resolution of the raw data.Secondly,compared with Wang's and the DUAC compensation results,the DCAC result tend to have a higher S/R and smoother spatial continuity without losing evident vertical resolution,which indicates that the proposed DCAC algorithm enjoys better compensation performance.Thirdly,we see from the spectra that the frequency components of Wang's and the DUAC results are slightly higher than that of the DCAC result,which may cause by the suppression of the signals inconsistent with the seismic dip in the DCAC algorithm.

    We also extract three reference traces from the original and the compensated data at CDP=50,CDP=150,and CDP=250 to further evaluate the compensation details of three compensation results.Note that we only show the seismic records from 2 s to 5 s and the original traces (green) are gained by a factor of 2 for a clearer comparison.As shown in Fig.11,all three methods improve the seismic energy in the deep reflections and they show similar performance in enhancing the vertical resolution.However,the DCAC result shows a relatively small amplitude compensation(see arrows) which may result from the suppression of the additional seismic noises inconsistent with the seismic dip.It indicates that the proposed DCAC method can not only recover the seismic events and improve the seismic resolution,but also keep the seismic noise at an acceptable level.

    4.Discussions

    Fig.8.The reference and the compensated traces and their spectra.(a)The reference and the compensated traces extracted from the reference data(Fig.5c)and the compensated data (Fig.7) at CDP=101,and (b) their spectra.

    Fig.9.The investigation of the influence of the dip constraint parameter μ on the compensation results.We fix the parameters λ=0.005 and respectively choose μ as(a)1,(b)0.1,(c) 0.01,and (d) 0.001.

    In the proposed DCAC method,there are two major factors affecting the compensation results,that is,the accuracy of the seismic dip estimation and the determination of the regularization parameters.Firstly,we discuss the influence of the seismic dip estimation.The structural dip information can be derived from many sources (e.g.geologic data,borehole data).In this paper,we derive the dip information directly from the seismic data.In the past few decades,many researchers have proposed or developed a large number of methods to estimate the seismic dip from the seismic data.For example,Bahorich and Farmer(1995)estimate the seismic dip based on the crosscorrelation of a set of windowed seismic data generated by using time lagging between nearby seismic traces.This algorithm is simple and efficient,but the stability of the seismic dip estimation is relatively poor.To improve the stability and accuracy of the dip estimation,many novel methods,including the semblance-based multiple window scanning methods(Marfurt et al.,1998),the structure-tensor-based methods(Fehmers and H¨ocker,2003;Wu,2017),and the plane-wave destruction (Fomel,2002),have been developed in recent years.Considering that the focus of this paper is not on developing a new method for the seismic dip estimation,therefore we select a relatively simple but efficient dip estimation algorithm proposed by Hamid and Pidlisecky (2016).For mid-and large-scale geologic structures,this method can provide a relatively reliable and stable result even if the seismic data are contaminated by some seismic noises.If the S/N of the seismic record is low,we have to increase the seismic noise attenuation before the seismic dip estimation to ensure the stability.In this case,the accuracy of the dip estimation will decrease.Actually,the seismic dip estimation is still an open and undergoing subject in geophysics and we will further research this problem in the future.In addition,we can incorporate any other dip estimation method into the proposed compensation algorithm by replacing the dip estimation method in this paper with the other selected method.This also shows that the proposed DCAC method has good flexibility.

    Fig.10.Seismic attenuation compensation for the field data.(a) The field data,(b) the compensation result from Wang's method,(c) the compensation result from the DUAC method,(d) the compensation result from the DCAC approach,(e) the amplitude spectra of the original data,Wang's result,the DUAC result and the DCAC result.

    Fig.11.Comparison of the compensation performance using three reference traces extracted from Figs.10a-d at(a)CDP=50,(b)CDP=150 and(c)CDP=250.The original traces(green) are displayed with a gain factor of 2.

    Then we discuss the strategy of selecting the good regularization parameters.In the DUAC algorithm,there is a regularization parameter λ to be determined.As discussed by Wang(2011),when the S/N of the seismic data is high,we can choose a small parameter λ to maximally enhance the bandwidth of the seismic data.While the S/N of the seismic data is lower,we should select a relatively larger value for suppressing the amplitude amplification of the high frequencies.This is the general principle for selecting a good regularization parameter λ.To quantitatively determine a suitable regularization parameter λ,Hansen and O'Leary(1993)propose a Lcurve technique,which applies a cross-plot of the data error versus the solution length as a function of λ.A good value for the parameter λ is the one located at the corner of the L-curve.For the proposed DCAC method,there are two regularization parameters,λ and μ,to be determined.The role of the parameter λ is the same as that in the DUAC algorithm,while the role of the parameter μ is to enhance the signals consistent with the seismic dip (or suppress the signals inconsistent with the seismic dip).Until now,the optimization for the hyperparameters functional remains a complicated problem(Du et al.,2018).In this paper,we select them by trial and error,but we apply a relatively elegant strategy.Firstly,we set the dip regularization parameter μ=0 (that means ignoring the dip constraint term)and use L-curve technique to choose a proper parameter λ.Secondly,we fix the parameter λ calculated in the first step and then search a rough value of the parameter μ from a wide range using a large search-step.Finally,we minorly adjust the value of λ and μ by trial and error to determine the most suitable parameters.

    In addition,the Q estimation (or Q analysis) is a basic work for the attenuation compensation.Wang (2004) has proposed two methods(the attenuation-based Q analysis and the compensationbased Q analysis)for Q estimation from the reflection seismic data.These methods are the single-trace Q eatimation methods.In the proposed DCAC algorithm,we require a 2D Q model as the input.But the Q analysis processing is not necessarily based on the multichannel model.We can use the attenuation-based Q analysis or the compensation-based Q analysis to estimate Q value trace-bytrace,and then form the single-trace results to a 2D Q model.As showed in Equation(10),the multichannel model can be expressed asActually,this expression can be rewritten asFor the first trace,the Q estimation result will be included in the matrix A1.Similarly,for the second trace,the Q estimation result will be included in the matrix A2.And for the M-th trace,the Q estimation result will be included in the matrix AM.Therefore,the Q analysis is not necessarily on the multichannel base and the compensation results would not be affected while using the single-trace Q eatimation method.

    5.Conclusions

    In this paper,we propose a novel DCAC method by incorporating the dip constraint into the inversion-based compensation algorithm.The dip regularization term,calculated from the rotated horizontal and vertical derivatives,plays an important role in improving the spatial continuity of the compensated result.Compared with Wang's and DUAC approaches,the proposed DCAC method has advantages in protecting the spatial continuity and maintaining a relatively higher S/R of the compensation data.The superior performance of the proposed approach is demonstrated by using both synthetic tests and field data.The synthetic tests indicate that the proposed method can not only enhance the vertical resolution but also improve the spatial continuity.The field data experiment further confirms the practicability and stability of the DCAC algorithm.

    Acknowledgments

    We thank the financial support provided by National Natural Science Foundation of China(42074141),the Strategic Cooperation Technology Projects of CNPC and CUPB (ZLZX2020-03),and National Key R &D Program of China (2018YFA0702504).We also gratefully acknowledge the helpful comments from the editors and anonymous reviewers,which greatly improved this manuscript.

    韩国高清视频一区二区三区| 欧美日韩综合久久久久久| 狂野欧美激情性bbbbbb| 亚洲欧洲日产国产| 最黄视频免费看| 99久久精品一区二区三区| 观看美女的网站| av免费观看日本| 成人毛片a级毛片在线播放| 久久精品国产亚洲av涩爱| 国产亚洲精品第一综合不卡 | 精品久久国产蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 边亲边吃奶的免费视频| 国产av一区二区精品久久| 秋霞在线观看毛片| 免费观看无遮挡的男女| av免费在线看不卡| 一级,二级,三级黄色视频| 国产一区二区三区综合在线观看 | 欧美 亚洲 国产 日韩一| 日韩强制内射视频| 精品一区二区免费观看| 亚洲精品视频女| 日韩精品免费视频一区二区三区 | 免费久久久久久久精品成人欧美视频 | 亚州av有码| 免费高清在线观看视频在线观看| 看十八女毛片水多多多| 日本欧美视频一区| 九九爱精品视频在线观看| 久久久久国产精品人妻一区二区| 91久久精品国产一区二区成人| av一本久久久久| 国产极品粉嫩免费观看在线 | 91aial.com中文字幕在线观看| 欧美三级亚洲精品| 久久精品国产亚洲av涩爱| 久久亚洲国产成人精品v| 国产在线一区二区三区精| 精品一区二区三区视频在线| 欧美国产精品一级二级三级| 日韩中文字幕视频在线看片| 日本av免费视频播放| 18+在线观看网站| 韩国av在线不卡| 日日爽夜夜爽网站| 亚洲一级一片aⅴ在线观看| 国产一区二区在线观看av| 日韩中字成人| 最新中文字幕久久久久| 亚洲精品日韩在线中文字幕| 国产在线免费精品| 亚洲av二区三区四区| 伦精品一区二区三区| 日本91视频免费播放| 国产日韩欧美亚洲二区| 特大巨黑吊av在线直播| 久久久久国产精品人妻一区二区| 99九九在线精品视频| 日韩av免费高清视频| 最新中文字幕久久久久| 中文欧美无线码| av网站免费在线观看视频| 十八禁网站网址无遮挡| 精品视频人人做人人爽| 欧美 日韩 精品 国产| 美女cb高潮喷水在线观看| 综合色丁香网| 亚洲欧美日韩卡通动漫| 一二三四中文在线观看免费高清| 国产精品欧美亚洲77777| 日韩大片免费观看网站| 午夜激情久久久久久久| 国产精品国产三级国产av玫瑰| 3wmmmm亚洲av在线观看| 亚洲精品自拍成人| 成人黄色视频免费在线看| 五月伊人婷婷丁香| 国产一区二区在线观看av| 日韩三级伦理在线观看| 免费播放大片免费观看视频在线观看| 汤姆久久久久久久影院中文字幕| 中文天堂在线官网| 亚洲av国产av综合av卡| 国产精品 国内视频| 夜夜看夜夜爽夜夜摸| 国产有黄有色有爽视频| 国产成人精品福利久久| 伊人久久精品亚洲午夜| 91精品一卡2卡3卡4卡| 午夜免费男女啪啪视频观看| 国产精品久久久久久精品古装| 女性被躁到高潮视频| 美女内射精品一级片tv| 综合色丁香网| 黄色毛片三级朝国网站| 91精品国产国语对白视频| 91精品国产国语对白视频| 高清不卡的av网站| 91久久精品国产一区二区成人| 18禁在线无遮挡免费观看视频| 亚洲国产精品成人久久小说| 一区二区三区免费毛片| 人人妻人人添人人爽欧美一区卜| 亚洲情色 制服丝袜| 国产国拍精品亚洲av在线观看| 成人国产av品久久久| 欧美少妇被猛烈插入视频| 菩萨蛮人人尽说江南好唐韦庄| 搡女人真爽免费视频火全软件| 久久精品久久久久久噜噜老黄| 日日撸夜夜添| 少妇的逼好多水| 男女边吃奶边做爰视频| 有码 亚洲区| 青春草亚洲视频在线观看| 国产高清有码在线观看视频| 中文字幕最新亚洲高清| 久久久久视频综合| 菩萨蛮人人尽说江南好唐韦庄| 这个男人来自地球电影免费观看 | 午夜免费观看性视频| 考比视频在线观看| 美女中出高潮动态图| 天天影视国产精品| 亚洲精品国产av蜜桃| 天堂俺去俺来也www色官网| 妹子高潮喷水视频| 妹子高潮喷水视频| 一级黄片播放器| 欧美精品亚洲一区二区| 丝袜脚勾引网站| 男人爽女人下面视频在线观看| 能在线免费看毛片的网站| 久久午夜福利片| 免费观看无遮挡的男女| 精品国产国语对白av| 99久国产av精品国产电影| 热99国产精品久久久久久7| 人妻少妇偷人精品九色| 免费日韩欧美在线观看| 少妇人妻精品综合一区二区| 日韩一区二区三区影片| 午夜精品国产一区二区电影| 91精品一卡2卡3卡4卡| 全区人妻精品视频| 制服丝袜香蕉在线| 男的添女的下面高潮视频| 亚洲精品自拍成人| 如日韩欧美国产精品一区二区三区 | 麻豆成人av视频| 国产精品.久久久| 国产永久视频网站| 美女大奶头黄色视频| 婷婷色综合大香蕉| 草草在线视频免费看| 观看av在线不卡| 一级毛片 在线播放| 最新中文字幕久久久久| av福利片在线| 国产一区二区三区av在线| 国产高清不卡午夜福利| 久久人妻熟女aⅴ| 亚洲成色77777| 五月开心婷婷网| 人人妻人人爽人人添夜夜欢视频| 精品视频人人做人人爽| 久久精品国产亚洲av涩爱| 久久久久国产网址| 成年美女黄网站色视频大全免费 | 亚洲精品av麻豆狂野| 亚洲av中文av极速乱| 最黄视频免费看| 一级二级三级毛片免费看| 在线 av 中文字幕| 久久99热6这里只有精品| 久久久久久久国产电影| 国产又色又爽无遮挡免| 国产白丝娇喘喷水9色精品| 波野结衣二区三区在线| 在线看a的网站| 亚洲精品日韩在线中文字幕| 久久精品国产自在天天线| 亚洲伊人久久精品综合| 精品午夜福利在线看| 九色成人免费人妻av| 99久久综合免费| 久久精品国产亚洲网站| 久久精品国产a三级三级三级| 夜夜爽夜夜爽视频| 欧美日本中文国产一区发布| 国产精品久久久久久久电影| 大话2 男鬼变身卡| 亚洲精品久久久久久婷婷小说| av在线观看视频网站免费| 婷婷色综合大香蕉| 中文字幕最新亚洲高清| 欧美少妇被猛烈插入视频| 中文字幕最新亚洲高清| 观看美女的网站| 成人18禁高潮啪啪吃奶动态图 | 国内精品宾馆在线| av电影中文网址| 午夜福利视频精品| av电影中文网址| 国产精品99久久99久久久不卡 | 亚洲欧美日韩另类电影网站| 欧美丝袜亚洲另类| 中文字幕人妻熟人妻熟丝袜美| 午夜免费观看性视频| 97超碰精品成人国产| 国产在线免费精品| 边亲边吃奶的免费视频| 国产 一区精品| 视频区图区小说| 人人妻人人澡人人看| 777米奇影视久久| 制服人妻中文乱码| 亚洲成人手机| 午夜激情av网站| 亚洲精品中文字幕在线视频| 黄色欧美视频在线观看| 性高湖久久久久久久久免费观看| 91久久精品电影网| 欧美人与性动交α欧美精品济南到 | 亚洲av中文av极速乱| 日韩中字成人| 中国国产av一级| 中文乱码字字幕精品一区二区三区| 美女视频免费永久观看网站| 一级爰片在线观看| 人妻制服诱惑在线中文字幕| 满18在线观看网站| 国产精品一国产av| 女人精品久久久久毛片| 99九九线精品视频在线观看视频| av在线播放精品| 最后的刺客免费高清国语| 国内精品宾馆在线| 男人爽女人下面视频在线观看| 亚洲精品中文字幕在线视频| 欧美精品亚洲一区二区| 日本黄色片子视频| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线| 日本色播在线视频| 天堂8中文在线网| 插阴视频在线观看视频| 一本一本综合久久| 特大巨黑吊av在线直播| 青春草视频在线免费观看| 三级国产精品欧美在线观看| 九九久久精品国产亚洲av麻豆| 国产一区亚洲一区在线观看| 欧美性感艳星| 一级毛片我不卡| 国产精品国产三级国产专区5o| 国产日韩欧美亚洲二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 蜜桃在线观看..| 九九在线视频观看精品| 汤姆久久久久久久影院中文字幕| 爱豆传媒免费全集在线观看| 天堂8中文在线网| a级毛片黄视频| 黑人巨大精品欧美一区二区蜜桃 | 亚洲不卡免费看| 不卡视频在线观看欧美| 成人毛片a级毛片在线播放| 亚洲av成人精品一二三区| 一级二级三级毛片免费看| 精品人妻熟女av久视频| 国产成人精品一,二区| 少妇高潮的动态图| 黑丝袜美女国产一区| 欧美精品人与动牲交sv欧美| 激情五月婷婷亚洲| 国产老妇伦熟女老妇高清| 老女人水多毛片| 亚洲人成77777在线视频| 五月伊人婷婷丁香| 亚洲情色 制服丝袜| 国产精品99久久99久久久不卡 | 国模一区二区三区四区视频| 欧美xxxx性猛交bbbb| av黄色大香蕉| 久久久久国产网址| 天天操日日干夜夜撸| 久久久精品94久久精品| 午夜激情av网站| 中文字幕制服av| 人妻 亚洲 视频| 久久99蜜桃精品久久| 天美传媒精品一区二区| 夜夜骑夜夜射夜夜干| 一二三四在线观看免费中文在| videos熟女内射| 久久人妻熟女aⅴ| 成人黄色视频免费在线看| aaaaa片日本免费| av国产精品久久久久影院| 欧美激情久久久久久爽电影 | 美国免费a级毛片| 国产日韩欧美亚洲二区| 欧美激情极品国产一区二区三区| 亚洲熟女精品中文字幕| 久久国产精品人妻蜜桃| 一边摸一边做爽爽视频免费| 亚洲熟妇熟女久久| 久久久国产成人免费| 黄色怎么调成土黄色| 午夜两性在线视频| 精品乱码久久久久久99久播| 国产黄频视频在线观看| 欧美精品av麻豆av| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品男人的天堂亚洲| 性色av乱码一区二区三区2| 99热国产这里只有精品6| 制服人妻中文乱码| 国产成+人综合+亚洲专区| 国产精品久久久久久精品电影小说| 日韩一区二区三区影片| 精品国产一区二区三区四区第35| 中文字幕另类日韩欧美亚洲嫩草| 夜夜夜夜夜久久久久| 色播在线永久视频| 日韩制服丝袜自拍偷拍| 国产精品一区二区在线观看99| 国产欧美日韩综合在线一区二区| 亚洲专区字幕在线| 乱人伦中国视频| 五月开心婷婷网| 国产有黄有色有爽视频| 青草久久国产| 菩萨蛮人人尽说江南好唐韦庄| 国产精品欧美亚洲77777| 下体分泌物呈黄色| 日本精品一区二区三区蜜桃| 久久天躁狠狠躁夜夜2o2o| 久久久精品94久久精品| 国产精品久久久久久精品古装| kizo精华| 男人操女人黄网站| 我的亚洲天堂| 日韩制服丝袜自拍偷拍| 国产三级黄色录像| 操出白浆在线播放| 咕卡用的链子| 丝瓜视频免费看黄片| 在线天堂中文资源库| 中国美女看黄片| 亚洲人成伊人成综合网2020| 国产成人免费无遮挡视频| 91字幕亚洲| kizo精华| 久久亚洲真实| 亚洲第一欧美日韩一区二区三区 | netflix在线观看网站| 精品国产亚洲在线| 国产欧美日韩精品亚洲av| 国产精品成人在线| svipshipincom国产片| 下体分泌物呈黄色| 精品国产一区二区三区四区第35| 亚洲国产精品一区二区三区在线| 法律面前人人平等表现在哪些方面| 侵犯人妻中文字幕一二三四区| 亚洲精品粉嫩美女一区| 99热网站在线观看| 黑人操中国人逼视频| 精品亚洲成a人片在线观看| 日韩熟女老妇一区二区性免费视频| 中文字幕人妻熟女乱码| 亚洲色图综合在线观看| 亚洲 国产 在线| 亚洲黑人精品在线| 9色porny在线观看| 女人被躁到高潮嗷嗷叫费观| av有码第一页| 亚洲免费av在线视频| 9热在线视频观看99| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 一本久久精品| 成年动漫av网址| 久久久久国内视频| 一区二区三区国产精品乱码| 极品少妇高潮喷水抽搐| 91成年电影在线观看| 亚洲免费av在线视频| 国产在视频线精品| 在线观看人妻少妇| 国内毛片毛片毛片毛片毛片| 一区二区av电影网| 香蕉久久夜色| 国产在线一区二区三区精| 久热这里只有精品99| 1024视频免费在线观看| 国产aⅴ精品一区二区三区波| 久久人人97超碰香蕉20202| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 999久久久精品免费观看国产| 肉色欧美久久久久久久蜜桃| 2018国产大陆天天弄谢| 妹子高潮喷水视频| 黑丝袜美女国产一区| 国产精品 欧美亚洲| 亚洲精品中文字幕在线视频| 久久久水蜜桃国产精品网| 国产极品粉嫩免费观看在线| 午夜日韩欧美国产| 免费少妇av软件| 亚洲av欧美aⅴ国产| 美女扒开内裤让男人捅视频| 一边摸一边抽搐一进一出视频| 精品一品国产午夜福利视频| 久久精品亚洲av国产电影网| 亚洲国产欧美网| 国产av一区二区精品久久| 免费观看a级毛片全部| 成人永久免费在线观看视频 | 日韩大码丰满熟妇| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 别揉我奶头~嗯~啊~动态视频| 国产成+人综合+亚洲专区| 日韩欧美一区二区三区在线观看 | kizo精华| 精品国产一区二区久久| 国产精品熟女久久久久浪| 国产精品久久久人人做人人爽| 国产黄频视频在线观看| 色播在线永久视频| 一区福利在线观看| 欧美日韩视频精品一区| 久久精品亚洲精品国产色婷小说| av在线播放免费不卡| 国产黄频视频在线观看| 99九九在线精品视频| 精品人妻在线不人妻| 后天国语完整版免费观看| 亚洲av成人一区二区三| 久久亚洲精品不卡| 后天国语完整版免费观看| 1024香蕉在线观看| 亚洲天堂av无毛| 丰满少妇做爰视频| 两个人免费观看高清视频| 成人手机av| 热re99久久国产66热| 一个人免费看片子| 最新的欧美精品一区二区| 亚洲 欧美一区二区三区| av天堂在线播放| 亚洲av欧美aⅴ国产| 黑人欧美特级aaaaaa片| 纯流量卡能插随身wifi吗| 国产亚洲精品久久久久5区| 国产一区二区 视频在线| 久久久国产欧美日韩av| 中文欧美无线码| 亚洲自偷自拍图片 自拍| 亚洲精品一二三| 亚洲成a人片在线一区二区| 亚洲欧美日韩另类电影网站| 五月开心婷婷网| 多毛熟女@视频| 成人av一区二区三区在线看| 免费女性裸体啪啪无遮挡网站| 亚洲精品久久午夜乱码| 日韩视频在线欧美| 中文字幕人妻熟女乱码| 成年人午夜在线观看视频| 欧美激情高清一区二区三区| 一级毛片精品| 欧美黄色片欧美黄色片| 国产激情久久老熟女| 色在线成人网| 黄色片一级片一级黄色片| 热re99久久精品国产66热6| 国产精品av久久久久免费| 9色porny在线观看| 国产真人三级小视频在线观看| 女警被强在线播放| 国产精品一区二区精品视频观看| 99精国产麻豆久久婷婷| 一级毛片精品| 日韩三级视频一区二区三区| 久久久精品94久久精品| 亚洲美女黄片视频| 这个男人来自地球电影免费观看| 成人永久免费在线观看视频 | 一本色道久久久久久精品综合| 久久婷婷成人综合色麻豆| 亚洲精品自拍成人| 在线十欧美十亚洲十日本专区| 男女床上黄色一级片免费看| 精品欧美一区二区三区在线| 久久人人97超碰香蕉20202| 电影成人av| 美国免费a级毛片| 90打野战视频偷拍视频| 中文字幕最新亚洲高清| 波多野结衣av一区二区av| 久久这里只有精品19| 妹子高潮喷水视频| 电影成人av| 亚洲va日本ⅴa欧美va伊人久久| 老汉色av国产亚洲站长工具| 国产高清videossex| 亚洲欧美色中文字幕在线| 视频区图区小说| 天堂俺去俺来也www色官网| 亚洲av国产av综合av卡| 欧美激情久久久久久爽电影 | 51午夜福利影视在线观看| 一区二区三区国产精品乱码| 国产伦理片在线播放av一区| 精品久久蜜臀av无| www.自偷自拍.com| 波多野结衣一区麻豆| 亚洲国产欧美日韩在线播放| 免费看十八禁软件| 国产精品99久久99久久久不卡| 51午夜福利影视在线观看| 亚洲av成人一区二区三| 亚洲中文av在线| 国产99久久九九免费精品| 18禁国产床啪视频网站| 免费久久久久久久精品成人欧美视频| 亚洲成人免费av在线播放| 国产麻豆69| 男女无遮挡免费网站观看| 人人妻人人澡人人看| 另类精品久久| 成年人黄色毛片网站| 侵犯人妻中文字幕一二三四区| 人成视频在线观看免费观看| 国产老妇伦熟女老妇高清| 亚洲精品中文字幕一二三四区 | 亚洲,欧美精品.| av天堂久久9| 午夜福利乱码中文字幕| av网站在线播放免费| 久久狼人影院| 久久久久久亚洲精品国产蜜桃av| 精品人妻1区二区| 99久久国产精品久久久| 成人国产一区最新在线观看| 后天国语完整版免费观看| 激情视频va一区二区三区| 亚洲av第一区精品v没综合| 美国免费a级毛片| 大陆偷拍与自拍| 精品少妇内射三级| 捣出白浆h1v1| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频| 国产高清videossex| 考比视频在线观看| 亚洲欧洲日产国产| 欧美变态另类bdsm刘玥| 久久午夜综合久久蜜桃| 国产精品国产av在线观看| 成年女人毛片免费观看观看9 | 91精品国产国语对白视频| 久久久水蜜桃国产精品网| 91老司机精品| 91麻豆精品激情在线观看国产 | 精品人妻在线不人妻| 搡老乐熟女国产| 在线观看免费视频网站a站| 欧美成人免费av一区二区三区 | 亚洲午夜精品一区,二区,三区| 亚洲七黄色美女视频| 免费在线观看视频国产中文字幕亚洲| 男女之事视频高清在线观看| av不卡在线播放| 男女免费视频国产| 国产精品国产高清国产av | 男女下面插进去视频免费观看| 国产精品久久久久久精品电影小说| 欧美在线一区亚洲| 黄色成人免费大全| 国产黄色免费在线视频| 亚洲三区欧美一区| 免费看十八禁软件| 麻豆乱淫一区二区| 国产色视频综合| 亚洲精品国产区一区二| 日韩有码中文字幕| 欧美+亚洲+日韩+国产| 99热网站在线观看| 丝袜喷水一区| 亚洲国产欧美在线一区| 国产淫语在线视频| 伦理电影免费视频| 亚洲久久久国产精品| 国产精品电影一区二区三区 | 日本精品一区二区三区蜜桃| 一级毛片女人18水好多| 99国产综合亚洲精品| a级毛片黄视频| 久久午夜亚洲精品久久| 久久精品国产a三级三级三级| 欧美日韩亚洲国产一区二区在线观看 | xxxhd国产人妻xxx| 国产一区二区三区在线臀色熟女 | 十八禁人妻一区二区| 丝袜在线中文字幕| 99热网站在线观看| 在线观看免费视频网站a站|