Recent Submissions

  • Imaging near-surface heterogeneities by natural migration of backscattered surface waves: Field data test

    Liu, Zhaolun; AlTheyab, Abdullah; Hanafy, Sherif M.; Schuster, Gerard T. (Society of Exploration Geophysicists, 2017-03-06)
    We have developed a methodology for detecting the presence of near-surface heterogeneities by naturally migrating backscattered surface waves in controlled-source data. The near-surface heterogeneities must be located within a depth of approximately one-third the dominant wavelength λ of the strong surface-wave arrivals. This natural migration method does not require knowledge of the near-surface phase-velocity distribution because it uses the recorded data to approximate the Green’s functions for migration. Prior to migration, the backscattered data are separated from the original records, and the band-passed filtered data are migrated to give an estimate of the migration image at a depth of approximately one-third λ. Each band-passed data set gives a migration image at a different depth. Results with synthetic data and field data recorded over known faults validate the effectiveness of this method. Migrating the surface waves in recorded 2D and 3D data sets accurately reveals the locations of known faults. The limitation of this method is that it requires a dense array of receivers with a geophone interval less than approximately one-half λ.
  • Opportunities and pitfalls in surface-wave interpretation

    Schuster, Gerard T.; Li, Jing; Lu, Kai; Metwally, Ahmed Mohsen Hassan; AlTheyab, Abdullah; Hanafy, Sherif (Society of Exploration Geophysicists, 2017-01-21)
    Many explorationists think of surface waves as the most damaging noise in land seismic data. Thus, much effort is spent in designing geophone arrays and filtering methods that attenuate these noisy events. It is now becoming apparent that surface waves can be a valuable ally in characterizing the near-surface geology. This review aims to find out how the interpreter can exploit some of the many opportunities available in surface waves recorded in land seismic data. For example, the dispersion curves associated with surface waves can be inverted to give the S-wave velocity tomogram, the common-offset gathers can reveal the presence of near-surface faults or velocity anomalies, and back-scattered surface waves can be migrated to detect the location of near-surface faults. However, the main limitation of surface waves is that they are typically sensitive to S-wave velocity variations no deeper than approximately half to one-third the dominant wavelength. For many exploration surveys, this limits the depth of investigation to be no deeper than approximately 0.5-1.0 km.
  • Imaging near-surface heterogeneities by natural migration of surface waves

    Liu, Zhaolun; AlTheyab, Abdullah; Hanafy, Sherif M.; Schuster, Gerard T. (Society of Exploration Geophysicists, 2016-09-06)
    We demonstrate that near-surface heterogeneities can be imaged by natural migration of backscattered surface waves in common shot gathers. No velocity model is required because the data are migrated onto surface points with the virtual Green's functions computed from the shot gathers. Migrating shot gathers recorded by 2D and 3D land surveys validates the effectiveness of detecting nearsurface heterogeneities by natural migration. The implication is that more accurate hazard maps can be created by migrating surface waves in land surveys.
  • Skeletonized inversion of surface wave: Active source versus controlled noise comparison

    Li, Jing; Hanafy, Sherif (Society of Exploration Geophysicists, 2016-07-14)
    We have developed a skeletonized inversion method that inverts the S-wave velocity distribution from surface-wave dispersion curves. Instead of attempting to fit every wiggle in the surface waves with predicted data, it only inverts the picked dispersion curve, thereby mitigating the problem of getting stuck in a local minimum. We have applied this method to a synthetic model and seismic field data from Qademah fault, located at the western side of Saudi Arabia. For comparison, we have performed dispersion analysis for an active and controlled noise source seismic data that had some receivers in common with the passive array. The active and passive data show good agreement in the dispersive characteristics. Our results demonstrated that skeletonized inversion can obtain reliable 1D and 2D S-wave velocity models for our geologic setting. A limitation is that we need to build layered initial model to calculate the Jacobian matrix, which is time consuming.
  • Extracting 220 Hz information from 55 Hz field data by near-field superresolution imaging

    Dutta, Gaurav; AlTheyab, Abdullah; Tarhini, Ahmad; Hanafy, Sherif; Schuster, Gerard T. (Oxford University Press (OUP), 2016-05-31)
    Field experiments are used to unequivocally demonstrate seismic superresolution imaging of subwavelength objects in the near-field region of the source. The field test is for a conventional hammer source striking a metal plate near subwavelength scatterers and the seismic data are recorded by vertical-component geophones in the far-field locations of the sources. Time-reversal mirrors (TRMs) are then used to refocus the scattered energy with subwavelength resolution to the position of the original source. A spatial resolution of lambda/10, where lambda is the dominant wavelength associated with the data, is seen in the field tests that exceeds the Abbe resolution limit of lambda/2.
  • Geotechnical Parameters from Seismic Measurements: Two Field Examples from Egypt and Saudi Arabia

    Khalil, Mohamed H.; Hanafy, Sherif M. (GeoScienceWorld, 2016-03-18)
    © 2016 EEGS. Geotechnical parameters were used to determine subsurface rock quality for construction purposes. We summarize the mathematical relationships used to calculate the geotechnical parameters from P- and S-wave velocities and density values. These relationships are applied to two field examples; the first is a regional seismic study in Egypt and the second is a 2-D seismic profile recorded in Saudi Arabia. Results from both field examples are used to determine the subsurface rock quality and locate zones that should be avoided during construction. We suggest combining all geotechnical parameters into one map using a normalized-weighted relation, which helps to locate the zones with high versus low rock quality for engineering purposes.
  • Controlled Noise Seismology

    Hanafy, Sherif M.; AlTheyab, Abdullah; Schuster, Gerard T. (Society of Exploration Geophysicists, 2015-08-19)
    We use controlled noise seismology (CNS) to generate surface waves, where we continuously record seismic data while generating artificial noise along the profile line. To generate the CNS data we drove a vehicle around the geophone line and continuously recorded the generated noise. The recorded data set is then correlated over different time windows and the correlograms are stacked together to generate the surface waves. The virtual shot gathers reveal surface waves with moveout velocities that closely approximate those from active source shot gathers.
  • Mapping the Qademah Fault with Traveltime, Surface-wave, and Resistivity Tomograms

    Hanafy, Sherif M. (Society of Exploration Geophysicists, 2015-08-19)
    Traveltime, surface-wave, and resistivity tomograms are used to track the buried Qademah fault located near King Abdullah Economic City (KAEC), Saudi Arabia. The fault location is confirmed by the 1) resistivity tomogram obtained from an electrical resistivity experiment, 2) the refraction traveltime tomogram, 3) the reflection image computed from 2D seismic data set recorded at the northern part of the fault, and 4) the surface-wave tomogram.
  • Qademah Fault Seismic Data Set - Northern Part

    Hanafy, Sherif M.; Lu, Kai; Hota, Mrinal Kanti; Guo, Bowen; Tarhini, Ahmad (2015-01)
    Objective: Is the Qademah fault that was detected in 2010 the main fault? We collected a long 2D profile, 526 m, where the fault that was detected in 2010 is at around 300 m. Layout: We collected 264 CSGs, each has 264 receivers. The shot and receiver interval is 2 m. We also collected an extra 48 CSGs with offset = 528 to 622 m with shot interval = 2 m. The receivers are the same as the main survey.
  • Iterative supervirtual refraction interferometry

    Al-Hagan, Ola; Hanafy, Sherif M.; Schuster, Gerard T. (Society of Exploration Geophysicists, 2014-05-02)
    In refraction tomography, the low signal-to-noise ratio (S/N) can be a major obstacle in picking the first-break arrivals at the far-offset receivers. To increase the S/N, we evaluated iterative supervirtual refraction interferometry (ISVI), which is an extension of the supervirtual refraction interferometry method. In this method, supervirtual traces are computed and then iteratively reused to generate supervirtual traces with a higher S/N. Our empirical results with both synthetic and field data revealed that ISVI can significantly boost up the S/N of far-offset traces. The drawback is that using refraction events from more than one refractor can introduce unacceptable artifacts into the final traveltime versus offset curve. This problem can be avoided by careful windowing of refraction events.
  • Qademah Fault Artificial Ambient Noise Test

    Hanafy, Sherif M.; AlTheyab, Abdullah (2014)
    This data set was collected on 7 Dec. 2014 by Sherif and Abdullah. The receiver layout is the same as that of the passive data test at the same location, which is described as follow: 288 receivers are used and arranged as follow - 12 lines, cross-line offset = 10 m - 24 receiver in each line, inline offset = 5 m - Additional 24 receivers are placed at line # 6, where the receiver interval is decreased to 2.5 m. Data Recording: We start recording at 10:10 am and stop recording at 11:25 am. Each record has total of 20 s, with time interval of 0.004 ms and around 2 s overlap between each two successive files. Source: We used a piece of wood attached to a pick-up truck to create the noise; we drove around the array of receivers in a rectangle-shape route during the recording time.
  • Qademah Fault 3D Survey

    Hanafy, Sherif M.; Lu, Kai; Hota, Mrinal Kanti; Guo, Bowen; Tarhini, Ahmad (2014)
    Objective: Collect 3D seismic data at Qademah Fault location to 1. 3D traveltime tomography 2. 3D surface wave migration 3. 3D phase velocity 4. Possible reflection processing Acquisition Date: 26 – 28 September 2014 Acquisition Team: Sherif, Kai, Mrinal, Bowen, Ahmed Acquisition Layout: We used 288 receiver arranged in 12 parallel lines, each line has 24 receiver. Inline offset is 5 m and crossline offset is 10 m. One shot is fired at each receiver location. We use the 40 kgm weight drop as seismic source, with 8 to 15 stacks at each shot location.
  • Qademah Fault Passive Data

    Hanafy, Sherif M.; Lu, Kai; Hota, Mrinal Kanti; Guo, Bowen; Tarhini, Ahmad (2014)
    OBJECTIVE: In this field trip we collect passive data to 1. Convert passive to surface waves 2. Locate Qademah fault using surface wave migration INTRODUCTION: In this field trip we collected passive data for several days. This data will be used to find the surface waves using interferometry and then compared to active-source seismic data collected at the same location. A total of 288 receivers are used. A 3D layout with 5 m inline intervals and 10 m cross line intervals is used, where we used 12 lines with 24 receivers at each line. You will need to download the file (rec_times.mat), it contains important information about 1. Field record no 2. Record day 3. Record month 4. Record hour 5. Record minute 6. Record second 7. Record length P.S. 1. All files are converted from original format (SEG-2) to matlab format P.S. 2. Overlaps between records (10 to 1.5 sec.) are already removed from these files