Removal of Intra-Array Statics in Seismic Arrays Due to Variable Topography and Positioning Errors

Featured Application: The proposed corrections can be used in two different seismic applications. The ﬁrst is to calculate accurate static corrections, which are applied to seismic reﬂection data. The second is high-resolution engineering applications. Abstract: A receiver array is an arrangement of geophones used to enhance the signal-to-noise ratio (S/N) of seismic data. However, deviations from ideal array conditions can lead to the non-optimal performance of the array. This study investigates, quantitatively, the array performance in the presence of topographic variations and positioning errors using 2D seismic data acquired in eastern Saudi Arabia. A receiver array was laid over a sand dune with variable topography underlain by a ﬂat sabkha that has a very shallow water table. The topographic variations and position errors were calculated from Differential Global Positioning System (D-GPS) measurements of source and receiver positions and elevations. The errors in receiver positions, measured relative to the ideal receiver spacing, gave a mean and standard deviation of about 0.35% and 1%, respectively. On the other hand, elevation errors (topographic variations) from a horizontal datum gave a mean and standard deviation of about 25% and 13%, respectively. The ideal array response was found by removing both elevation and position errors. The ﬁrst-arrival energy of the array was calculated after removal of elevation and position errors separately and compared to the ideal-array energy. Comparison showed a 64% enhancement in the ﬁrst-arrival energy after correcting for elevation errors alone and almost no enhancement after correcting for position errors alone. The proposed approach can be used to calculate accurate static corrections for seismic reﬂection processing and to generate high-resolution subsurface images for engineering applications.


Introduction
A seismic array is a group of seismic receivers and/or sources planted in a systematic arrangement.In general, seismic arrays are used to suppress noise in the field.Many of the existing studies on seismic arrays assume ideal array conditions with perfect element (geophone) weights, perfect implementation conditions, and simple near-surface layers.However, these assumptions are rarely satisfied in the field, hence limiting the applicability of the results of these studies to field data.In particular, it is difficult to have similar ground coupling of all array elements, leading to variable element weights.Furthermore, element spacing is practically estimated by survey crew, leading to often slightly variable element positions.Local variations in the near-surface layer, such as topography and heterogeneity, can also lead to more non-ideal array conditions.All these conditions may decrease the effectiveness of the array's multi-channel processing and its waveform [1].Originally, seismic arrays were developed to monitor nuclear explosions [2].The studies of nuclear explosion identification have been employed to distinguish between a natural earthquake and nuclear explosion [3,4].Ref. [5] shows that seismic arrays improve the velocity models of the Earth's interior using large-number and dense spatial sampling of the seismic wavefield at several locations on the Earth's surface.More recently, seismic arrays have been implemented in high-resolution imaging [6,7], shallow geology applications [8][9][10], monitoring induced seismicity [11], seismic imaging using ambient seismic noise [12], and seismic signal processing using distributed acoustic sensing (DAS) [13].Ref. [14] proposed an efficient algorithm for ambient seismic noise beamforming without the need for explicit cross-correlations, while [15] presented a case study using a 10-month ambient DAS data set acquired on a fiber array to investigate long-term seasonal velocity variations.Seismic arrays are often integrated with other geophysical methods for a better understanding of the subsurface conditions [10,[16][17][18] In terms of reflection seismology, the main purpose of an array of geophones is to increase the signal-to-noise ratio (S/N) by enhancing the desired signals (reflections) and attenuating seismic noise during the acquisition of seismic data in the field.As the desired seismic signals (reflections) travel across the array with a much higher horizontal velocity, the signals produce a longer wavelength, as compared to that of the unwanted signals [19].The desired signals are usually those that have been reflected from deep geologic horizons approaching the array at almost normal incidence.The unwanted signals are incoherent noise or horizontally traveling ground roll.Thus, a geophone array can also be considered as a directional filter [20].
Seismic array response is the arithmetic sum of traces recorded by geophones (elements) in the array.Ref. [21] presented general concepts of arrays, where they used harmonic analysis to explain the array response disrupted by variations in velocity, dip, and elevation within the array length.Ref. [22] described a method for computing the expected waveform from a closely spaced array of shots or geophones when the individual signal is a pulse and studied the effects of various array weighting functions on the performance of seismic arrays.Ref. [23] presented a study on the effects of errors in position and weighting on the array response, where they used plane harmonic waves with different weights.Ref. [24] investigated effects of errors in gain and position on the array response.Ref. [25] explained the wavelet response due to errors in weights and positions of the elements of the geophone arrays with respect to several variables.These variables included the error amount, number of array elements, and weighting function.They found that the degradation of the wavelet response increased nonlinearly with the error's value, number of elements, and weighting function complexity.A similar case was investigated by [26], wherein they evaluated the effects on the array response caused by topographic variations (i.e., deviations from the datum) within receiver arrays.They found that these variations had more negative effects on the array response than position errors because of their relatively high error amount in most environments.In addition, topographic variations directly affect the vertically traveling signals as opposed to position errors, which mainly affect noise.Finally, topographic variations leave very little human control to reduce their amount, whereas position errors can be reduced using more accurate instruments and measurement procedures.
Ref. [27] studied the performance of an array in the presence of lateral near-surface heterogeneities.They used a model with two adjacent layers to study the array performance as a function of the number of affected elements and the velocity ratio between the two adjacent layers.They found that the array performance decreased with the number of affected elements and amount of deviation of the velocity ratio from unity.Ref. [28] studied the change in the wavelet response of seismic receiver arrays in the presence of a dipping near-surface layer.They used a model with different dip and incidence angles as well as various types of signal and noise incident on the array.Their results show that incidence and dip angles have opposite effects on the array response.Refs.[29,30] studied the array response in the presence of laterally varying thickness of the weathering layer.They generated several models and calculated their corresponding array responses in the ideal and perturbed cases.They found that horizontally traveling waves were less affected than vertically traveling waves in all studied models and errors.Ref. [31] studied the seismic array response in the presence of variations in elements' elevations, positions, and weights using a zero-phase Ricker wavelet.They found that topography had the largest effects on the array response and that different error types did not combine linearly.
Almost all previous studies have addressed the effects of non-ideal array conditions using analytical (or numerical) modeling approaches without testing their approaches on field data.The problem with field data studies is the difficulty in measuring the exact position and elevation of all receivers.In this paper, we summarize, modify, and update equations required for correcting non-ideal array conditions; then, we study the array response in the presence of topographic variations (i.e., elevation errors) and positioning errors within the array using field data from a shallow seismic survey acquired in eastern Saudi Arabia.We quantify the elevation and position errors (relative to a datum and nominal geophone spacing, respectively) at every geophone location using their Differential Global Positioning System (D-GPS) coordinates.Furthermore, we quantify enhancements in the array response using the stacked energy of first arrivals in each error case normalized by its corresponding error-free case.
This paper starts with the calculation of the velocities of the direct and refraction events in the presence of elevation and position errors.Next, the elevation and position errors are calculated using the survey geometry information.This is followed by calculation of the array response in terms of its trace energy after elevation error removal, position error removal, and removal of both errors.We finally compare them with the trace response of the ideal case (i.e., after removing both elevation and position errors).

Direct Event
The direct P-wave can be defined as the P-wave that travels through the ground surface directly from the source to the receiver.To calculate the velocity (V 1 ) of the incidence (uppermost) medium, Equation (1) is used [32]: where T d is the arrival travel time of the direct wave to a receiver at offset X.The arrival travel times of the direct wave are picked (estimated) from the traces (records), while offsets X are given by the survey geometry information.

Refraction Event
The head wave is generated when the velocity of the refraction medium (V 2 ) is larger than that of the incidence medium and the incidence angle is equal to the critical angle θ c .Here, the transmitted wave travels just below and parallel to the interface.To find the velocity of the refraction medium, the following time-distance equation is used: where T h is the arrival traveltime of the head wave at offset X and H is the thickness of the uppermost layer.Equation (2) describes a straight line relating T h and X with an intercept of 2H cos θ c /V 1 and a slope of 1/V 2 .The refraction velocity in this case can be calculated by finding the inverse of the slope of the best-fit line to the T h − X plot of the head wave.
In general, when a wave encounters an interface between two media, Snell's law is used to relate incident and transmitted angles [32]: where θ 1 and θ 2 are the ray incident and transmitted angles in the upper (V 1 ) and lower (V 2 ) media, respectively.The critical angle θ c takes place when θ 2 = 90 o , which can be calculated as: In the case of a dipping refractor, the critical angle θ c can be calculated from the apparent refractor velocities in the up-dip (V 2u ) and down-dip (V 2d ) directions using Equation (5) [32]: The apparent refractor velocities can be determined by picking arrival travel times from shots in the up-dip and down-dip sides.Once θ c is calculated, the refractor medium velocity V 2 can be determined as [32]:

Position and Elevation Error 2.3.1. Position Error
Ideally, there is a desirably fixed nominal spacing between successive receivers, which we call the ideal position (spacing) ∆X.However, in practice, receiver spacing might differ from the ideal spacing due to field conditions and human errors (Figure 1).The actual receiver spacing at the nth receiver (∆X cn ) can be calculated from D-GPS measurements acquired in the field as follows: where x n and x n−1 are the D-GPS Easting and y n and y n−1 are the Northing values of the nth receiver and its neighbor closer to the source.To calculate the position's error at the nth receiver, Equation ( 8) is used: which is the difference between the ideal and calculated positions of the nth receiver.
where and are the ray incident and transmitted angles in the upper ( ) and lowe ( ) media, respectively.The critical angle takes place when = 90 , which can be calculated as: In the case of a dipping refractor, the critical angle can be calculated from the apparent refractor velocities in the up-dip ( ) and down-dip ( ) directions using Equa tion (5) [32]: The apparent refractor velocities can be determined by picking arrival travel time from shots in the up-dip and down-dip sides.Once is calculated, the refractor medium velocity can be determined as [32]:

Position Error
Ideally, there is a desirably fixed nominal spacing between successive receivers which we call the ideal position (spacing) ∆ .However, in practice, receiver spacing might differ from the ideal spacing due to field conditions and human errors (Figure 1) The actual receiver spacing at the nth receiver (∆ ) can be calculated from D-GPS meas urements acquired in the field as follows: where and are the D-GPS Easting and and are the Northing values o the nth receiver and its neighbor closer to the source.To calculate the position's error a the nth receiver, Equation ( 8) is used: which is the difference between the ideal and calculated positions of the nth receiver.

Elevation Error
Recording seismic data at sites with rough terrain or big elevation differences (Figure 2) is challenging.The big difference in elevation affects the data results, particularly re garding accurate travel time picking.To overcome this problem, we first calculate the el evation error, then remove it from the recorded data.Equation ( 9) is used to calculate the elevation error of the nth receiver:

Elevation Error
Recording seismic data at sites with rough terrain or big elevation differences (Figure 2) is challenging.The big difference in elevation affects the data results, particularly regarding accurate travel time picking.To overcome this problem, we first calculate the elevation error, then remove it from the recorded data.Equation ( 9) is used to calculate the elevation error of the nth receiver: where Z n and Z D are the elevations of the nth receiver and the datum (reference plane), respectively, from mean sea level, and are available from the D-GPS measurements.
Appl.Sci.2022, 12, 12810 5 of 1 where and are the elevations of the nth receiver and the datum (reference plane) respectively, from mean sea level, and are available from the D-GPS measurements.

Perturbed Case
The perturbed case is defined in this study as the raw data acquired in the field with out applying any time corrections.The travel time of the nth receiver for the perturbed case ( ) is the picked first arrival traveltime with no elevation and/or position corrections

Position's Error Removal
The travel time difference due to position error (∆ ) is the travel time of the nth receiver calculated using Equation (10): To calculate the first arrival travel time of the nth receiver for the position remova error case ( ), ∆ must be removed from the perturbed case traveltime ( ) as fol lows:

Elevation Error Removal
The travel time difference due to the elevation error (∆ ) of the nth receiver i calculated using Equation ( 12): Similar to the travel time calculation of the position error removal case, the first arri val travel time of the elevation's error removal case ( ) can be calculated by removing ∆ from the perturbed case traveltime using Equation (13): The ideal case can be defined as the case after applying both position and elevation corrections on the raw data acquired in the field.The first arrival travel time of the idea case for the nth receiver ( ) can be calculated from the perturbed case travel time using Equation ( 14):  The perturbed case is defined in this study as the raw data acquired in the field without applying any time corrections.The travel time of the nth receiver for the perturbed case (T n ) is the picked first arrival traveltime with no elevation and/or position corrections.

Position's Error Removal
The travel time difference due to position error (∆Te xn ) is the travel time of the nth receiver calculated using Equation (10): To calculate the first arrival travel time of the nth receiver for the position removal error case (T xn ), ∆Te xn must be removed from the perturbed case traveltime (T n ) as follows:

Elevation Error Removal
The travel time difference due to the elevation error (∆Te zn ) of the nth receiver is calculated using Equation ( 12): Similar to the travel time calculation of the position error removal case, the first arrival travel time of the elevation's error removal case (T zn ) can be calculated by removing ∆Te zn from the perturbed case traveltime T n using Equation (13):

Ideal Case
The ideal case can be defined as the case after applying both position and elevation corrections on the raw data acquired in the field.The first arrival travel time of the ideal case for the nth receiver (T cn ) can be calculated from the perturbed case travel time using Equation ( 14): Appl.Sci.2022, 12, 12810 6 of 16

Wavelet Response and Trace Energy
To calculate the trace energy in the first arrival event, the following steps are followed: 1.
Window the first wavelet of the first arrival head-wave event while maintaining the timing information (T n ) of the first arrival on each trace.

2.
Shift each nth windowed trace such that its first sample is at t = 0.This produces the ideal trace wavelet response (G cn (t)).Sum G cn (t) across n to produce the ideal array response (G c (t)) and calculate its corresponding array trace energy (E c ) as follows: where M is the number of samples in G c (t).
The decibel (dB)-scaled normalized array energy with position errors only E Nx (dB) is calculated as follows: where E Nx is the normalized array trace energy with position errors only defined as: Starting with G cn (t) again, shift each G cn (t) by its corresponding ∆Te zn to produce the trace wavelet response with elevation errors only (G zn (t)).Sum G zn (t) across n to produce the array response with elevation errors only (Gz(t)) and calculate its corresponding array trace energy (E z ) as follows: The dB-scaled normalized array energy with elevation errors only E NZ (dB) is calcu- lated as follows: where E Nz is the normalized array trace energy with elevation errors only defined as: The dB-scaled normalized array trace energy for the perturbed case E NxZ (dB) is calculated as follows: where E Nxz is the normalized array trace energy for the perturbed case defined as:

Geological Setting
The field data were recorded over a barchan dune near Dammam city (26 • 20 57.2 N 49 • 49 51.4 E), Saudi Arabia (Figure 3a,b).The thickness of the dune ranges from 0 to 7 m and it overlays a flat Sabkha layer.One borehole was drilled up to 10 m deep at the midslope of the dune, and was used as ground truth and for result confirmation.A complete core was extracted from the top 3.5 m, then a core sample of only 0.7 m was extracted at each 1.5 m depth range.The well penetrated the entire dune, as well as the top 3 m of the underlying sabkha.
The dB-scaled normalized array trace energy for the perturbed case ( ) is calculated as follows: where is the normalized array trace energy for the perturbed case defined as:

Geological Setting
The field data were recorded over a barchan dune near Dammam city (26°20′57.2″N 49°49′51.4″E), Saudi Arabia (Figure 3a,b).The thickness of the dune ranges from 0 to 7 m and it overlays a flat Sabkha layer.One borehole was drilled up to 10 m deep at the midslope of the dune, and was used as ground truth and for result confirmation.A complete core was extracted from the top 3.5 m, then a core sample of only 0.7 m was extracted at each 1.5 m depth range.The well penetrated the entire dune, as well as the top 3 m of the underlying sabkha.Based on the samples extracted from the borehole, the dominant mineral in the dune is quartz.We also found small amounts of feldspar and traces of calcite.The sabkha layer contains quartz, feldspar, calcite, and significant amounts of halite and gypsum.The top 2.5 m of the sabkha was composed of fine-grained to medium-grained sand, while the bottom 0.5 m contained muddy sand with a noticeable fraction of fine grains.Sediments Based on the samples extracted from the borehole, the dominant mineral in the dune is quartz.We also found small amounts of feldspar and traces of calcite.The sabkha layer contains quartz, feldspar, calcite, and significant amounts of halite and gypsum.The top 2.5 m of the sabkha was composed of fine-grained to medium-grained sand, while the bottom 0.5 m contained muddy sand with a noticeable fraction of fine grains.Sediments in the top 7.6 m of the drilled interval showed partial saturation, whereas the bottom 2.4 m was fully water saturated.

Data Acquisition
We demonstrate our methodology using a shallow seismic survey acquired in eastern Saudi Arabia (Figure 3a) at the approximate GPS coordinates (26.351575N, 49.829184 E). Figure 3b shows the survey location where the geological setting was a sand dune over Sabkha.The acquired data consisted of a 2D seismic line with 24 vertical receivers having a nominal spacing of 4 m.Five shots were produced by a 10 kg sledgehammer located at offsets of −4, 0, 46, 92, and 96 m, where the first receiver was placed at offset 0 m.Positions and elevations of all receiver and shot points were measured using D-GPS.

First Arrival Picking
For illustration purposes, we applied the proposed method to the 12 receivers between shots 1 and 3 (Figure 3b,c).Far-offset receivers were not used because the signal became too weak to pick confidently.Figure 4 shows the first and third shot gathers of the real data set.Arrival travel times were picked manually for both direct and head waves.For the direct wave, the first arrival travel time was picked only on the two receivers closest to the source since direct waves were observed only on these traces.The first picked arrival travel times of the direct wave from all shots are shown in Table 1, whereas the picks of the head wave from the down-dip (shot 1) and up-dip (shot 3) are shown in Table 2.
m was fully water saturated.

Data Acquisition
We demonstrate our methodology using a shallow seismic survey acquired in e ern Saudi Arabia (Figure 3a) at the approximate GPS coordinates (26.351575 N, 49.829 E). Figure 3b shows the survey location where the geological setting was a sand dune o Sabkha.The acquired data consisted of a 2D seismic line with 24 vertical receivers hav a nominal spacing of 4 m.Five shots were produced by a 10 kg sledgehammer locate offsets of −4, 0, 46, 92, and 96 m, where the first receiver was placed at offset 0 m.Posit and elevations of all receiver and shot points were measured using D-GPS.

First Arrival Picking
For illustration purposes, we applied the proposed method to the 12 receivers tween shots 1 and 3 (Figure 3b,c).Far-offset receivers were not used because the sig became too weak to pick confidently.Figure 4 shows the first and third shot gathers of real data set.Arrival travel times were picked manually for both direct and head wa For the direct wave, the first arrival travel time was picked only on the two receivers c est to the source since direct waves were observed only on these traces.The first pic arrival travel times of the direct wave from all shots are shown in Table 1, whereas picks of the head wave from the down-dip (shot 1) and up-dip (shot 3) are shown in T 2.

Velocity Calculation
The average dune velocity is V 1 = 313 m/s, calculated from the direct-wave velocities in the five shot records (Table 1).Head-wave velocities V 2 are calculated from the slopes of the best-fit lines to the first arrival traveltimes of the up-dip and down-dip directions (Figure 5).The critical angle and true head-wave velocity are calculated by Equations ( 5) and ( 6) to be θ c = 9.2 • and V 2 = 1949 m/s, respectively (Table 2).

Error Calculation
To calculate position errors in the study area, Equation ( 8) is used, where the ideal spacing ∆X = 4 m is assumed and the calculated spacing ∆X cn is measured from D-GPS Easting and Northing coordinates for each receiver (Table 3).Results of the position errors show that the mean and standard deviation of the errors are 0.014 and 0.041 m (0.35% and 1.025% of ideal spacing), respectively, which are negligible.Moreover, Equation ( 9) is used to calculate the elevation errors where a horizontal datum is chosen at the first shot elevation (i.e., datum = 9.4 m above mean sea level).Elevation errors are shown in Table 4.In comparison to position errors, elevation errors have a mean and standard deviation of 2.36 and 1.22 m (25.1% and 13% relative to the datum elevation), respectively, which are relatively considerable.Therefore, it is expected that most effects on the seismic array response are attributed to elevation errors in this case.The first arrival travel times of the perturbed and ideal cases are shown in Figure 6a.

Error Calculation
To calculate position errors in the study area, Equation ( 8) is used, where the ideal spacing ∆ = 4 m is assumed and the calculated spacing ∆ is measured from D-GPS Easting and Northing coordinates for each receiver (Table 3).Results of the position errors show that the mean and standard deviation of the errors are 0.014 and 0.041 m (0.35% and 1.025% of ideal spacing), respectively, which are negligible.Moreover, Equation ( 9) is used to calculate the elevation errors where a horizontal datum is chosen at the first shot elevation (i.e., datum = 9.4 m above mean sea level).Elevation errors are shown in Table 4.In comparison to position errors, elevation errors have a mean and standard deviation of 2.36 and 1.22 m (25.1% and 13% relative to the datum elevation), respectively, which are relatively considerable.Therefore, it is expected that most effects on the seismic array response are attributed to elevation errors in this case.The first arrival travel times of the perturbed and ideal cases are shown in Figure 6a.

Responses Calculation
After calculating the position and elevation errors, we calculate the array responses in various cases using the procedure described in Section 2.4 above.The ideal trace wavelet responses (G cn (t)) are shown in Figure 7a and their corresponding array wavelet response (G c (t)) is shown in Figure 7b.The array trace energy for the ideal case (i.e., after removal of both position and elevation errors) is E c = 2237.The perturbed trace wavelet responses (G xzn (t)) are shown in Figure 7c and their corresponding array wavelet response (G xz (t)) is shown in Figure 7d.The array trace energy for the perturbed case is E xz = 797.The normalized array trace energy of the perturbed case is calculated using Equation (23) to be E Nxz = −9 dB.These results demonstrate the considerable amount of degradation in the array trace energy due to position and elevation errors in this survey.To understand the separate effects of position and elevation errors on the seismic array response, the first arrival travel times after separately removing each of these errors is calculated.Starting with elevation error removal, Equation ( 13) is used to calculate the first arrival travel time for this case.Comparison between the travel times of the ideal case and the elevation error removal case shows that elevation errors have large contributions in all picked head-wave first arrival travel times (Figure 6b). Figure 7e shows the trace wavelet responses after application of elevation errors to the ideal trace responses (G zn (t)) and Figure 7f shows the corresponding array response (G z (t)).Now, the array trace energy is calculated to be E z = 795, which is much smaller than the array trace energy of the ideal case (i.e., E c = 2237).The normalized array trace energy is E Nz = −9 dB, which is similar to that in the perturbed case, indicating that most degradation is attributed to elevation errors.
The same steps are carried out for the removal of only the position errors using Equation (11) and, hence, update the first arrival travel times.Figure 6c shows a comparison of the head-wave travel times between the ideal and position error removal cases.This comparison indicates that position errors have a very small effect on the first arrival travel times.The trace wavelet responses after position error removal and their corresponding array response are shown in Figure 7g,h, respectively.The array trace energy after the removal of only position errors is calculated to be E x = 2237, which is almost equal to that of the ideal case, leading to a normalized array trace energy E Nx = 0 dB.Similar results were also reported by previous studies [30,31], as discussed in the next section.

Comparison with Previous Studies
We compare our results to those expected by [31] by applying their method to calculate array trace responses involving position errors only, elevation errors only, and combined position and elevation errors using the parameters of our data.The head-wave arrival in our data has a dominant frequency of 65 Hz and is incident at 9.2 • on a 12-element array that has an ideal spacing of 4 m laid on a dune with an average velocity of 313 m/s.wavelet responses after application of elevation errors to the ideal trace responses ( and Figure 7f shows the corresponding array response ( ( )).Now, the array trace ergy is calculated to be = 795, which is much smaller than the array trace energ the ideal case (i.e., = 2237).The normalized array trace energy is = −9 , w is similar to that in the perturbed case, indicating that most degradation is attribute elevation errors.The same steps are carried out for the removal of only the position errors using E tion (11) and, hence, update the first arrival travel times.Figure 6c shows a compariso the head-wave travel times between the ideal and position error removal cases.This c parison indicates that position errors have a very small effect on the first arrival tr times.The trace wavelet responses after position error removal and their correspon array response are shown in Figure 7g,h, respectively.The array trace energy after Similar to [31], we first generated the ideal array response due to a 65 Hz zero phase Ricker wavelet incident at 9.2 • on a 12-element equally weighted array and calculated the corresponding trace energy E c = 1060.57.Then, we calculated array responses with various error types using the following procedure: 1.
Normalize the standard deviation of position errors (Tables 3 and 4) by the ideal geophone spacing.2.
Generate 12 random position errors from a normal distribution with zero mean and normalized standard deviation of position errors calculated in step 1.

3.
Generate the array response with position errors only.

4.
Calculate the corresponding trace energy.
We repeat these steps for the cases of elevation errors only and combined position and elevation errors.For the case of combined position and elevation errors, we assume normally distributed independent position and elevation errors and calculate a combined standard deviation of 0.3052 [33].Table 5 presents calculations using this procedure and compares them to their counterparts calculated using steps described in Sections 2.2 and 2.3 of this paper.Similar to the observed data, the array response due to position errors shows only a negligible degradation (0.009%) when calculated using the method of [31].On the other hand, observed degradations due to elevation errors only (64.5%) and combined position and elevation errors (64.4%) are smaller than those expected by the method of [31] (88.0% and 87.9%, respectively).Nevertheless, observed degradations due to elevation errors only and combined position and elevation errors are comparable to each other and much higher than the degradation due to position errors only, which is also observed by the method of [31].Differences in our results from those calculated using the method of [31] may be attributed to the presence of environmental noise in our data.In addition, imperfect geophone couplings within our array deviate it from the perfect equally weighted array used in [31].Finally, our hammer source generates a slightly variable minimum-phase wavelet in contrast to the exactly repeatable zero-phase Ricker wavelet used by [31].

Conclusions
The objective of this paper is to investigate the effects of both elevation and positioning errors on the seismic array response.One field data set was recorded using 24 receivers and five different shot gathers.Processing and interpreting the recorded data show two subsurface layers with velocities of 313 and 1949 m/s for the first and second layers, respectively, and a critical angle of θ c = 9.2 • .Knowing the exact source and receiver positions and elevations, the elevation errors and position errors are calculated.Our results show that elevation has higher errors relative to position, where the mean errors of elevation and position are 2.4 and 0.014 m, respectively.The reason for the higher errors in the elevation case is the big difference between the receiver elevations and the horizontal datum elevation, which reached a maximum of 4 m.The first arrival travel times and the array trace energy are calculated for all studied cases.The first case is the perturbed case involving the first arrival travel times of the head wave without applying any correction for position or elevation errors.The second case is the ideal case involving the first arrival travel times of the head wave after applying corrections for both position and elevation errors.Results show that there is a considerable difference between the two cases, where the array trace energy of the ideal case is enhanced by almost three times in comparison with the array trace energy of the perturbed case.Afterwards, the first arrival travel times and the array trace energy of position and elevation error removal cases are calculated separately and compared to the ideal case.Results show a considerable decrease in the array trace energy after elevation error removal, and almost no change after position error removal, compared to the array trace energy of the ideal case.These results assure us that almost all signal energy loss is due to the presence of elevation errors in this study area.Therefore, it is necessary to correct for this error to minimize energy degradation of the seismic array, which is very important for many steps of seismic data processing such as static corrections, velocity analysis, and imaging.
The two main applications of the proposed technique are; (1) processing deep seismic reflection data and (2) high-resolution imaging for engineering applications.In both applications, it is essential to obtain accurate source and receiver positioning.Although the elevation errors have larger effects on the final data, both position and elevation need to be accurately measured.
The main limitation of the proposed technique is the need for D-GPS to obtain accurate position and elevation of receiver points.This could be an obstacle in some cases; however, the more accurate the position/elevation readings, the more accurate the final results.

Figure 1 .
Figure 1.Illustration of error in receiver positions.Black and gray triangles indicate ideal and actua positions of the nth receiver, respectively.

Figure 1 .
Figure 1.Illustration of error in receiver positions.Black and gray triangles indicate ideal and actual positions of the nth receiver, respectively.

Figure 2 .
Figure 2. Illustration of errors in receiver elevations.Gray triangles indicate actual elevations o receivers, while black triangles indicate ideal elevations of receivers (after datum correction).

Figure 2 .
Figure 2. Illustration of errors in receiver elevations.Gray triangles indicate actual elevations of receivers, while black triangles indicate ideal elevations of receivers (after datum correction).

Figure 3 .
Figure 3. (a) Satellite view showing Saudi Arabia and the study area marked with red stars; (b) image showing the study area with the seismic line (red line); and (c) the seismic line including 12 receivers (blue triangles) and two sources (red stars) with their actual positions and elevations using D-GPS data.

Figure 3 .
Figure 3. (a) Satellite view showing Saudi Arabia and the study area marked with red stars; (b) image showing the study area with the seismic line (red line); and (c) the seismic line including 12 receivers (blue triangles) and two sources (red stars) with their actual positions and elevations using D-GPS data.

16 Figure 5 .
Figure 5. Raw first arrival travel times of shot gathers 1 (circles) and 3 (crosses) picked manually for the 12 traces used in this study.

Figure 5 .
Figure 5. Raw first arrival travel times of shot gathers 1 (circles) and 3 (crosses) picked manually for the 12 traces used in this study.

Figure 6 .
Figure 6.Comparison between head-wave first arrival travel times of the 10 studied traces: (a) perturbed case; (b) elevation error removal case; and (c) position error removal case.

Figure 6 .
Figure 6.Comparison between head-wave first arrival travel times of the 10 studied traces: (a) perturbed case; (b) elevation error removal case; and (c) position error removal case.

Figure 7 .
Figure 7. Trace wavelet responses (G) for: (a) ideal case; (c) perturbed case; (e) elevation erro moval case; and (g) position error removal case.The corresponding array wavelet response of t cases are: (b) ideal case; (d) perturbed case; (f) elevation error removal case; and (h) position removal case.

Figure 7 .
Figure 7. Trace wavelet responses (G) for: (a) ideal case; (c) perturbed case; (e) elevation error removal case; and (g) position error removal case.The corresponding array wavelet response of these cases are: (b) ideal case; (d) perturbed case; (f) elevation error removal case; and (h) position error removal case.
3. Shift each G cn (t) by its corresponding ∆Te xn to produce the trace wavelet response with position errors only (G xn (t)).Sum G xn (t) across n to produce the array response with position errors only (G x (t)) and calculate its corresponding array trace energy (E x ) as follows: . Starting with G cn (t) again, shift each G cn (t) by its corresponding (∆Te xn + ∆Te zn ) to produce the perturbed wavelet response of the trace with combined position and elevation errors (G xzn (t)).Sum G xzn (t) across n to produce the array response with both position and elevation errors (G xz (t)) and calculate its corresponding array trace energy (E xz ) as follows:

Table 1 .
First arrival travel times of direct wave and corresponding offsets and velocities from all shots of the survey.The two traces from shot gather 3 (marked by *) have been excluded from averaging because they were considered outliers.

Table 2 .
First arrival travel times of the head wave from shot gather 1 and 3 with the values of ground-surface dip angle (α), critical angle (θ c ) and refraction medium velocity (V 2 ).

Table 3 .
Measured position coordinates using D-GPS and position errors of the 12 studied traces.Position errors are measured with respect to the ideal spacing of 4 m giving a relative mean and standard deviation of about 0.35% and 1%, respectively.

Table 3 .
Measured position coordinates using D-GPS and position errors of the 12 studied traces.
Position errors are measured with respect to the ideal spacing of 4 m giving a relative mean and standard deviation of about 0.35% and 1%, respectively.Trace X (m) Y (m) Calculated Distance (m) e xn (m)

Table 4 .
Elevation values and errors of the 12 studied traces.Errors are measured with respect to the datum elevation at 9.41 m giving a relative mean and standard deviation of about 25% and 13%, respectively.

Table 5 .
[31]arison of results of this study to those of[31].