Shale gas plays a main role in the recent hydrocarbons exploration industry. North part of Saudi Arabia is one of the most promising unconventional gas reserves in shale gas reservoir geologically and known by Qusaiba Formation. The Qalibah Group is located in this part of the Arabian shield and consists of two Formations: the (upper) Silurian Sharawra Formation and the (lower) Qusaiba shale Formation. Silurian Sharawra consists mainly of sandstone and siltstone with beds of shale. The important of Silurian Sharawra is that it is laying over one the biggest rock source of hydrocarbon in the region Qusaiba Formation. The Formation is presented from outcrop exposures in subsurface in this part of Saudi Arabia . The second Formation of Surilian Qalibah Paleozoic is Qusaiba Formation; it is classified as an organic-rick shale. The Formation mostly consists of shale and claystone embedded with small amount of siltstone and sandstone . They reported the Qusaiba Formation occurrence in two sub-basins in Central and Eastern of Saudi Arabia. Moreover, the Qusaiba Formation is eroded recording to most recent reports in the area between the Central Arch Block anticline oil and gas North Gawar field, and North horst graben oil and gas of Saudi Arabia . The expectation of oil and gas existences in Qusaiba could reach 37 billion barrels and 808 trillion cubic feet, respectively .
To explore the unconventional shale gas many studies and field works were conducted in several locations in Saudi Aribia. One of the interesting studies on unconventional shale gas was done in North West of Saudi Arabia to investigate pore pressure prediction for the hot shale of the Qusaiba formation and Sarah Formation using multi well logs data in the area supported with 3D seismic data . They found the thickness of Qusaiba is diverse across the basin and that hot shale presents partially as seal formation for Sarah formation. The study was focused on the characterization of well logs and the pore pressure changes of Qusaiba hot shale without mentioned to the true depth and thickness of the formation. Jones and Stump  indicated that the thick part of Qusaiba formation is located South of Ar Riyad in the Southern basin. In contrast, at the North of Saudi Arabia, the Qusaiba Formation becomes thinner between 150 m to 300 m at the study area . At the base of Qusaiba shale formation, there is a presence of black “hot shales” that is rich in organic matter. It is considered as a proven source rock of the Paleozoic hydrocarbons ranges in thickness between 30 - 70 m  -.
Here, we aim to utilize the integration of high resolution reflection seismic and gravity data to image the sequence of subsurface layers and structures of near surface shale gas formations and verify of the reservoir characterization depth and thickness for unconventional gas fields exploration and production in the region.
In Tayma area in Saudi North part, determining the rock source Qusaiba Formation and detecting the top depth of the Formation will be most beneficial to verify the depth of hot shale gas and mitigate the uncertainties of source rock Surillian Qusaiba and Sarah Formation Figure 1. Furthermore, to study the differences in depths and dipping of the reservoir in the North part of Arabian Peninsula comparing with amount of study what have done on the basin in Central
Figure 1. Map showing the location of study area at North part of Saudi Arabia. The black line in enlarge map indicates the seismic survey location which is identical to the gravity survey as well. In zoomed map, Ga Bohan is showed below the survey line which has a very hard surface.
and South Arch host of oil and gas.
2. Method and Fieldwork
Using the seismic imaging was so critical because the Qusaiba shale in North West is very shallow comparing to its depths in Central and South part. In this study, 2D seismic profile was acquired using 96 channel and 408 shots. The intervals of sources and receivers were 10m with the minimum and maximum offsets 5 m and 950 m, respectively. We utilize one Mini IVI Vibrosies source to shot 408 linear upsweep with bandwidth between 20 - 200 Hz and stacking 3 numbers of sweeps. The field Parameters were selected based on testing shots that showed the best parameters for the nature geological phenomena and area accessibility. The layout of seismic survey in study area was spreading on the top of Sharawra shale layer. It was suggested to have the seismic survey on the surrefce of Qa Buhan, but the idea changes because of some field observations Figure 2. Figure 2(A) shows the direction of the survey line which is placed on the hill instead what suggested on the surface Figure 2(C). The reason to change the location of survey line refers to hardness of Qa Buhan surface which makes spreading spike geophones into ground more difficult. In addition, Gravity data was acquired along the seismic line with 4100 m length and interval 100 m between stations to detect and map an unexplored lower Silrurian Qusaiba-Plaeozoic Formation. We mainly used gravity method to measure the distribution of density of the subsurface, mapping geologic model based on external constraints, and adding more understanding of the geological subsurface
Figure 2. Shows the environment of the field observations. In photo (A), it can be clear the high of the hill which represents an outcrop of the Sharawra formation. The arrow in photo (B)refers to the direction of the survey location. Photo (D) shows the elevation of the hill and the arrow refers to other side of the area which is called Qa Buhan. Photo (C) shows the hard surface of the study area which makes the spreading of the seismic survey layout is an impossible especially with spike geophones.
structures. In the gravity interval, recording locations were identical with the seismic survey and exactly over the locations of the seismic source shots.
3. Results and Discussions
Commonly, Shale is heterogeneous and brittleness undersurface. The surface Sharawra formation is mainly composed of sandstone with many beds siltstone and shale. The strength of our source (Mini IVI Vibrosies source) is prop for near surface application but according to the lithology of the subsurface, we recognized absence of some high frequencies underneath and appear a wide of the ringing waves. These ringing and discontinuous horizons may refer to the lithology diversity between litstone and sandstone beds and brittleness shales layers inside Sharawra Formation in shallow part, which are marked by the red square in seismic section Figure 3. Sparitosh and D, Thomas  found that a strong seismic expression and multiples ringing in the continuously of the P-wave gathers were existed in the shallow reservoir of sandstone surrounding by shale beds causes. That may be caused by the differences in the acoustic impedance between sandstone and shale that gives a weak contrast. Moreover, the stress and orientation of shale beds with lower and upper litstone and sandstone may cause more attenuation of high frequencies signals.
In Figure 3, seismic section also shows strong horizon events at depth 230 and at depths between 370 - 390 m. The horizons denote the tops of the Qusaiba “hot
Figure 3. Shows the prestack time migration of the area with gather intervals 10 m and maximum offsets 950 m, for 409 shots. The red square indicates some shallow shale beds and demonstrates the evidences of the Sharawra shale formation on the surface. Also in the section the detection of the top of Qusaiba formations (TOQF) is imaged on the 230 m with thickness between 70 - 90 m along the profile, and followed by the top of the Sarah formation (TOSF) at 300 to 320. Moreover, the interpretation of the section shows some shallow faults, which are marked by blue lines.
shale” (TOQF) and Sarah formation (TOSF) with some interactions of shallow geological structures and faults appearing in the middle of the section. According to high frequencies signals attenuation caused by sandstone surrounding Sharawra shale beds with lower and upper huge amplitude energy were loose and caused absent of the deep formation beneath the Sarah formation.
Figure 4 shows the diversity of both sandstone and shale interbeds along the outcrop in the area which is continued throughout subsurface Sharawra formation and that gives an evidence of the appearance of the near surface ringing waves in seismic section.
Gravity data was acquired along the seismic line with 4100 m length and interval 100 m between stations to detect and map an unexplored lower Silrurian Qusaiba-Plaeozoic Formation. Figure 5(A) shows the processed gravity data utilizing GM-SYS Profile modelling platform. The accurate processing and modeling results usually help to achieve more accurate and detailed subsurface novel and decrease the uncertainties in the interpretation modeling. The geological modelling observed from the gravity data shows brilliant mapping of TOQF and TOSF with roughly same detections of the depth and dipping in seismic imaged result. The both modeling gives the same depth of TOQF along the profile at depth roughly 230 m with change in range of 40 - 60 m. Note, the seismic result, here, was cropped in both sides to remove the affection of the migration imaging.
Figure 4. (A) Shows the interbeds shale and sandstone beds of Sharawra formation and (B) Shows the outcrop of the formation.
Figure 5. Shows the gravity processing profile and interpretation modeling result of the gravity survey of the area. The gravity survey is identical with the seismic along and recorded on the locations of the seismic source shots. (A) Shows a result of gridding 2D profile with some high anomalies in the end of the survey area. (B) GM-SYS Profile Modelling was utilized to map an accurate geological model by picking observed measurements of diverse subsurface gravity recording (D).
Therefore, the gravity geological model seems longer the seismic data Figure 6.
Mapping the subsurface of organic-rich shale gas of the Lower Silurian Qusaiba formation and Sarah formation are performed using the integration of seismic and gravity data. The investigation of the unconventional Qusaiba hot shale reservoir in North of Saudi Arabia is necessary to increase the possibilities of exploring the shale gas field in the area and to decrease the uncertainties of subsurface depth and thickness of the rock source. Imaging of the migrated 2D seismic data and modeling of 2D profile gravity data are performed in order to distinguish the subsurface formations and detect their depths and thicknesses. The present integration modeling and imaging results are showing identical
Figure 6. Shows results of the area of the geological gravity model (A) and the seismic image (B). The comparison between models shows a robust detection for both formations TOQF and TOSF. The depth scale and the length of the two models are different.
detections of three formations: Sharawra (on surface) followed by Qusaiba and Sarah formations. The depth of the shale gas Qusaiba formation is diverse and the continuously is interpret with some shallow geological faults. The result of migrated seismic prove the laying of Sharawra formation on the surface which is ascertained as well by the field surface scenes and evidences of the outcrop. The lithology of the Sharawra formation exhibition with a series of sandstone interfering with thin interbed shale along the Formation cause a wide of multiple ringing in shallow seismic waves. The work result adds a respected contribution in shale gas unconventional reservoir exploration and development in North part of the Saudi Arabia, but it needs more data for verification and covered field works to discover the area significantly. Moreover, it is highly important to increase the income profits for Saudi Arabia in shale gas.
The author would like to acknowledge King Abdulaziz City for Science and Technology for supporting this work. Many thanks shall extend the geophysical team of KACST for conducting some activities related to this study.
 Guo, R. and Weng, Y. (2017) Analysis of the Positive Effect from the Typhoon Saomai to the Hydrothermal Environment of Shanghai. Journal of Geoscience and Environment Protection, 5, 221-234. https://doi.org/10.4236/gep.2017.58018
 Ding, D. and Li, Y. (2011) A Study on Typhoon-Induced Rainfalls over Beijing: Statistics and Case Analysis. Journal of Meteorological Research, 25, 742-753.
 Mojgan, G., Mehdi, M. and Reza, B. (2017) The Trend of Changes in Surface Wind in the Indian Ocean, in the Period from 1981 to 2015, Using Reanalysis Data, NCEP/NCAR. Open Journal of Marine Science, 7, 445-457.
 Pattanaik, D.R., Mukhopadhyay, B. and Kumar, A. (2012) Monthly Forecast of Indian Southwest Monsoon Rainfall Based on NCEP’s Coupled Forecast System. Atmospheric and Climate Sciences, 2, 479-491. https://doi.org/10.4236/acs.2012.24042
 Zhang, X., Zhang, L., Fu, J. and Zhang, L. (2017) Analysis of Characteristics of the Forecast Jump in the NCEP Ensemble Forecast Products. Atmospheric and Climate Sciences, 7, 151-159. https://doi.org/10.4236/acs.2017.71011
 Xiong, X., Deng, W., Hu, J., et al. (2017) Design and Implementation of Real-Time Monitoring and Warning System for Regional Disastrous Weather Based on CIMISS. Meteorological Science & Technology, 45, 453-459.
 Berman, F., Casanova, H., Chien, A., et al. (2005) New Grid Scheduling and Rescheduling Methods in the GrADS Project. International Journal of Parallel Programming, 33, 209-229. https://doi.org/10.1007/s10766-005-3584-4
 Duan, B., Zhang, W., Liu, H., et al. (2017) The Spatial and Temporal Distributions of Warm Sector Rainfall and Frontal Rainfall for the Torrential Rain Event in Beijing on 21 July 2012. Torrential Rain & Disasters, 36, 108-117.