static corrections methods in the processing of deep reflection seismic...

10
Journal of Earth Science, Vol. 25, No. 2, p. 299–308, April 2014 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-014-0422-x Zhu, X. S., Gao, R., Li, Q. S., et al., 2014. Static Corrections Methods in the Processing of Deep Reflection Seismic Data. Journal of Earth Science, 25(2): 299–308, doi:10.1007/s12583-014-0422-x Static Corrections Methods in the Processing of Deep Reflection Seismic Data Xiaosan Zhu*, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu, Haiyan Wang Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China ABSTRACT: Statics are big challenges for the processing of deep reflection seismic data. In this paper several different statics solutions have been implemented in the processing of deep reflection seismic data in South China and their corresponding results have been compared in order to find proper statics solutions. Either statics solutions based on tomographic principle or combining the low-frequency components of field statics with the high-frequency ones of refraction statics can provide reasonable statics solutions for deep reflection seismic data in South China with very rugged surface topography, and the two statics solutions can correct the statics anomalies of both long spatial wavelengths and short ones. The surface-consistent residual static corrections can serve as the good compensations to the several kinds of the first statics solutions. Proper statics solutions can improve both qualities and reso- lutions of seismic sections, especially for the reflections of Moho in the upmost mantle. KEY WORDS: deep reflection seismic data, static correction, field static, refraction static, tomographic static, residual static correction. 1 INTRODUCTION Correcting near-surface velocity and elevation variations with statics is an essential step and static corrections are very important in the processing of land data, which can improve the qualities of subsequent processing steps and are related to the quality and resolution of final imaged section (Li L et al., 2011; Deere, 2009; Laak and Zaghloul, 2009; Li P et al., 2009a; Raef, 2009; Stein et al., 2009; Han et al., 2008; Vossen and Trampert, 2007; Yan et al., 2006; Criss and Cunningham, 2001). Static corrections are defined as (Cox, 1999; Sheriff, 1991): correc- tions applied to seismic data to compensate for the effects of variations in elevation, weathering thickness, weathering veloc- ity, or reference to a datum. The objective is to determine the reflection arrival times which would have been observed if all measurements had been made on a (usually) flat plane with no weathering or low-velocity material present. Hence it leads to the concept of surface-consistent corrections, which are de- pendent on the location of the source (or receiver) but are in- dependent of the source to receiver offset or time of the record data (Deere, 2009; Cox, 1999). There are many issues which are associated with the near surface and related with the variations of velocity and thickness in the near-surface layers. Field statics can compensate the data with some of the problems mentioned above and there are many papers which are focus on it (Luo et al., 2010; Li et al., 2009b; Huang et al., 2008). There are lots of static correction *Corresponding author: [email protected] © China University of Geosciences and Springer-Verlag Berlin Heidelberg 2014 Manuscript received August 21, 2013. Manuscript accepted January 15, 2014. methods based on seismic refraction principle, which can be used to resolve velocities of shallow layers using head waves, such as slope (or intercept) method (Knox, 1967), delay time method (Coppens, 1985), reciprocal method (Palmer, 1980), least square method (Chang et al., 2002; Simmons and Backus, 1992) and turn-rays method (Henley, 2009; Criss and Cun- ningham, 2001). Tomographic static correction methods have been developed by many researchers (Liu et al., 2010; Li et al., 2009b; Yordkayhun et al., 2009; Zhu et al., 2008; Taner et al., 1998) to obtain the static corrections, which use the tomo- graphic velocity models based on the first-arrival information to predict static corrections. These statics methods require a large number of rays going through the model areas evenly with different ray angles. Ray tomography methods have been used to build near-surface velocity models using first-arrival information and to estimate the static corrections (Zhang et al., 2009; Ke et al., 2007). Many residual static correction methods have been developed in order to compensate for the time delays during the last several decades, such as traveltime inversion based method (Hatherly et al., 1994), stack-power maximiza- tion method (Ronen and Claerbout, 1985), nonstationary resid- ual statics method (Henley, 2012) and sparsity maximization method (Gholami, 2013). In real world, there are many factors which cause the static corrections and residual static corrections difficult to be han- dled. These factors are including rugged surface acquisition topography, non-planar refractors, near-surface low-velocity layers, lateral variant velocities of weathering layers and varia- tions of underground water tables (Li et al., 2009b; Wang, 1999). Errors in static corrections lead to the losses of seismic resolutions, both temporal and spatial, and bring the difficulties and confusions during the interpretations of seismic sections. In this paper, we study the static correction methods using

Upload: danglien

Post on 31-Jan-2018

231 views

Category:

Documents


5 download

TRANSCRIPT

Page 1: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Journal of Earth Science, Vol. 25, No. 2, p. 299–308, April 2014 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-014-0422-x

Zhu, X. S., Gao, R., Li, Q. S., et al., 2014. Static Corrections Methods in the Processing of Deep Reflection Seismic Data. Journal of Earth Science, 25(2): 299–308, doi:10.1007/s12583-014-0422-x

Static Corrections Methods in the Processing of Deep Reflection Seismic Data

Xiaosan Zhu*, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu, Haiyan Wang

Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China

ABSTRACT: Statics are big challenges for the processing of deep reflection seismic data. In this paper several different statics solutions have been implemented in the processing of deep reflection seismic data in South China and their corresponding results have been compared in order to find proper statics solutions. Either statics solutions based on tomographic principle or combining the low-frequency components of field statics with the high-frequency ones of refraction statics can provide reasonable statics solutions for deep reflection seismic data in South China with very rugged surface topography, and the two statics solutions can correct the statics anomalies of both long spatial wavelengths and short ones. The surface-consistent residual static corrections can serve as the good compensations to the several kinds of the first statics solutions. Proper statics solutions can improve both qualities and reso-lutions of seismic sections, especially for the reflections of Moho in the upmost mantle. KEY WORDS: deep reflection seismic data, static correction, field static, refraction static, tomographic static, residual static correction.

1 INTRODUCTION

Correcting near-surface velocity and elevation variations with statics is an essential step and static corrections are very important in the processing of land data, which can improve the qualities of subsequent processing steps and are related to the quality and resolution of final imaged section (Li L et al., 2011; Deere, 2009; Laak and Zaghloul, 2009; Li P et al., 2009a; Raef, 2009; Stein et al., 2009; Han et al., 2008; Vossen and Trampert, 2007; Yan et al., 2006; Criss and Cunningham, 2001). Static corrections are defined as (Cox, 1999; Sheriff, 1991): correc-tions applied to seismic data to compensate for the effects of variations in elevation, weathering thickness, weathering veloc-ity, or reference to a datum. The objective is to determine the reflection arrival times which would have been observed if all measurements had been made on a (usually) flat plane with no weathering or low-velocity material present. Hence it leads to the concept of surface-consistent corrections, which are de-pendent on the location of the source (or receiver) but are in-dependent of the source to receiver offset or time of the record data (Deere, 2009; Cox, 1999).

There are many issues which are associated with the near surface and related with the variations of velocity and thickness in the near-surface layers. Field statics can compensate the data with some of the problems mentioned above and there are many papers which are focus on it (Luo et al., 2010; Li et al., 2009b; Huang et al., 2008). There are lots of static correction

*Corresponding author: [email protected] © China University of Geosciences and Springer-Verlag Berlin Heidelberg 2014 Manuscript received August 21, 2013. Manuscript accepted January 15, 2014.

methods based on seismic refraction principle, which can be used to resolve velocities of shallow layers using head waves, such as slope (or intercept) method (Knox, 1967), delay time method (Coppens, 1985), reciprocal method (Palmer, 1980), least square method (Chang et al., 2002; Simmons and Backus, 1992) and turn-rays method (Henley, 2009; Criss and Cun-ningham, 2001). Tomographic static correction methods have been developed by many researchers (Liu et al., 2010; Li et al., 2009b; Yordkayhun et al., 2009; Zhu et al., 2008; Taner et al., 1998) to obtain the static corrections, which use the tomo-graphic velocity models based on the first-arrival information to predict static corrections. These statics methods require a large number of rays going through the model areas evenly with different ray angles. Ray tomography methods have been used to build near-surface velocity models using first-arrival information and to estimate the static corrections (Zhang et al., 2009; Ke et al., 2007). Many residual static correction methods have been developed in order to compensate for the time delays during the last several decades, such as traveltime inversion based method (Hatherly et al., 1994), stack-power maximiza-tion method (Ronen and Claerbout, 1985), nonstationary resid-ual statics method (Henley, 2012) and sparsity maximization method (Gholami, 2013).

In real world, there are many factors which cause the static corrections and residual static corrections difficult to be han-dled. These factors are including rugged surface acquisition topography, non-planar refractors, near-surface low-velocity layers, lateral variant velocities of weathering layers and varia-tions of underground water tables (Li et al., 2009b; Wang, 1999). Errors in static corrections lead to the losses of seismic resolutions, both temporal and spatial, and bring the difficulties and confusions during the interpretations of seismic sections.

In this paper, we study the static correction methods using

Page 2: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Xiaosan Zhu, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu and Haiyan Wang

 

300

the deep reflection seismic data (Liu et al., 2010; Xu et al., 2005; Huang and Gao, 2001) along a long survey in South China (Fig. 1) and compare the statics solutions of different methods in order to find the proper statics solutions for the processing of deep reflection seismic data. The outline of the paper is the following. Firstly, we briefly summarize the theo-ries of field statics, refraction statics, tomographic statics and residual static corrections, separately. We then study on the statics solutions of deep reflection seismic data based on these methods mentioned above and compare the results of different statics solutions in detail. Finally, we summarize this study. 2 METHODOLOGY OF STATIC CORRECTIONS

In this section, the brief summarizations of four static cor-rection methods have been presented and these methods are field statics, refraction statics, tomographic statics and residual static corrections. 2.1 Field Statics

The source and receiver can be replaced on a reference datum with the datum static corrections according to the infor-mation of both elevation and near-surface velocity distribution from the uphole survey and the near-surface refraction data (Cox, 1999). The datum static corrections are including the weathering corrections for removing the effects of near-surface layers and the elevation corrections for moving from the base of these near-surface layers up to (or down to) a reference da-tum. The assumption of static corrections is that a simple time shift of an entire seismic trace which will yield the seismic record being observed if the geophone had been displaced ver-tically downward to the reference datum and the assumption is not strictly true in most cases. Strictly, the elevation correction can be used only in those areas there are no weathered layers and lateral velocity changes in low-velocity layers (Luo et al., 2010). If the velocity variations only affect the high-frequency components of the datum static corrections, then the elevation corrections can be used companying with residual static correc-tions. 2.2 Refraction Statics

Refraction methods allow us to derive estimates of the thicknesses and velocities of the near-surface layers by analyz-ing the first-breaks of the seismic records (Luo et al., 2010; Wu et al., 2009; Duan, 2006; Lin et al., 2006; Pan et al., 2003). According to the Huygens’ Principle, that is, every point on an advancing wavefront can be regarded as the source of a secon-dary wave and that a later wavefront is the envelope tangent to all the secondary waves (Cox, 1999). The important concept in seismic refraction is that when a seismic ray crosses a boundary between two formations of different velocities, then the ray is bent according to Snell’s law which defines that the sine of refracted angle is equal to the ratio of the velocities of the two formations. Therefore, the static correction based on refraction survey acquires the information of the first-arrival time of wavefield from refractor and the refractor velocity. Hence, there are two basic conditions for refraction survey, that is, a relative stable refraction interface between the two formations and the acknowledged near-surface velocity distribution (Bridle

and Aramco, 2009; Liu, 1998). Applying the static corrections based on refraction survey

can ensure structural integrity in the processed section. Refrac-tion statics are effective for correcting long spatial wavelength anomalies and compensating for the weathering layers. Actually, refraction statics are also effective against short spatial wave-length anomalies (Liu, 1998). 2.3 Tomographic Statics

Tomographic statics are commonly used during the proc-essing of seismic data, especially in the areas with rapid veloc-ity variations in laterally (Hao et al., 2011; Luo et al., 2010; Han et al., 2008; Wang, 2005; Yang et al., 2005). The definition of tomography is that (Sheriff, 1991) a method for finding the velocity and reflectivity distribution from a multitude of obser-vations using combinations of source and receiver locations. The tomographic inversion approaches use the first arrival in-formation of wavefront to inverse the velocity distribution of near-surface without the assumption of layer structure in order to produce a near-surface velocity model which best fits the observed minimum arrival times. Space is divided into cells and the data are expressed as line integrals along raypaths through the cells. Iterated adjusting and updating the near-surface velocity model, until the differences between arri-val times of model and those of the observed data reach ac-ceptable levels or are unchanged between iterations (Becerra et al., 2009; Henley, 2009; Li et al., 2009b; Vossen and Trampert, 2007; Chang et al., 2002). Tomographic methods include the algebraic reconstruction technique (ART) (Henley, 2009), the simultaneous reconstruction technique (SIRT) (Aster et al., 2005; Emily and Bradford, 2002) and Gauss-Seidel method (Taner et al., 1998).

The static solutions based on tomography principle need a large number of different ray paths which go through each of cells with a wide-angle coverage and constrains of indirect regularization during the inversion. The methods provide proper corrections for long and middle spatial wavelength components under most of situations with rugged surface to-pography and rapidly changed velocities in near-surface layers. However, there are still some disadvantages of static correc-tions based on tomographic techniques and the uncertainties in tomographic velocity models have also been qualified from a 2D seismic line acquired in Colombia through a variety of nu-merical techniques (Becerra et al., 2009).

2.4 Residual Static Corrections

The residual static corrections are time shifts applied to traces in order to compensate for time delays and the statics model is a function of time and space (Henley, 2012; Li et al., 2011; Sheriff, 1991). Residual static corrections are defined as a subset of the static corrections (Cox, 1999). Data-smoothing statics methods assume that patterns of irregularities which most events have in common result from near-surface varia-tions and hence static correction trace shifts should be such as to minimize those irregularities. Sheriff (1991) describes that the concept of static correction is the assumption that a simple time shift of an entire seismic trace will yield the seismic re-cords which would have been observed if the geophones had

Page 3: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Static Corrections Methods in the Processing of Deep Reflection Seismic Data

 

301

been displaced vertically downward to a reference datum. The time shift approximation means that static corrections are surface-consistent and independent of reflection times and trace offsets.

Due to the near-surface model for statics solutions is a simplification of the geology resulting in a tradeoff between thicknesses and velocities which result in inexact static correc-tions and these corrections are the approximations for more complex problems, the applications of field statics, or refraction statics, never leave the seismic data completely free of static anomalies (Yin et al., 2004; Jing, 2003). Therefore, it’s definite necessary that the residual static anomalies should be handled properly. In reality, residual static anomalies are compensated for using statistical correlation techniques. Usually the residual static corrections are extracted from integrated seismic sections and designed to correct small inaccuracies in the near-surface velocity model, which seek to enhance the qualities and resolu-tions of stacked seismic sections. 3 STATIC CORRECTIONS OF DEEP REFLECTION SEISMIC DATA

The deep reflection seismic data used in this paper are along a two-dimensional survey in South China (Fig. 1). There are very rugged surface topographies and rapid variant veloci-ties of near-surface layers in both laterally and vertically due to the variations of compaction and lithology along the survey. The acquisition line is very long (around 550 km) and the ele-vation along the acquisition line of survey is shown in Fig. 2a. It’s difficult to deal with the static solutions of the deep reflec-tion seismic data from this area and the velocity model of near-surface layers along the survey (Fig. 2b) is obtained using a ray-tracing method during the procedure of tomographic stat-ics. If static corrections are not properly handled during the processing of seismic data, then a whole catalog of problems will affect the interpretations of the seismic sections, including lines with variable elevations, false structural anomalies re-

maining in the sections, false events being created out of noises and the data qualities not being optimized. Therefore, proper statics solutions are definite desirable for obtaining high-resolution sections which can be used for both strati-graphic and lithologic interpretations. As for the deep reflection seismic data, proper statics solutions are very important in or-der to obtain the final clear and accurate images of the crust and upper mantle.

Several static correction methods including field statics, refraction statics, tomographic statics and in the wave of resid-ual statics are used during the data processing of deep reflection seismic data. We also combine the low-frequency components of field statics solutions with the high-frequency ones of refrac-tion statics solutions in order to obtain more reasonable statics solutions and this procedure is shown in Fig. 3a. In order to compare the results of these different statics solutions, a raw shot profile and those profiles applied with different statics solutions are shown in Figs. 3b–3f. The static corrections of these methods for all receivers and sources of the survey are shown in Figs. 4a and 4b, respectively. Four common middle point (CMP) stacked profiles of deep reflection seismic data corresponding to applying these statics solutions, that is, field statics, refraction statics, tomogrpahic statics and combining of the former two statics solutions, are shown in Figs. 5a, 5b, 6a and 6b, respectively. Comparing these solutions shown in Figs. 3c–3f, the results of refraction statics (Fig. 3d) are slightly bet-ter than those of field ones (Fig. 3c). Both the results of tomo-graphic statics solutions (Fig. 3e) and those of the combining solutions of field statics with those of refraction ones (Fig. 3f) can provide good qualities of shot profiles, and the reflection events have better continuities than those in the other two pro-files. The similar conclusions can be derived from the CMP stacked profiles after applied the four kinds of statics solutions mentioned above which are shown in Figs. 5 and 6. Comparing the three kinds of statics solutions using separately for both receivers (Fig. 4a) and sources (Fig. 4b), the field static

Figure 1. Location map showing the deep reflection seismic profile in South China. The red line indicates the location of sur-vey and its length is around 550 km. The red square at the right upper corner shows the location of study area and the black lines indicate the locations of faults. The survey is a symmetric survey with 700 receivers distributed at two sides of source and the total source number is 2 269, receiver spacing is every 40 m and source spacing is every 280 m, near offset is 140 m from the source, 7 499 time gates are recorded with 4 ms spacing.

Page 4: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Xiaosan Zhu, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu and Haiyan Wang

 

302

Figure 2. (a) Elevation along the survey line of the deep reflection seismic profile in South China. (b) Velocity model of near-surface layers along the survey which is obtained using ray-tracing method (the black line indicates the ray bottom and the velocity below the line is not creditable).

Figure 3. (a) Static corrections of all the receivers of the survey by combining the low-frequency components of field statics solutions with the high-frequency components of refraction ones; (b) raw shot profile; (c) field static corrections applied; (d) refraction static corrections applied; (e) tomographic static corrections applied; (f) applying the combining low spatial wave-length components of field statics solutions with the high ones of refraction statics solutions. The white rectangles show the areas with great improvements after applied statics solutions.

Page 5: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Static Corrections Methods in the Processing of Deep Reflection Seismic Data

 

303

Figure 4. Profiles of static corrections of all the receivers (a) and sources (b) of the survey for field statics, refraction statics, tomographic statics and combining field statics with refraction ones.

Figure 5. CMP stacked sections illustrating the results of field static corrections (a) and those of refraction static corrections (b).

Page 6: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Xiaosan Zhu, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu and Haiyan Wang

 

304

Figure 6. CMP stacked sections illustrating the results of tomographic static corrections (a) and those of combining the low-frequency components of field statics solutions with the high-frequency ones of refraction statics solutions (b). corrections are slightly big and the refraction ones are some-what small, however the tomographic ones are the best solu-tions among them. The combining statics solutions of the for-mer two can also provide reasonable solutions in this case.

Implement of field statics is very fast and need only small amount of computation time. During the processing of seismic data, field static corrections are usually served as a basic stan-dard of quality control in order to obtain some basic informa-tion for both the parameters of static corrections and its pre-liminary stack section of the deep reflection seismic data. The statics solutions based on refraction principles work well in the region with mild topography and well behaved weathering layers. However, the refraction model does not match the geo-logic reality of complex terrains in most cases and the refrac-tion statics cannot properly handle the conditions with inverse velocity distribution layers where the low-velocity layers locate under the high-velocity ones and hidden layers which are too thin to be recognized. It’s a good choice that combing the ad-vantages of both field statics and refraction ones, which is shown in Fig. 3a. The reason is that the combining statics solu-tions can correct the statics anomalies of both long spatial wavelengths and short ones. However, the problem for this procedure is that it's difficult to handle the ratios of the static corrections of field statics versus those of refraction ones be-cause the velocity distributions of near-surface layers varying strongly. The statics solutions based on tomographic techniques

can provide proper corrections for long and middle spatial wavelengths components in most cases and work well for those areas with rapidly changed thicknesses, strong velocity varia-tions in laterally and vertically of near-surface layers and com-plex subsurface geology as long as enough first-arrival infor-mation from the relative small offset has been used during in-version procedure. Tomographic statics can be used to obtain the model of near-surface low-velocity layers (Fig. 2b) and it is very useful for the velocity updating in the subsequent proc-essing steps. However, there are still some shortcomings of this method. The solution of tomographic inversion is not unique and unstable usually, and it is sensitive to the initial velocity model and the picking accuracies of first-arrivals. Zhang et al. (2005) have developed a hybrid optimization inversion method to calculate large statics by integrating stack-power maximiza-tion, simulated annealing and genetic algorithm in complex terrains. In this method, large statics are corrected using a spe-cial smooth filtering operator which can eliminate the pseudo-static corrections from long spatial wavelengths com-ponents to short ones iteratively. Therefore, the method has some combined abilities of field statics, refraction statics and tomographic statics.

Due to the complex geological structures and most of static correction methods are based on simplified models, it’s difficult to obtain the accurate velocity model of near-surface layers no matter what kinds of methods being used, and there

Page 7: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Static Corrections Methods in the Processing of Deep Reflection Seismic Data

 

305

are always some residual statics anomalies remaining in seis-mic sections. Residual static corrections can enhance the quali-ties of stacked traces using statistical correlation methods after applying those first statics solutions (i.e., field statics, refraction statics, tomographic statics and combining of the former two statics solutions). In reality, the residual static corrections are

used iteratively for obtaining close to free of statics anomalies of stacked traces. In this implement, we iteratively apply the residual static corrections three times. The first and third resid-ual static corrections for all the receivers and sources of the survey are shown in Figs. 7a–7d. From these figures, we see that both the residual statics solutions of receivers and source

Figure 7. Profiles of the first residual statics solutions of both receivers (a) and sources (b) after applying those first static corrections (i.e., field static corrections, refraction static corrections, tomographic static corrections and combining the low-frequency components of field statics solutions with the high-frequency ones of refraction statics solutions). Profiles of the third residual statics solutions of both receivers (c) and sources (d) of the survey after applying the second residual statics solutions.

Page 8: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Xiaosan Zhu, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu and Haiyan Wang

 

306

Figure 8. CMP stacked sections illustrating the changes in apparent structure after residual static corrections applied. (a) Field static corrections applied; (b) residual static corrections applied after field static corrections; (c) refraction static corrections applied; (d) residual static corrections applied after refraction static corrections. The white rectangles show the areas with great im-provements of the reflections of Moho after applied re-sidual statics solutions.

Figure 9. CMP stacked sections illustrating the changes in apparent structure after residual static corrections applied. (a) Tomographic static corrections applied; (b) residual static corrections applied after tomographic static corrections; (c) combining the statics solutions of field statics solutions with those of refraction statics solu-tions applied; (d) residual static corrections applied after the combining statics solutions. The white rectangles show the areas with great improvements of the reflec-tions of Moho after applied residual statics solutions.

Page 9: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Static Corrections Methods in the Processing of Deep Reflection Seismic Data

 

307

become smaller and smaller with the iterated applying the re-sidual statics corrections. Therefore, we can reduce the statics anomalies in the deep reflection seismic sections iterated using the same procedure.

In order to compare the results of residual statics solutions, four CMP stacked profiles of deep reflection seismic data are shown in Figs. 8 and 9. Due to the reflections of Moho are very important for deep reflection seismic profiles, these profiles are shown corresponding to the special time sections of these re-flections of Moho in CMP stacked profiles. Figures 8 and 9 show that most of the reflections of Moho in those profiles applied residual statics solutions have better continuities than those without applying this procedure, and the resolutions of them have been greatly improved after applying residual static corrections. The residual static corrections are represented usu-ally by the short spatial wavelength components of the profiles and not the long ones (Li et al., 2011). However, the serious shortcoming of residual static corrections is that there will exit “cycle skipping” in the sections when the statics errors are greater than half the length of seismic wavelet, which looks like faults in the data and may be misaligned by a whole cycle of the seismic wavelet.

Theoretically, all the first static corrections (i.e., field stat-ics, refraction statics, tomographic statics and combining the former two) and residual static corrections should be treated as not only frequency dependent time shifts but also time and phase components. Therefore, the challenging rugged surface topography and large magnitude statics need the processing procedures with iterative statics calculations, noises attenuation and the velocity updating of near-surface model for obtaining high-resolution seismic section. 4 CONCLUSIONS

In this paper, both the field statics solutions and refraction ones used separately in the processing of deep reflection seis-mic data along the long survey in South China cannot derive reasonable statics solutions anymore due to the rugged surface topography, low-velocity layers and the velocities of near-surface layers varying strongly in laterally and vertically along the survey. However, statics solutions based on tomo-graphic principle can provide proper solutions for this kind of situation. Combining the low-frequency components of field statics solutions with the high-frequency ones of refraction statics solutions can also provide reasonable solutions for the deep reflection seismic data in South China. The surface- consistent residual static corrections are good compensations to the procedures of the first statics solutions studied in the paper and can leave the deep reflection seismic data close to free of statics anomalies. ACKNOWLEDGMENTS

This work was supported by the Foundation of Institute of Geology, Chinese Academy of Geological Sciences (No. J1315), the 3D Geological Mapping Project (No. D1204) and the SinoProbe-02 project of China. The authors thank Hongqiang Li and Gong Deng for providing wonderful com-ments and suggestions.

REFERENCES CITED Aster, R. C., Borchers, B., Thurber, C. H., 2005. Parameter

Estimation and Inverse Problems. Elsevier Academic Press

Becerra, C., Agudelo, W., Guevara, S., 2009. Uncertainty Analysis in Statics Corrections Obtained by Tomographic Inversion: Application in a Mountainous Zone in Catatumbo (Colombia). The Leading Edge, 28(2): 212–215

Bridle, R., Aramco, S., 2009. Delay-Time Refraction Methods Applied to a 3D Seismic Block. The Leading Edge, 28(2): 228–237

Chang, X., Liu, Y. K., Wang, H., et al., 2002. 3-D Tomography Static Correction. Geophysics, 67(4): 1275–1285

Coppens, F., 1985. First Arrival Picking on Common-Offset Trace Collections for Automatic Estimation of Static Cor-rections. Geophys. Prosp., 33: 1212–1231

Cox, M., 1999. Static Corrections for Seismic Reflection Sur-veys: Society of Exploration Geophysicists Publ., Tulsa, Oklahoma. 1–531

Criss, D. E., Cunningham, D., 2001. Turning-Ray Tomography for Statics Solution. EAGE 63rd Conference and Techni-cal Exhibition-Amsterdam, The Netherlands. 11–15

Deere J., 2009. Introduction to This Special Section—Statics. The Leading Edge, 28(2): 190–191

Duan, Y. Q., 2006. Residual Static Corrections Based on Re-fraction Survey. OGP, 41(1): 32–35 (in Chinese with Eng-lish Abstract)

Emily, A. H., Bradford, J. H., 2002. Ground-Penetrating-Radar Reflection Attenuation Tomography with an Adaptive Mesh. Geophysics, 75(4): WA251–WA261

Gholami, A., 2013. Residual Statics Estimation by Sparsity Maximization. Geophysics, 78(1): 11–19

Han, X. L., Yang, C. C., Ma, S. H., et al., 2008. Static of To-mographic Inversion by First Breaks in Complex Areas. Progress in Geophysics, 23(2): 475–483 (in Chinese with English Abstract)

Hao, J., Yang, R. J., Wu, J., et al., 2011. Processing of Static Correction Problems of Seismic Data in the Complex Sur-face. Complex Hydrocarbon Reservoirs, 4(3): 34–37 (in Chinese with English Abstract)

Hatherly, P. J., Urosevic, M., Lambourne, A., et al., 1994. A Simple Approach to Calculating Refraction Statics Cor-rections. Geophysics, 59(1): 156–160

Henley, D. C., 2009. Raypath interferometry: Statics in Diffi-cult Places. The Leading Edge, 28(2): 202–205

Henley, D. C., 2012. Interferometric Application to Static Cor-rections. Geophysics, 77(1): Q1–Q13

Huang, L. Y., Gao, R., 2001. The Sub-Dataset of Deep Reflec-tion Seismic Data. Acta Geoscientia Sinica, 22(6): 491–496 (in Chinese with English Abstract)

Huang, M. Z., Feng, Z. Y., Zhou, D. T., 2008. Directly Iterated Static Corrections Method in Offset Domain and Its Ap-plication. Progress in Exploration Geophysics, 31(2): 122–128 (in Chinese with English Abstract)

Jing, X. L., 2003. Two Steps Solution Method for Big Residual Static Corrections. OGP, 38(1): 22–26 (in Chinese with English Abstract)

Page 10: Static Corrections Methods in the Processing of Deep Reflection Seismic ...en.earth-science.net/PDF/20140410092534.pdf · Static Corrections Methods in the Processing of Deep Reflection

Xiaosan Zhu, Rui Gao, Qiusheng Li, Ye Guan, Zhanwu Lu and Haiyan Wang

 

308

Ke, B., Zhang, J., Chen, B., et al., 2007. Fat-Ray First Arrival Seismic Tomography and Its Application. 77th Annual International Meeting, SEG. 3354–3358 (Expanded Ab-stracts)

Knox, W. A., 1967. Multilayer Near-Surface Refraction Com-putations. In: Musgrave, A. W., ed., Seismic Refraction Prospecting. Soc. Expl. Geophys., 197–216

Laake, A., Zaghloul, A., 2009. Estimation of Static Corrections from Geologic and Remote-Sensing Data. The Leading Edge, 28(2): 192–196

Li, P., Zhou, H., Yan, Z., 2009a. Deformable Layer Tomostatics: 2D Examples in Western China. The Leading Edge, 28(2): 206–210

Li, P., Feng, Z., Li, Z., et al., 2009b. Static Correction Tech-nology and Applications in Complex Areas of Western China. The Leading Edge, 28(2): 1384–1386

Li, L., Chen, X. J., Jing, X. L., 2011. Multiscale Inversion Al-gorithm for Seismic Residual Static Correction and Its Application. Xinjing Petroleum Geology, 32(4): 402–405 (in Chinese with English Abstract)

Lin, B. X., Sun, J. M., Xu, Y., et al., 2006. Discussion on Sev-eral Common Static Correction Methods. Geophysical Prospecting for Petroleum, 45(4): 367–373 (in Chinese with English Abstract)

Liu, L. S., 1998. Constrained First-Arrival Pickup and First-Break Residual Static Correction. OGP, 33(5): 604–610 (in Chinese with English Abstract)

Liu, J. K., Kuang, C. Y., Gao, R., et al., 2010. Data Processing Test and Research on the Deep Seismic Reflection Profile in Polymetallic Deposits Area: Taking an Example of Lu-zong Ore Concentrated Area. Acta Petrological Sinica, 26(9): 2561–2576 (in Chinese with English Abstract)

Luo, Y. W., Yang, J., Duan, W. X., et al., 2010. Comparing Be-tween Several Static Corrections Methods. Petroleum In-struments, 24(5): 41–43 (in Chinese with English Ab-stract)

Pan, H. X., Fang, W. B., Wu, Y. S., et al., 2003. An Improved Relative Refraction Statics Technique. Geophysical Pros-pecting for Petroleum, 42(2): 208–211 (in Chinese with English Abstract)

Palmer, D., 1980. The Generalized Reciprocal Method of Seismic Refraction Interpretation. Society of Exploration Geophysicists. 1–104

Raef, A., 2009. Land 3D-Seismic Data: Preprocessing Quality Control Utilizing Survey Design Specifications, Noise Properties, Normal Moveout, First Breaks, and Offset. Journal of Earth Science, 20(3): 640–648, doi:10.1007/s12583-009-0053-9

Ronen, J., Claerbout J., 1985. Surface-Consistent Residual Statics Estimation by Stack-Power Maximization, Geo-physics, 50(2): 2759–2767

Sheriff, R. E., 1991. Encyclopedic Dictionary of Exploration

Geophysics. Society of Exploration Geophysicists.1–323 Simmons, J. L., Backus, M. M., 1992. Linearized Tomographic

Inversion of First-Arrival Times. Geophysics, 57: 1482–1492

Stein, J. A., Langston, T., Larson, S. E., 2009. A Successful Statics Methodology for Land Data. The Leading Edge, 28(2): 222–226

Taner, M. T., Wagner, D. E., Baysal, E., et al., 1998. A Unified Method for 2-D and 3-D Refraction Statics. Geophysics, 63: 260–274

Vossen, R. V., Trampert, J., 2007. Full-Waveform Static Cor-rections Using Blind Channel Identification. Geophysics, 72(4): U55–U66

Wang, S. D., 2005. Static Corrections of Complex Topography Based Wave Equation Datuming. OGP. 40(1): 31–34 (in Chinese with English Abstract)

Wang, J. H., 1999. Thanking about the Normal Moveout Cor-rections and Static Corrections. OGP, 34(Suppl.): 18–26 (in Chinese with English Abstract)

Wu, K. F., Zhang, X. Q., Zheng, G. Y., et al., 2009. Review of Converted-Wave Statics Correction Method Based on Body-Wave. Chinese Journal of Engineering Geophysics, 6(6): 768–774 (in Chinese with English Abstract)

Xu, M. C., Gao, J. H., Rong, L. X., et al., 2005. Exploration Technique Combing Tomographic Statics with High-Resolution Seismic Data. Geology and Prospecting, 41(4): 83–87 (in Chinese with English Abstract)

Yan, X., Zhong, G. F., Li, Q. Y., et al., 2006. Stratal Carbonate Content Inversion Using Seismic Data and Its Applica-tions to the Northern South China Sea. Journal of China University of Geosciences, 17(4): 320–325

Yang, W. J., Duan, Y. Q., Jiang, W. C., et al., 2005. Tomo-graphic Statics. Geophysical and Geochemical Explora-tion, 29(1): 41–43 (in Chinese with English Abstract)

Yin, C., Xiong, X. J., Zhang, B. L., et al., 2004. The Study of Using the Fourth Accumulated Component in Residual Static Corrections. Gas Industry, 24(12): 48–50 (in Chi-nese with English Abstract)

Yordkayhun, S., Tryggvason, A., Norden, B., et al., 2009. 3D Seismic Traveltime Tomography Imaging of the Shallow Subsurface at the CO2SINK Project Site, Ketzin, Ger-many. Geophysics, 74(1): G1–G15

Zhang, J., Zhao, B., Zhou, H., 2009. Fat Ray Tomography with Optimal Relaxation Factor. 79th Annual International Meeting, SEG. 4044–4048 (Expanded Abstracts)

Zhang, Z. J., Lin, Y. H., Liu, E. R., 2005. Large Static Correc-tion Using a Hybrid Optimization Method in Complex Terrains: Some Experience Learnt From China. Journal of Seismic Exploration, 13(4): 337–352

Zhu, X., Valasek, P., Roy, B., et al., 2008. Recent Applications of Turning-Ray Tomography. Geophysics, 73(5): VE243–VE254