<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">IJG</journal-id><journal-title-group><journal-title>International Journal of Geosciences</journal-title></journal-title-group><issn pub-type="epub">2156-8359</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ijg.2017.81003</article-id><article-id pub-id-type="publisher-id">IJG-73691</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Optimization on the High Precision 3D Seismic Acquisition Parameters in Leijia Tight Sandstone Oil Area
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Haibo</surname><given-names>Wang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wei</surname><given-names>Zhang</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bing</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Qiwei</surname><given-names>Zou</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Guanhua</surname><given-names>Xiao</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Meibin</surname><given-names>Tang</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hong</surname><given-names>Zhu</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wenjing</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>School of Geophysics and Information Technology, China University of Geosciences, Beijing, China</addr-line></aff><aff id="aff3"><addr-line>Liaohe Oil field, CNPC, Panjin, China</addr-line></aff><aff id="aff2"><addr-line>BGP Inc., CNPC, Panjin, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>wanghaibo004@qq.com(HW)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>15</day><month>01</month><year>2017</year></pub-date><volume>08</volume><issue>01</issue><fpage>30</fpage><lpage>45</lpage><history><date date-type="received"><day>October</day>	<month>31,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>January</month>	<year>19,</year>	</date><date date-type="accepted"><day>January</day>	<month>22,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The main targets of seismic exploration research in Leijia carbonatite tight sandstone oil area of Liaohe depression are thin reservoirs prediction and minor faults identification, which is one of the important representatives of complex exploration objects in Liaohe depression. High precision 3D seismic exploration has significantly improved the ability of thin reservoirs prediction and minor faults identification of this area. Reducing the cost of high precision 3D seismic exploration through optimizing the acquisition parameters is very important for the next step exploration and development of Liaohe depression and similar areas. Based on high precision 3D seismic acquisition data in Leijia tight sandstone oil area, multiple sub-geometries are obtained with different bin sizes, different folds, different aspect ratio, different line intervals by extracting receiver points and shot points, and PSTM processing is performed respectively, obtained PSTM datasets of the sub-geometries, extract time slices, amplitude slices along the layer, coherent slices and so on. We evaluate the data results of the sub-geometries from the aspects of signal-to-noise ratio, thin reservoirs resolution, acquisition footprint and so on. Considering the exploration cost and data effect of each sub-geometry, the optimal direction of the main parameters of high precision seismic exploration in Liaohe depression is put forward, and the acquisition effect of adjacent area by the optimized parameters is given.
 
</p></abstract><kwd-group><kwd>High Precision 3D</kwd><kwd> Acquisition Parameters Optimization</kwd><kwd> Fold</kwd><kwd> Bin Size</kwd><kwd> Aspect Ratio</kwd><kwd> Line Interval</kwd><kwd> Acquisition Footprint</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>“Two Wide and One High” (Wide azimuth, wide band and high density) seismic exploration technology is a high precision seismic exploration push out for solving complicated area exploration technology in recent years, which has carried out trial in the Caspian Sea region abroad [<xref ref-type="bibr" rid="scirp.73691-ref1">1</xref>] , and in Junggar basin, Tarim basin, Liaohe depression of China [<xref ref-type="bibr" rid="scirp.73691-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.73691-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.73691-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.73691-ref5">5</xref>] , and achieved good results. But sample density of the seismic acquisition geometry becomes higher and higher, such as the bin size becomes smaller and smaller, the line interval becomes smaller and smaller, and the fold becomes higher and higher. It is generally believed that the high sampling density will certainly be better for imaging, which may result in a waste of cost. It is necessary to study whether the continuous enhancement of the geometry parameters has much effect on the imaging of complex geological targets and whether there is the possibility of weakening.</p><p>Domestic and foreign scholars using the actual acquired seismic data to do degenerative processing analysis to study the impact of 3D acquisition parameters on the data imaging are rare. Wloszczowski [<xref ref-type="bibr" rid="scirp.73691-ref6">6</xref>] carried out degradation processing about maximum offset, folds using a real data from the North Sea, to evaluate the impact of geometry parameters changing on the influence of signal to noise ratio, wavelet consistency, target attribute maps and deep imaging. P. I. Pecholcs [<xref ref-type="bibr" rid="scirp.73691-ref7">7</xref>] , decimating the line intervals of the original volume from 100 m to 200 m and 500 m using a real high density data from Saudi Arabia, clearly shows how seismic signal to noise ratio and resolution improves with finer spatial sampling. Yun Ling [<xref ref-type="bibr" rid="scirp.73691-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.73691-ref9">9</xref>] carried out degradation processing and comprehensive interpretation using a full azimuth acquisition data from Junggar basin of western China, indicates that the wide azimuth (aspect ratio is greater than 0.5) has better spatial imaging resolution than narrow azimuth on lithologic exploration. Wei Zhang [<xref ref-type="bibr" rid="scirp.73691-ref10">10</xref>] carried out degradation processing about bin size, the number of shot points and array length from a high resolution acquisition data from Junggar basin of western China, to evaluate the impact of these parameters on the final stack profile. Ying Xu [<xref ref-type="bibr" rid="scirp.73691-ref11">11</xref>] and Yunhong Zhu [<xref ref-type="bibr" rid="scirp.73691-ref12">12</xref>] carried out degradation processing from high precision acquisition data from Tahe oilfiled and Tazhong oilfiled of western China respectively, they both proposed suggestions on key acquisition parameters for the high-precision 3D seismic exploration of their own area.</p><p>Bohai Bay Basin in eastern China is an important oil and gas production base. The sedimentary tectonics of Bohai Bay Basin is totally different from western China, and the seismic and geological conditions vary greatly. The literatures on degradation processing of high precision seismic data acquired in eastern China were not found. Based on the high precision 3D data of Leijia area in Liaohe depression of Bohai Bay Basin, this paper carried out degradation processing such as bin size, fold, aspect ratio, line interval, etc, systematically compared the difference of signal to noise ratio, thin layer imaging, fault imaging, data continuity caused by the geometry parameter variations from prestack time migration profile, time slice, amplitude slices along the layer, coherence slices, etc., put forward the optimization direction of high precision seismic exploration key acquisition parameters in Liaohe depression, provides the guidance for the optimization of high precision seismic acquisition in Bohai Bay Basin similar regions.</p></sec><sec id="s2"><title>2. Trials and Analysis</title><p>We performed degradation processing of bin size, fold, aspect ratio and line interval on the base of Leijia high precision acquisition data, discusses the acquisition parameters with high cost performance under the premise of satisfying the geological task requirements, economic and technological integration requirements, put forward the optimization of the geometry, to guide the next step acquisition solutions development.</p><sec id="s2_1"><title>2.1. Original Geometry Analysis</title><p>Leijia high precision 3D seismic acquisition geometry is 32L &#215; 10S &#215; 352T orthogonal (<xref ref-type="fig" rid="fig1">Figure 1</xref>), bin size: 10 m &#215; 10 m, fold: 256, receiver line interval/shot line interval: 200 m/220 m, the maximum offset: 4743 m; aspect ratio: 0.91. Shooting with charge.</p><p>This geometry has the characteristics of small bin size, high fold, large aspect ratio, high fold, the azimuth and offset is evenly distributed, and can meet the requirement of degradation processing of bin size, fold, aspect ratio and so on. Through analysis the influence of different geometry parameters on the effect of prestack time migration processing to optimize the geometry parameters in the complex area of Liaohe depression. In order to meet the requirements of geological task under the premise of saving seismic acquisition cost, to achieve economic and technological integration.</p></sec><sec id="s2_2"><title>2.2. Degradation Processing Solutions</title><p>According to the characteristics of the 3D acquisition geometry in the experimental area, design four sets of 6 trial geometries (<xref ref-type="table" rid="table1">Table 1</xref>). S0 is the original geometry. S1 to S5 are sub-geometries.</p><p>The first group of geometry is composed of S0 and S1 in the table, which used</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Geometry template (a) and rose diagram (b)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x2.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Degradation processing trial geometry parameters table</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Parameters Data set</th><th align="center" valign="middle" >geometry</th><th align="center" valign="middle" >Bin size (m)</th><th align="center" valign="middle" >receiver line interval (m)</th><th align="center" valign="middle" >shot line interval (m)</th><th align="center" valign="middle" >fold</th><th align="center" valign="middle" >trace density (10<sup>4</sup>/km<sup>2</sup>)</th><th align="center" valign="middle" >aspect ratio</th></tr></thead><tr><td align="center" valign="middle" >S0 (original)</td><td align="center" valign="middle" >32 &#215; 10S &#215; 352T</td><td align="center" valign="middle" >10 &#215; 10</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >256</td><td align="center" valign="middle" >256</td><td align="center" valign="middle" >0.91</td></tr><tr><td align="center" valign="middle" >S1</td><td align="center" valign="middle" >32L &#215; 5S &#215; 176T</td><td align="center" valign="middle" >20 &#215; 20</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >256</td><td align="center" valign="middle" >64</td><td align="center" valign="middle" >0.91</td></tr><tr><td align="center" valign="middle" >S2</td><td align="center" valign="middle" >32L &#215; 10S &#215; 352T</td><td align="center" valign="middle" >10 &#215; 10</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >440</td><td align="center" valign="middle" >128</td><td align="center" valign="middle" >128</td><td align="center" valign="middle" >0.91</td></tr><tr><td align="center" valign="middle" >S3</td><td align="center" valign="middle" >16L &#215; 10S &#215; 352T</td><td align="center" valign="middle" >10 &#215; 10</td><td align="center" valign="middle" >400</td><td align="center" valign="middle" >440</td><td align="center" valign="middle" >64</td><td align="center" valign="middle" >64</td><td align="center" valign="middle" >0.91</td></tr><tr><td align="center" valign="middle" >S4</td><td align="center" valign="middle" >28L &#215; 10S &#215; 352T</td><td align="center" valign="middle" >10 &#215; 10</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >224</td><td align="center" valign="middle" >224</td><td align="center" valign="middle" >0.80</td></tr><tr><td align="center" valign="middle" >S5</td><td align="center" valign="middle" >16L &#215; 10S &#215; 352T</td><td align="center" valign="middle" >10 &#215; 10</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >220</td><td align="center" valign="middle" >128</td><td align="center" valign="middle" >128</td><td align="center" valign="middle" >0.51</td></tr></tbody></table></table-wrap><p>for the bin size trial. Compared the PSTM processing results of different bin size (10 m &#215; 10 m, 20 m &#215; 20 m) with the same fold (256), to determine the appropriate bin size of this area. The bin size of S1 is 20 m &#215; 20 m, extract one for two shot points in the crossline direction, and extract one for two receiver points in the inline direction, that is the shot point interval and receiver point interval from 20 m into 40 m.</p><p>The second group of geometries is composed of S0, S2 and S3 in <xref ref-type="table" rid="table1">Table 1</xref>, which used for fold trial. Compared the PSTM processing results of different fold (256, 128, 64) with the same bin size (10 m &#215; 10 m), to determine the appropriate fold of this area. S2 is based on the original geometry, extract one for two shot points in the inline direction, and the receiver point remains unchanged, that is the shot line interval from 220 m into 440 m. S3 is based on S2, extract one for two receiver line, that is the shot line interval from 200 m into 400 m.</p><p>The third group of geometries is composed of S1 and S3 in <xref ref-type="table" rid="table1">Table 1</xref>. Compared the PSTM processing results of different fold (256, 64) and different bin size (20 m &#215; 20 m, 10 m &#215; 10 m) with the same trace density (64 &#215; 10<sup>4</sup>/km<sup>2</sup>), to determine the appropriate fold of this area. S2 is based on the original geometry, extract one for two shot points in the inline direction and the receiving point remains unchanged, that is the shot line interval from 220 m into 440 m. S3 is base on S2, extract one for two receiver lines, that is the receiver line interval from 200 m into 400 m. The aim was to determine whether it would be more appropriate to use large bin size and high fold or small bin size and small fold acquisition under essentially the same workload.</p><p>The fourth group geometries is composed of S0, S4 and S5 in <xref ref-type="table" rid="table1">Table 1</xref>, and is used for aspect ratio trial, comparing with the same surface element (10 m &#215; 10 m), and the processing effect of different aspect ratio (0.91, 0.80, 0.51).</p><p>The fifth group geometries is composed of S0, S4 and S5 in <xref ref-type="table" rid="table1">Table 1</xref>, and is used to compared different receiver lines with the same bin size.</p></sec><sec id="s2_3"><title>2.3. Degradation Processing Flow</title><p>Trial processing and analysis flow shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, mainly divided into fine</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Trial processing and analysis flow</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x3.png"/></fig><p>prestack CMP gather processing and geometries degradation processing. The prestack processing is based on all the data of S0 geometry (bin size 10 m &#215; 10 m, fold 256) for a detailed processing parameters and flow testing. Include the following steps, (1) amplitude preserved denoising, cross domain adaptive surface wave attenuation to remove surface wave , fast and parallel frequency division to remove the abnormal amplitude, coherent noise attenuation to eliminate linear interference, and sub-azimuth high precision multi-wave attenuation in Radon domain; (2) the near surface correction, curved surface static correction, surface Q absorption compensation, not diagenetic medium absorption compensation; (3) Prestack spectral expansion of low frequency protection, low frequency protection based on wavelet transform; (4) the residual static correction and velocity analysis, improve the high frequency residual static correction accuracy, meanwhile, velocity field for stack and poststack migration trial is established, and the PSTM processing and multiple iterations of velocity analysis are processed with bin size 10 m &#215; 10 m is processed for all the prestack data after the prestack CMP gather processing, and established accurate PSTM velocity field, and all sub-geometries data are conducted PSTM processing using the above velocity field.</p><p>The pre-stack CMP gathers are sorted according to the various trial geometries in <xref ref-type="table" rid="table1">Table 1</xref>, and the prestack time migration is processed and analyzed.</p></sec><sec id="s2_4"><title>2.4. Degradation Processing Effect Analysis</title><sec id="s2_4_1"><title>2.4.1. Comparison of Different Bin Size Effect with the Same Fold</title><p>The data were compared with fold 256 (geometries S0 and S1 of <xref ref-type="table" rid="table1">Table 1</xref>). From the PSTM profiles (<xref ref-type="fig" rid="fig3">Figure 3</xref>), compared with the profile with large bin size and low trace density, the profile with small bin size and high trace density has</p><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> PSTM profiles of bin size 10 m &#215; 10 m (a) and 20 m &#215; 20 m (b) with fold 256.</title></caption><fig id ="fig3_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x4.png"/></fig></fig-group><p>obviously higher SNR and resolution, and has better continuity for interlayer weak reflection and steep dip layer.</p><p>From the PSTM time slices at 2.5 s (<xref ref-type="fig" rid="fig4">Figure 4</xref>), compared with the time slice with large bin size and low trace density, the time slice with small bin size and high trace density has obviously higher SNR in the plane, stronger fault recognition ability, better seismic event continuity (as indicated by ellipses in the <xref ref-type="fig" rid="fig2">Figure 2</xref>.4-2), better lateral resolution (as indicated by arrows in the figure).</p><p>From the PSTM amplitude slices along the layer Es<sub>32</sub> (<xref ref-type="fig" rid="fig5">Figure 5</xref>), compared with the amplitude slices along the layer with large bin size and low trace density, the amplitude slices along the layer with small bin size and high trace density has obviously higher SNR in the plane, better seismic event continuity (as indicated by ellipses in the <xref ref-type="fig" rid="fig5">Figure 5</xref>), better lateral resolution of the geological body.</p><p>From the PSTM coherent slices (<xref ref-type="fig" rid="fig6">Figure 6</xref>), compared with the coherent slice with large bin size and low trace density, the coherent slice with small bin size and high trace density has obviously higher SNR, clearer fault characterization, more reasonable for fracture distribution.</p><p>In addition, the PSTM profiles of 10 m &#215; 10 m, 20 &#215; 10 m and 30 &#215; 10 m are compared and analyzed. Based on the analysis of the well data, for bin size 10 m &#215; 10 m and 20 m &#215; 10 m, the continuity of the carbonatite thin layer of Du3 and gaosheng oil group is equivalent, and the lateral resolution is also equivalent (as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>). For bin size 30 m &#215; 10 m, the continuity and the lateral resolution of the thin layer decreased obviously. For bin size 10 m &#215; 10 m when the area is 10 m &#215; 10 m, the super-geological phenomena at the bottom of buried hills is clear. When the dome is 20 m &#215; 10 m and the area is 30 m &#215; 10 m, the interpretation of the buried hilltop will appear. This further shows that the smaller bin size can improve the thin layer resolution and breakpoint recognition accuracy.</p><p>According to the above analysis, considering the influence of cost, the bin size</p><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> PSTM time slices of bin size 10 m &#215; 10 m (a) and 20 m &#215; 20 m (b) with fold 256.</title></caption><fig id ="fig4_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x5.png"/></fig></fig-group><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> PSTM amplitude slices along the layer (Es32) of bin size 10 m &#215; 10 m (a) and 20 m &#215; 20 m (b) with fold 256.</title></caption><fig id ="fig5_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x6.png"/></fig></fig-group><fig-group id="fig6"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> PSTM coherent slices of bin size 10 m &#215; 10 m (a) and 20 m &#215; 20 m (b) with fold 256.</title></caption><fig id ="fig6_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x7.png"/></fig></fig-group><p>10 m &#215; 10 m surface can be weakened appropriately in meeting the requirements of geological task. In the Liaohe Depression, the bin size of the second three-dimensional acquisition is 25 m &#215; 25 m, and it is suggested to optimize the</p><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Localgraph of the prestack time migration section of different bin sizes (a) 10 m &#215; 10 m; (b) 20 m &#215; 10 m; (c) 30 m &#215; 10 m</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x8.png"/></fig><p>bin 10 m &#215; 10 m to 12.5 m &#215; 12.5 m from the new and old data fusion.</p></sec><sec id="s2_4_2"><title>2.4.2. Comparison of Different Fold Effect with the Same Bin Size</title><p>The comparison and analysis of prestack time migration is carried out using three dataset of different folds and trace densities with the same bin size 10 m &#215; 10 m (geometries S0, S2 and S3 in <xref ref-type="table" rid="table1">Table 1</xref>).</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> shows prestack time migration profiles of fold 256, 128, and 64. When the fold reduced from 256 to 128, the signal-to-noise ratio of the target layer Es4 decreased a little, the continuity of seismic event became worse and the lateral resolution decreased a little. When the fold is further reduced to 64, the signal-to-noise ratio of the target layer Es4 decreased obviously, the seismic event became disordered, and the lateral resolution is significantly reduced obviously.</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref> shows time slices of fold 256, 128, and 64. When the fold reduced from 256 to 128, the signal-to-noise ratio of the target layer Es4 decreased a little, the continuity of seismic event became worse and the lateral resolution decreased a little. When the fold is further reduced to 64, the signal-to-noise ratio of the target layer Es4 decreased obviously, the seismic event became disordered, and the lateral resolution is significantly obviously.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows amplitude slices (Es32) of fold 256, 128, and 64. When the fold reduced from 256 to 128, the signal-to-noise ratio of the target layer Es4 decreased a little, the fault recognition ability reduced a little. When the fold is further reduced to 64, the signal-to-noise ratio of the target layer Es4 decreased obviously, Clearness of the fault was reduced obviously.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1 shows PSTM coherent slices of fold 256, 128, and 64. When the fold reduced from 256 to 128, the signal-to-noise ratio of the target layer Es4 decreased a little, the fault recognition ability reduced a little. When the fold is further reduced to 64, the signal-to-noise ratio of the target layer Es4 decreased obviously, clearness of the fault was reduced obviously.</p><p>According to the above analysis, in order to meet the requirements of geological task, considering the influence of cost, the fold can be reduced on the basis of 256, but higher than 128.</p></sec><sec id="s2_4_3"><title>2.4.3. Comparison of Big Bin Size and High Fold with Small Bin Size and Low Fold</title><p>Under the condition of the trace density as 640,000/km<sup>2</sup>, the comparison and</p><fig-group id="fig8"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> PSTM profiles of different folds 256 (a), 128 (b), 64 (c) with the same bin size 10 m &#215; 10 m.</title></caption><fig id ="fig8_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x9.png"/></fig></fig-group><fig-group id="fig9"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> PSTM time slices of different folds 256 (a), 128 (b), 64 (c) with the same bin size 10 m &#215; 10 m.</title></caption><fig id ="fig9_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x10.png"/></fig></fig-group><fig-group id="fig10"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> PSTM amplitude slices along the layer (Es32) of different folds 256 (a), 128 (b), 64 (c) with the same bin size 10 m &#215; 10 m.</title></caption><fig id ="fig10_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x11.png"/></fig></fig-group><fig-group id="fig11"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> PSTM coherent slices of different folds 256 (a), 128 (b), 64 (c) with the same bin size 10 m &#215; 10 m.</title></caption><fig id ="fig11_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x12.png"/></fig></fig-group><p>analysis of prestack time migration is carried out using two datasets (geometries S1 and S3 in <xref ref-type="table" rid="table2">Table 2</xref>.4-1) of big bin size (20 m &#215; 20 m) &amp; high fold (256) and small bin size (10 m &#215; 10 m) &amp; low fold (64).</p><p>From the PSTM profiles (<xref ref-type="fig" rid="fig1">Figure 1</xref>2), compared with the profile with small bin size (10 m &#215; 10 m) &amp; low fold (64), the profile with big bin size (20 m &#215; 20 m) &amp; high fold (256) has higher signal to noise ratio than the profile with small bin size (10 m &#215; 10 m) &amp; low fold (64).</p><p>From the PSTM time slices at 2.5 s (<xref ref-type="fig" rid="fig1">Figure 1</xref>3), compared with the profile with small bin size (10 m &#215; 10 m) &amp; low fold (64), the profile with big bin size (20 m &#215; 20 m) &amp; high fold (256) has higher signal to noise ratio, better seismic event continuity (as indicated by ellipses in <xref ref-type="fig" rid="fig1">Figure 1</xref>3).</p><p>From the PSTM amplitude slices of Es32 (<xref ref-type="fig" rid="fig1">Figure 1</xref>4), compared with the profile with small bin size (10 m &#215; 10 m) &amp; low fold (64), the profile with big bin size (20 m &#215; 20 m) &amp; high fold (256) has higher signal to noise ratio, better seismic event continuity (as indicated by ellipses in <xref ref-type="fig" rid="fig1">Figure 1</xref>4).</p><p>From the PSTM coherent slices (<xref ref-type="fig" rid="fig1">Figure 1</xref>5), compared with the profile with small bin size (10 m &#215; 10 m) &amp; low fold (64), the profile with big bin size (20 m</p><fig-group id="fig12"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>2</label><caption><title> PSTM profiles of different bin size and fold bin size 20 m &#215; 20 m &amp; fold 256 (a), bin size 10 m &#215; 10 m &amp; fold 64 (b).</title></caption><fig id ="fig12_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x13.png"/></fig></fig-group><fig-group id="fig13"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>3</label><caption><title> PSTM time slices of different bin size and fold, bin size 20 m &#215; 20 m &amp; fold 256 (a), bin size 10 m &#215; 10 m &amp; fold 64 (b).</title></caption><fig id ="fig13_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x14.png"/></fig></fig-group><fig-group id="fig14"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>4</label><caption><title> PSTM amplitude slices along the layer (Es32) of different bin size and fold, bin size 20 m &#215; 20 m &amp; fold 256 (a), bin size 10 m &#215; 10 m &amp; fold 64 (b).</title></caption><fig id ="fig14_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x15.png"/></fig></fig-group><fig-group id="fig15"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>5</label><caption><title> PSTM coherent slices along the layer (Es32) of different bin size and fold bin size 20 m &#215; 20 m &amp; fold 256 (a), bin size 10 m &#215; 10 m &amp; fold 64 (b).</title></caption><fig id ="fig15_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x16.png"/></fig></fig-group><p>&#215; 20 m) &amp; high fold (256) has higher signal to noise ratio, better seismic event continuity (as indicated by ellipses in the <xref ref-type="fig" rid="fig1">Figure 1</xref>5).</p></sec><sec id="s2_4_4"><title>2.4.4. Comparison of Different Aspect Ratio with the Same Bin Size</title><p>The comparison and analysis of prestack time migration is carried out using three datasets of different aspect ratio 0.91, 0.80 and 0.51 with the same bin size 10 m &#215; 10 m (geometries S0, S4 and S5 in <xref ref-type="table" rid="table1">Table 1</xref>).</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>6 shows prestack time migration profiles of different aspect ratio 0.91, 0.80, and 0.51. With the increase of the aspect ratio, the breakpoint imaging is gradually improved, and the deep signal-to-noise ratio is also improved.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>7 shows time slices of aspect ratio 0.91, 0.80, and 0.51. The difference is small between the data of wide azimuth (aspect ratio 0.91, 0.80). Compared with the profile with narrower azimuth (aspect ratio 0.51), the data of wide azimuth (aspect ratio 0.91, 0.80) has higher signal to noise ratio and better seismic event continuity.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>8 shows amplitude slices (Es32) of aspect ratio 0.91, 0.80, and 0.51. The difference is small between the data of wide azimuth (aspect ratio 0.91, 0.80). Compared with the profile with narrower azimuth (aspect ratio 0.51), the data of wide azimuth (aspect ratio 0.91, 0.80) has higher signal to noise ratio and better seismic event continuity.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>9 shows PSTM coherent slices of different aspect ratio 0.91, 0.80, and0.51. The difference is small between the data of wide azimuth (aspect ratio 0.91, 0.80). Compared with the profile with narrower azimuth(aspect ratio 0.51), the data of wide azimuth (aspect ratio 0.91, 0.80) has better thin layer continuity, better fault continuity and better fault identification ability.</p></sec><sec id="s2_4_5"><title>2.4.5. Comparison of Different Shot Line Interval with the Same Bin Size and Fold</title><p>The comparison and analysis of prestack time migration is carried out using two</p><fig-group id="fig16"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>6</label><caption><title> PSTM profiles of different aspect ratio with the same bin size 10 m &#215; 10 m 0.91 (a), 0.80 (b), 0.51 (c).</title></caption><fig id ="fig16_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x17.png"/></fig></fig-group><fig-group id="fig17"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>7</label><caption><title> PSTM time slices of different aspect ratio with the same bin size 10 m &#215; 10 m 0.91 (a), 0.80 (b), 0.51 (c).</title></caption><fig id ="fig17_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x18.png"/></fig></fig-group><fig-group id="fig18"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>8</label><caption><title> PSTM amplitude slices along the layer (Es32) of different aspect ratio with the same bin size 10 m &#215; 10 m 0.91 (a), 0.80 (b), 0.51 (c).</title></caption><fig id ="fig18_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x19.png"/></fig></fig-group><fig-group id="fig19"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>9</label><caption><title> PSTM coherent slices of different aspect ratio with the same bin size 10 m &#215; 10 m 0.91 (a), 0.80 (b), 0.51 (c).</title></caption><fig id ="fig19_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x20.png"/></fig></fig-group><p>datasets (geometries S2 and S5 in <xref ref-type="table" rid="table1">Table 1</xref>), which has the same bin size 10 m &#215; 10 m, the same fold and the same trace density, but different shot line interval (440 m and 220 m).</p><p>For the survey with the same fold, the fold of the middle-shallow Layer is high for the survey with small line interval, the distribution of middle-shallow layer fold is uniform, and the footprint is weak. For the survey with large line interval, the distribution of middle-shallow layer fold is uneven and it will produce a significant acquisition footprint.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>0 shows prestack time migration profiles of shot interval 440 m, 220 m. For the survey with shot line interval 440 m, the noise of middle-shallow Layer is serious, and the acquisition footprint is obvious, the signal to noise ratio of deep layer is low, but for the survey with shot line interval 220 m, the signal to noise ratio of middle-shallow Layer is high, the acquisition footprint is weak, the signal to noise ratio of deep layer is high, the seismic event continuity is good.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>1 shows time slices of different line intervals. For the survey with small shot line interval (220 m), the signal to noise ratio is high, and the amplitude attribute is uniform, but for the survey with big shot line interval (440 m), seismic event continuity is interrupted significantly, the amplitude attribute uniformity is poor, and this is caused by acquisition footprints.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>2 shows amplitude slices of Es32 of different line intervals. For the survey with small shot line interval (220 m), the signal to noise ratio is high, the seismic event is clear and the seismic event continuity is good, but for the survey with big shot line interval (440 m), data acquisition footprints is serious.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>3 shows PSTM coherent slices at 2.5 s. Compared with the coherent slice with big line interval (shot line interval 440 m), the coherent slice with small line interval (shot line interval 220 m) has better noise suppression effect, more clear fault imaging.</p><p>In order to ensure the amplitude attribute fidelity and uniformity and improve the accuracy of reservoir prediction, we should adopt a smaller line interval to minimize the impact of the acquisition footprint on the data.</p></sec></sec></sec><sec id="s3"><title>3. Effect of the Seismic Data</title><p>Based on the above research, the optimization direction of high precision seismic in Liaohe depression is determined, which provide guidance for the establish of the following 3D seismic acquisition geometry. HX 3D acquisition adopted the suggested geometry parameters proposed in this paper, bin size: 12.5 m &#215; 12.5 m, fold: 396 (shooting with vibroseis, and the fold is twice over shooting with charge), receiver line interval/shot line interval: 200 m/200 m, aspect ratio: 0.82, and achieved good results. The signal-to-noise ratio of new</p><fig-group id="fig20"><label><xref ref-type="fig" rid="fig2">Figure 2</xref>0</label><caption><title> PSTM profiles of different shot line interval 440 m (a), 220 (b).</title></caption><fig id ="fig20_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x21.png"/></fig></fig-group><fig-group id="fig21"><label><xref ref-type="fig" rid="fig2">Figure 2</xref>1</label><caption><title> PSTM time slices of different shot line interval 440 m (a), 220 m (b).</title></caption><fig id ="fig21_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x22.png"/></fig></fig-group><fig id="fig22"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref>2</label><caption><title> PSTM amplitude slices along the layer (Es32) of different shot line interval 440 m (a), 220 m (b)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x23.png"/></fig><p>acquisition high precision 3D data is high, and the resolution is improved obviously [<xref ref-type="bibr" rid="scirp.73691-ref13">13</xref>] . For the new acquisition data, it is blank weak reflection on west side of Jia Zhang Si fault, and is disorder reflection on the east side of Jia Zhang Si fault, further confirmed on both sides of the fault have the characteristics of lithology mutation. For the old data, it is disorder reflection on both sides of Jia Zhang Si fault, which does not reflect lithology mutations (<xref ref-type="fig" rid="fig2">Figure 2</xref>4).</p></sec><sec id="s4"><title>4. Conclusions</title><p>Theoretically, the higher the fold, the smaller the bin size, and the higher the trace density, the better the seismic data quality. In fact, optimizing the cost-effective acquisition parameters under the premise of meet the requirements of geological tasks, is an important objective of geometry parameter optimization. Comprehensive consideration of thin-layer identification, fault imaging and reduce the acquisition cost, the geometry parameters are recommended as follows: bin size: 12.5 m &#215; 12.5 m, fold: 200, aspect ratio: more than 0.80. The above recommended parameters are validated in the other 3D</p><fig-group id="fig23"><label><xref ref-type="fig" rid="fig2">Figure 2</xref>3</label><caption><title> PSTM coherent slices of different shot line interval 440 m (a), 220 m (b).</title></caption><fig id ="fig23_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x24.png"/></fig></fig-group><fig id="fig24"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref>4</label><caption><title> Conventional 3D acquisition profile (a) and high precision 3D prestack migration spectrum (b)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-2801362x25.png"/></fig><p>acquisition area in Liaohe depression.</p><p>Based on the high precision 3D acquisition data, degradation processing and comparative analysis are the effective means to optimize the geometry parameters. And the best observation system parameters can be found to meet the geological task requirements and reduce the cost. The method is also suitable for other regions.</p></sec><sec id="s5"><title>Cite this paper</title><p>Wang, H.B., Zhang, W., Liu, B., Zou, Q.W., Xiao, G.H., Tang, M.B., Zhu, H. and Li, W.J. (2017) Optimization on the High Precision 3D Seismic Acquisition Parameters in Leijia Tight Sandstone Oil Area. International Journal of Geosciences, 8, 30-45. http://dx.doi.org/10.4236/ijg.2017.81003</p></sec></body><back><ref-list><title>References</title><ref id="scirp.73691-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Wang, H.B., Liu, B., Zou, Q.W., et al. (2016) 3D Seismic Acquisition Technology and Effect in HX Volcanic Area in Liaohe Depression. International Journal of Geosciences, 7, 944-955.</mixed-citation></ref><ref id="scirp.73691-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, Y.H., Chen, X.Q. and Wang, N.J. (2015) Optimization of High Precision 3D Seismic Acquisition Geometry in Tazhong Oilfield. China Petroleum Institute Geophysical Technology Conference, 105-108.</mixed-citation></ref><ref id="scirp.73691-ref3"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Xu</surname><given-names> Y. </given-names></name>,<etal>et al</etal>. (<year>2014</year>)<article-title>Optimization on the High-Precision 3D Seismic Acquisition Parameters in Tahe Oilfield</article-title><source> Geophysical Prospecting for Petroleum</source><volume> 1</volume>,<fpage> 68</fpage>-<lpage>77</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.73691-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Wei, Z., Yin, C., Tian, J.D., et al. (2007) Degradation Processing Experiment Based on 3-D Geometry Parameters. CPP, 46, 69-73.</mixed-citation></ref><ref id="scirp.73691-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ling, Y., et al. (2005) Analysis and Appreciation of Wide/Narrow Azimuth Exploration Cases(II). Oil Geophysical Prospecting, 40, 423-427.</mixed-citation></ref><ref id="scirp.73691-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Ling, Y., Gao, J., Sun, D.S., et al. (2005) Analysis and Appreciation of Wide/Narrow Azimuth Exploration Cases (1). OGP, 40, 305-308, 317.</mixed-citation></ref><ref id="scirp.73691-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Pecholcs, P.I., Zhang, Y. and Lafon, S.K. (2010) Distance Separated Custom Slip-Sweep—A New High-Productivity Method. Eage Workshop on Developments in Land Seismic Acquisition for Exploration, 41-53.</mixed-citation></ref><ref id="scirp.73691-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Wloszczowski, D., Gou, Y. and Faraj, A. (1998) 3D Acquisition Parameters: A Cost-Saving Study. Expanded Abstracts of 68th Annual International SEG Mtg, 70-73.</mixed-citation></ref><ref id="scirp.73691-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Wang, H.B., Liu, Y.K. and Zou, Q.W. (2016) Single Point High Density 3D Seismic Acquisition Technical Studies in Leijia Tight Sandstone Oil Reservoir of Liaohe Depression. Progress in Geophysics, 31, 0782-0787. (In Chinese)</mixed-citation></ref><ref id="scirp.73691-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, B.Q., Jin, S.T., Zeng, T.J., et al. (2015) Application of “Two High and One Wide” Seismic Exploration Technology in the Eastern Margin of the Caspian Sea Basin. China Petroleum Institute Geophysical Technology Conference, 82-85.</mixed-citation></ref><ref id="scirp.73691-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Li, P.M., Kang, N.C., Zou, X.F., et al. (2013) The Key Technology of High Precision Seismic Exploration in “Two-Wide One-High” High Precision Seismic Exploration. Memoir of 2013 China Geophysical Nineteenth Sub Assembly Hall. Petroleum Industry Press, Beijing.</mixed-citation></ref><ref id="scirp.73691-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kong, D.Z., Yu, M.J., Liu, X.W., et al. (2016) Application and Effects of the “Two-Wide One-High” Seismic Acquisition Technique in Complex Piedmont. Xinjiang ShiYou TianRan Qi, 12, 33-38.</mixed-citation></ref><ref id="scirp.73691-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Wang, X.J., Yu, B.L., Zhao, X.H., et al. (2015) Development and Application of “2W1H” Technique in Oil and Gas Exploration. China Petroleum Exploration, 20, 41-53.</mixed-citation></ref></ref-list></back></article>