Annals n.6/2003 ok 23/04 1325 ANNALS OF GEOPHYSICS, VOL. 46, N. 6, December 2003 Key words geomagnetic survey – sub-surface structures – aeromagnetic data – seismic lines 1. Introduction The studied area lies in the southern part of the Sinai Peninsula, Gulf of Suez and northern part of the Eastern Desert (fig. 1) between lati- tudes 27°50′ and 30°05′ north and longitudes 32°48′ and 34°23′ east covering an area of about 24 000 km2. It is a Tertiary cratonic rift south of Egypt’s Sinai Peninsula between North-Eastern Africa and the Arabian Peninsu- la. The rifting essentially ceased (Early-Middle Miocene, 18-14 Myr) when continental separa- tion became more oblique because of the dom- inant movement on the left-lateral transform fault that extends through the Gulf of Aqaba north-eastward to the Dead Sea (Patton et al., 1994; Lindquist, 1998). The tectonics of this area has been construct- ed using different geophysical tools (e.g., aero- magnetic, land magnetic gravity, and seismic). These tectonics were analyzed from deep depths to show depths along the studied area. 2. Geology of the area The Red Sea, Gulf of Suez basins consists of a high faulted block more than 100 000 km2 in size that contains block-faulted trapping con- figurations in the Gulf of Suez. This unit in- Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data Taha Rabeh National Research Institute of Astronomy and Geophysics (NRIAG), Geophysical Department, Helwan, Cairo, Egypt Abstract This work deals with the evaluation of the structural set-up of the southern part of Sinai Peninsula, Gulf of Suez and western part of Gulf of Suez from magnetic, gravity, and seismic data. The utilised techniques including the Least Squares separation method suggest NW, NE, and E-W trends. The trend analysis shows north 35°-45° west, north 15°-25° east and E-W which may be related to the Gulf of Suez and Red Sea stresses.The Euler deconvolution illustrates that the area is highly affected by these trends. Depths range from 1 km to more than 3 km below sea level and its magnetic susceptibility ranges between 1 to 3 SI units. The 2.5D magnetic model- ling and analytical signal techniques confirm the depths to the magnetic sources deduced by the Euler method, whereas the depth to the basement rocks ranges between 0 km to about 3 km indicating that it is subjected to strong tectonic activities. In addition, two seismic sections (EG-31 and MP-70), compiled by the Egyptian Gen- eral Petroleum Cooperation (EGPC), were interpreted together with a geologic cross section. The studied area may be divided into several major blocks along the Gulf of Suez area. It can be concluded that the sedimentary was affected by basement tectonics as revealed by the two seismic sections. Mailing address: Dr. Taha Rabeh, National Research Institute of Astronomy and Geophysics (NRIAG), Geophy- sical Department, Helwan, Cairo, Egypt; e-mail: tahara- beh@yahoo.com 1326 Taha Rabeh cludes coastal areas of the Red Sea and the Southern Gulf of Suez. The geology of the area ranges from Pre- cambrian basement rocks to the Quaternary de- posits. According to the Egyptian Geological Survey map (EGSMA, 1993) (fig. 2), the Qua- ternary deposits cover the basement rocks along the Gulf of Suez and Gulf of Aqaba. The basement rocks (Precambrian) occupy the southern part of the studied area, along the Gulf of Suez and Gulf of Aqaba. On the other hand, the Palaeozoic rocks are present in the middle part while the Mesozoic rocks covered the eastern and western parts of the area. Furthermore, a geologic cross section was constructed (fig. 4), using the most updated borehole data. This section was performed crossing Gulf of Suez taking a northeast-south- west direction. It indicates that the depth to the basement rocks varies from 0 km to more than 3 km. The basement rocks consist of majorFig. 1. Location map of the studied area. Fig. 2. Geological map of the studied area (EGSMA, 1993). 1327 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data Fig. 3. Location of the drilled wells, the two seismic lines and profiles of land magnetic survey (L-L′). Fig. 4. Geologic cross section along profile G-G′ crossing the Gulf of Suez. 1328 Taha Rabeh Fig. 5. RTP aeromagnetic anomaly map of the studied area (EGPC, 1990). Fig. 6. Bouguer anomaly map of the studied area (after IPRG, ca. 1982). Fig. 7. Least square residual RTP aeromagnetic map (2nd order). Fig. 8. Least square residual Bouguer anomaly map (2nd order). 1329 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data the input magnetic and gravity data. This method was applied using the software Surfer v. 7.00 (Surfer, 1997). The correlation coefficients be- tween successive residual maps were computed to determine the optimum order of the regional surface to be used. Results were r21 = 0.854, r21 = 0.96 and r34 = 0.546. (4.1) According to Davis (1973), a correlation fac- tor higher than 0.7 indicates a good similarity between maps, the regional field in this area can be represented by a second-order surface, and the polynomial used has the form (x, y)2 = − 7.947 − 0.4778y + 0.065x + + 0.08144y − 0.03575xy − 0.002528x. (4.2) The resulting least-squared residual maps of second-order (figs. 7 and 8) show that most of the anomalies are the same trends trending north west-south east, north east-south west and east- west directions. Also, it can be concluded that most of these anomalies are related to different blocks with different polarities and structures. Furthermore these maps represent the best fitted maps to the original RTP aeromagnetic and Bouguer anomaly maps. 5. Trend analysis technique A two lineaments map (fig. 9) was con- structed from application of the trend analysis technique (Grant et al., 1965; Parasnis, 1997) for the RTP aeromagnetic and Bouguer anom- alies maps. This map was subjected to statisti- cal calculating for the lengths percentages of the faults measured clockwise from the north. The deduced lineation of different directions are grouped every 10° around the north for their length percentage (L%). A statistical procedure was used to illustrate the predominant fault trends affecting the studied area. The results of azimuth distribution of both surface and subsurface elements are presented in the form of a rose diagram (fig. 10). The re- sults show that the most predominant direction is north 35°-45° west, north 45°-65° east, east- west and Aqaba trends. It is clear that the main faulted blocks along the Gulf of Suez. The depth to the carrier reservoir bed (Belayim Fm.) varies from about 2 km to about 3 km. 3. Data analysis and interpretations The land magnetic survey was performed using two magnetic Protons of Geometric type. One instrument is used as a base station for di- urnal corrections for each profile, and the sec- ond is used to measure the observed magnetic data. The land survey was performed in the form of profiles along the Gulf of Suez area, (fig. 3). The distance between the stations ranges from 200 to 300 m according to the changes in the recorded geomagnetic data. The aeromagnetic map was prepared by the EGPC (1990) with a flight elevation of 1 km. The aeromagnetic data was digitized into a computer using a suitable Surfer Program v. 8.00 (Surfer, 2002). Both land magnetic data and the aeromag- netic data were reduced to the North Magnetic Pole (Hansen and Pawlowski, 1989), to obtain the RTP magnetic anomaly map (fig. 5). Both the RTP aeromagnetic (fig. 5) and Bou- guer Anomaly maps (fig. 6) indicate that most anomalies are aligned to NW, NE, and EW which may be related to the Gulf of Suez, and Gulf of Aqa- ba. The negative anomalies in the central part may be due to lithology variation of the basement. The Bouguer Anomaly map (IPRG, ca. 1982) was digitized and contoured at 1 mGal (fig. 6). This map conforms to the observed RTP anom- alies (fig. 5). It reveals a gravity lows in the central and southern parts. These alternative negative and positive anomalies along the Gulf of Suez may be due to the faulted blocks or presence of different basins with different thickness of sedimentary se- quence in the area. This leads to the conclusion that the shallow parts are extending along the two Gulfs and southern part of Sinai where the base- ment rocks are outcropping (fig. 2). 4. Least squared residual technique The Least Square method of Henderson (1966) was applied to the RTP magnetic Bouguer maps using first, second, and third order fitted to 1330 Taha Rabeh Fig. 10. Rose diagrams illustrate the distribution of the lengths percentage (L%) of faults against azimuth. Fig. 9. Tectonic lines as deduced from the RTP aeromagnetic and Bouguer anomaly maps. tectonic trend is north 35°-45° west prevailing for the large scale in both the Bouguer and RTP aeromagnetic maps. This may be due to a very strong and active force resulting from the open- ing of the Gulf of Suez. The second predomi- nant tectonic force is north 45°-65° east that is probably due to the Syrian Arc system. The least predominant is the north 15°-25° east due to tectonics related to the Gulf of Aqaba. 6. Euler deconvolution method The two-dimensional Euler technique is a method that uses the horizontal and vertical de- rivatives of the total magnetic field to estimate the location and depth to magnetic sources (Thomp- son, 1982; Reid et al., 1990). This method has been applied to the magnetic profiles along the area (figs. 11 and 12). Applying Euler deconvolution to the vertical gradient of the total magnetic field data provides improved source resolution (Stavrev, 1997; Hsu, 2002) . x z f x x x z f z z z f BSI n n n n n n z 0 0 2 2 2 2 2 2 2 2 2 2 - + - = = +∆ J L K K e ^ e ^ e N P O O o h o h o (6.1) where n is the order (which is not necessarily an integer (Cooper and Cowan, 2003)) of the gra- dient used. It may similarly be applied to the horizontal gradient . x z f x x z x f z z x f BSI n n n n n n x 0 0 2 2 2 2 2 2 2 2 2 2 - + - = = +∆ J L K K e ^ e ^ e N P O O o h o h o (6.2) In general the magnetic field falls off as r – N with r = (x2 + z2)1/2. The parameter x is the location of the magnetic source along the magnetic profile and z is the depth to the top of the magnetic source (Thompson, 1982). The value of N, called the structural index, depends on the shape of the mag- netic source. For example, N = 3 corresponds to a three-dimensional body such as a sphere while a 1331 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data Fig. 11. Euler deconvolution along L1-L1′ of the land magnetic data. Fig. 12. Euler deconvolution along L3-L3′ of the land magnetic data. 1332 Taha Rabeh two-dimensional dike corresponds to N = 2. Val- ues of N between 1 and 0.5 correspond to two-di- mensional faults and contacts. The depth to the structures depends on structural index N and the size of the applied sliding window. The Euler deconvolution method has been applied to the aeromagnetic using 0 magnetic contact index: to determine the depths to the ir- regular magnetic bodies, 0.5 magnetic step in- dex in order to indicate the depth to the base- ment rocks and their structures and 1 dike index for delineating the depth to the anomalies caused by vertical thin dykes. Most the results derived from the Euler de- convolution method indicate that area is highly affected by the regional forces related to the Gulf of Suez, and Red Sea tectonics. It is indicated by the different position and depths of the sources of the subsurface magnetic bodies (figs. 11 and 12), which consist of dissected bodies along these profiles. Also, it can be stated that these magnet- ic sources are restricted by tectonic trends have NW-SE and NE-SW directions (fig. 13). These structures were confirmed by the analytical signal applied to the same profiles. The magnetic susceptibility of the rocks rang- es between 1 to 3 SI units whereas their depths range between 1 km to about 4 km. These wide ranges in the depth to the subsurface basement rocks are due to the highly faulted areas. 7. Application of the analytical signal and depth estimation methods The analytical signal was performed along the same profiles interpreted by the Euler method by Reid et al. (1990) using Nabighian’s technique (Nabighian, 1972, 1984). It is applied to confirm the results deduced from the Euler method (figs. 12 and 13). The spikes of the analytical signal profile are an indication for the subsurface struc- tures. The presence of one spike indication of subsurface fault structure, and the presence of two closed spikes are due to the presence of sub- surface dike, the high amplitude spike for the near edge of the dike and the small amplitude spike is due to the lower edge of the dike. Also, the depth to the magnetic sources was applied to the analytical signal using (Atchuta et al., 1981; Roest et al., 1992). The depth to the magnetic sources equals half of the ampli- tude of the analytical signal. This method had been applied to amplitudes along the profiles. The results confirm the parameters deduced from the Euler method. The depth to these bod- ies ranges between 0.5 km to more than 3 km. 8. Application of 2.5D modelling Computations of the magnetic effects for 2.5D models with complex geometry were car- ried out using the GM-SYS computer program produced by Northwest Geophysical Associates Inc. (GM-SYS, 1995), for the arbitrary two-di- mensional polygon. This type of modeling is a re- verse model that depends on the changes in the magnetic susceptibility, or the density contrasts of the sedimentary rocks and the depth parame- ters till complete matching between the calculat- ed and observed profiles (fig. 14). The resulting parameters represent the main target of this meth- od (e.g., depth, magnetic susceptibility, etc.). This 2.5D model confirms the depths to the magnetic sources deduced by the Euler method. It can be noticed from (fig. 14) that the depth to the Fig. 13. Applications of Euler deconvolution to the aeromagnetic map. 1333 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data Fig. 14. 2.5D modelling along L1-L1′ of the land magnetic profile. Fig. 15. The depth to the basement rocks along the Gulf of Suez area, derived from the wells drilled in the area. 1334 Taha Rabeh basement rocks ranges between 0 km to about 3 km while the magnetic susceptibility reaches about 0.0005 CGS unit. Also, it is obvious that this result is correlated with the results deduced by the Euler deconvolution and analytical signal method. 9. Depth to the basement rocks using well logging data The depth to the basement rocks along the Gulf of Suez from north to south direction was traced, using the most recent wells data (EGPC, 2000) drilled to the basement rocks in the stud- ied area. Figure 15 indicates that the depth to the basement rocks varies widely from less than 2 km at the north western part to more than 3 km at the eastern part of the Gulf. This may be due to the controlling structures related to Gulf of Suez tectonics. These results conform with the RTP magnetic, gravity and geologic cross- section showing the wide variations of the depth along the Gulf of Suez. Thus it is clear that the Gulf of Suez consists of different fault- ed blocks forming the uplifted basement, and the deep basins. This emerges from the results Fig. 16. Interpretation along seismic line EG-31 western side of Gulf of Suez. 1335 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data of the depth estimation methods, negative and positive gravity anomalies. These structures have a great importance for oil production and explo- ration processes. 10. Seismic interpretations Two seismic sections EG-31 and MP-70 were prepared by EGPC trending north west- south east and north east-south west direc- tions respectively (figs. 16 and 17). These lines were correlated with the vertical veloci- ty logs to determine the exact horizons (mark- er beds) and their depth in the studied area. These seismic sections were used to illus- trate the upward continuation and extension of the previously interpreted subsurface struc- tures in the sedimentary sequences. 10.1. Interpretation along line EG-31 This line (fig. 16) lies in the western side of the Gulf of Suez and extends for about 20 km in a northwest-southeast direction. It can Fig. 17. Interpretation along seismic line MP-70 along Gulf of Suez. 1336 Taha Rabeh be noticed that most of the interpreted struc- tures are older than South Gharib Fm. where they extended upward and crossed the older formations (Rudies and Kareem formations), whereas the older structures did not extend upward to cross the younger formations (Be- layim and South Gharib formations). The pres- ence of these structures in large numbers indi- cate that the studied area is affected by active tectonic forces. These structures have been re- juvenated as indicated by their extensions in different time rock units. 10.2. Interpretation along line MP-70 The length of this line (fig. 17) reaches about 12 km and extends in a north east-south west direction. This line is affected by four fault lines. It seems that Zeit Fm. (top of the sedimentary sequence) is affected at its lower surface by the interpreted faults, while the older formations (Rudeis, Kareem, Belayim and South Gharib formations) are crossed by these faults. It is clear that the studied area was affect- ed by the tectonics related to opening of the Gulf of Suez, and the Red Sea. The faults di- vide the studied area into several numbers of major blocks. These faults are very important in sealing oil traps (e.g., Middle Miocene for- mations). The depth to these blocks ranges from 0 km to about 3 km. 11. Discussion and conclusions This study was carried out using different geophysical tools (aeromagnetic, land magnet- ic, gravity and seismic methods). The data were analysed using the most advanced and suitable techniques. The results of tectonic trend analy- sis indicate that most of the predominant direc- tions are north 35°-45° west direction that relat- ed to Gulf of Suez and Red Sea tectonics, and the least predominant direction is north 10°-25° east, that related to Gulf of Aqaba tectonics. The application of Euler deconvolution with a suit- able sliding window to estimate the position of the subsurface intruded magnetic bodies indi- cates that the depth to these bodies ranges from 0 km to about 3 km and the magnetic suscepti- bility ranges between 1 to 3 SI units. Further- more, application of the analytical signal method confirms the result deduced by the Euler deconvolution technique. The 2.5D tech- nique prints out a complete view for the geome- try of the basement surface where the depth ranges between 0 km to more than 3 km. For the upper investigations of the sedimentary se- quence, a geological cross-section was con- structed using the well logging data drilled in the studied area as well as two seismic lines. The results indicate that the structures of the area are related to Gulf of Suez, and Red Sea tectonics. These structures extend from the basement rocks upward into the sedimentary sequences and divide the area into several major faulted blocks that have great important in seal- ing the oil traps. This work is considered a guide line for further exploration process and geo- physical works. REFERENCES ATCHUTA RAO, D., H.V. RAN BABU and P.V. SANKER NARAYAN (1981): Interpretation of magnetic anom- alies due to dikes. The complex gradient method, Geophysics, 46, 1572-1578. COOPER, G.R.J. and D.R. COWAN (2003): Application of fractional calculus to potential field data, Explor. Geophys., 34, 51-56. DAVIS, J.C. (1973): Statistical and Data Analysis in Geol- ogy (Jhon Wily and Sons Inc.), pp. 550. EGPC (Egyptian General Petroleum Cooperation) (1990): Aeromagnetic Map for Sinai Peninsula (scale 1:100 000). EGSMA (Egyptian Geological Survey and Mining Au- thority) (1993): Geologic Map Sinai, Egypt (scale 1:100 000). GM-SYS (1995): GM-SYS®, Gravity and Magnetic Mod- elling Version 3.6 (Northwest Geophysical Associa- tion Inc., Corvallis, Oregon, U.S.A.). GRANT, F.S. and G.F. WEST (1965): Interpretation Theory in Applied Geophysics (McGraw-Hill Book Co., New York), 179-191. HANSEN, R.O. and R.S. PAWLOWSKI (1989): Reduction to the pole at low latitudes by weiner filtering, Geo- physics, 54, 1607-1613. HENDERSON, R.G. (1966): A comprehensive system of au- tomatic computer in magnetic and gravity interpreta- tion, Geophysics, 25, 569-585. HSU, S.-K. (2002): Imaging magnetic sources using Euler’s equation, Geophys. Prospect., 50, 15-25. IPRG (Institute for Petroleum Research and Geophysics) (ca. 1982): Bouguer Anomaly Map of Sinai (scale 1:100 000). 1337 Structural set-up of Southern Sinai and Gulf of Suez areas indicated by geophysical data LINDQUIST, S.J. (1998): The Red Sea Basin Province: Su- dr-Nubia(!) and Maqna(!) petroleum systems, U.S. Geol. Surv. Open-File Rep. 99-50-A. NABIGHIAN, M.N. (1972): Additional comment of the ana- lytic signal of two-dimensional magnetic bodies with polygonal cross-section, Geophysics, 39, 85-92. NABIGHIAN, M.N. (1984): Toward the three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms. Fundamental rela- tions, Geophysics, 53, 957-966. PARASNIS, D.S. (1997): Principles of Applied Geophysics (Chapman and Hall Publisher, London, U.K.), 5th edition, 21-25. PATTON, T.L., A.R. MOUSTAFA, R.A. NELSON and S.A. ABDINE (1994): Tectonic evolution and structural setting of the Suez Rift, in Interior Rift Basins: Tulsa, OK, edited by S.M. LANDON, Am. Ass. Pet. Geol. Mem., 59, 9-55. REID, A.B., J.M. ALLSOP, H. GRANSER, A.J. MILLETT and I.W. SOMERTON (1990): Magnetic interpretation in three dimensions using Euler deconvolution, Geo- physics, 55 (1), 80-91. ROEST, W.E., J. VERHOEF and M. PILKINGTON (1992): Mag- netic interpretation using 3D analytic signal, Geo- physics, 57, 116-125. STAVREV, P.Y. (1997): Euler deconvolution using differen- tial similarity transformations of gravity or magnetic anomalies, Geophys. Prospect., 45, 207-246. SURFER (1999): Surfer® Version 7.00, Feb 28, 2000 (Gold- en Software Inc., Golden, CO, U.S.A.). SURFER (2002): Surfer® Version 8.00, Feb 11, 2002 (Gold- en Software Inc., Golden, CO, U.S.A.). THOMPSON, D.T. (1982): EULDPH – A new technique for making computer-assisted depth estimates from mag- netic data, Geophysics, 47, 31-37. (received June 17, 2003; accepted December 3, 2003)