Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Simulation of photoelectric log in oil-bearing formation using artificial neural network■Simulation of photoelectric log in oil-bearing formation using artificial neural network■12082208FAAliMoradzadehFaculty of Mining and Geophysics, Shahrood University of Technology, Shahrood, IranElhamBakhshiFaculty of Mining and Geophysics, Shahrood University of Technology, Shahrood, IranJournal Article20050226Estimating of various petrophysical parameters and determining of subsurface geology is very important in petroleum reservoir evaluation. Exploration drilling and various well logs normally provide this sort of information. Among the various well logs, the photoelectric (PEF) log is very important as it able to determine the lithology of the reservoir precisely. Therefore for those wells for which this log is not available it is necessary to predict them somehow. In this study it is aimed to use artificial neural network (ANN) ability to tackle this problem. To achieve the goals, a back-propagation ANN (BP-ANN) is planned to model the interrelationships between seven different well logs, and PEF logs. Data from three wells in the Ahvaz oil field (Asmari reservoir) are organized into training, testing and validation data sets for BP-ANN modeling. Data of the fourth and the fifth wells in the same field are retained as independent data sets for evaluating the ability of the network PEF prediction. Once the designed network has been trained properly, its performance has also been tested. When it has been found that the performance is satisfactory the data set of the fourth and the fifth wells are applied to the trained network. The results of the ANN modeling show that the designed network with three layers and architecture of 7-10-1 can produce the precise PEF log that compares well with the measured PEF logs. This means that the designed network is capable enough to predict the PEF logs for the required wells in the same area.Estimating of various petrophysical parameters and determining of subsurface geology is very important in petroleum reservoir evaluation. Exploration drilling and various well logs normally provide this sort of information. Among the various well logs, the photoelectric (PEF) log is very important as it able to determine the lithology of the reservoir precisely. Therefore for those wells for which this log is not available it is necessary to predict them somehow. In this study it is aimed to use artificial neural network (ANN) ability to tackle this problem. To achieve the goals, a back-propagation ANN (BP-ANN) is planned to model the interrelationships between seven different well logs, and PEF logs. Data from three wells in the Ahvaz oil field (Asmari reservoir) are organized into training, testing and validation data sets for BP-ANN modeling. Data of the fourth and the fifth wells in the same field are retained as independent data sets for evaluating the ability of the network PEF prediction. Once the designed network has been trained properly, its performance has also been tested. When it has been found that the performance is satisfactory the data set of the fourth and the fifth wells are applied to the trained network. The results of the ANN modeling show that the designed network with three layers and architecture of 7-10-1 can produce the precise PEF log that compares well with the measured PEF logs. This means that the designed network is capable enough to predict the PEF logs for the required wells in the same area.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Dating loess by red thermoluminescence of quartz■Dating loess by red thermoluminescence of quartz■213382210FAMortezaFattahiInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, Iran / OUCE, University of Oxford, Mansfield Road, OX1 3TB0000-0003-2719-5859Journal Article20050718The loess deposits in central Asia, China, Europe and the USA provide records of past climate, which permits the study of paleoclimatology. Although luminescence dating has had an important role in dating loess, significant challenges still remain.
UV-blue emission from polymineral and feldspar suffer from anomalous fading. UV-blue emission from quartz saturates at c. <500 Gy. Fattahi (2001) reported the ability of red thermoluminescence for dating volcanic quartz to c.1.3 Ma.
This paper investigates the possibility of applying single aliquot regeneration red thermoluminescence (SAR RTL) protocol for dating loess and demonstrating that:
1- It is possible to detect red thermoluminescence from loess.
2- Red thermoluminescence from loess is bleachable.
3- SAR RTL is able to recover a known laboratory dose.
4- Equivalent doses around 1000 Gy are accessible.
This information suggests that SAR RTL is a powerful protocol for dating loess which can extend the time range of luminescence dating of unburnt Aeolian sediments.The loess deposits in central Asia, China, Europe and the USA provide records of past climate, which permits the study of paleoclimatology. Although luminescence dating has had an important role in dating loess, significant challenges still remain.
UV-blue emission from polymineral and feldspar suffer from anomalous fading. UV-blue emission from quartz saturates at c. <500 Gy. Fattahi (2001) reported the ability of red thermoluminescence for dating volcanic quartz to c.1.3 Ma.
This paper investigates the possibility of applying single aliquot regeneration red thermoluminescence (SAR RTL) protocol for dating loess and demonstrating that:
1- It is possible to detect red thermoluminescence from loess.
2- Red thermoluminescence from loess is bleachable.
3- SAR RTL is able to recover a known laboratory dose.
4- Equivalent doses around 1000 Gy are accessible.
This information suggests that SAR RTL is a powerful protocol for dating loess which can extend the time range of luminescence dating of unburnt Aeolian sediments.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Optimization of the source and receiver array in 2D acquisition of the Hoveizeh oil field■Optimization of the source and receiver array in 2D acquisition of the Hoveizeh oil field■354082215FAMohammadMorshediInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, IranMohammad AliRiahiInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, IranJournal Article20060118In seismic data acquisition, source and receiver array is used for attenuation of surface waves caused by the source. On the other hand regarding the spatial sampling, to prevent spatial aliasing, the wavelength should not be shorter than twice the trace interval. Therefore by designing a suitable source and receiver array one may suppress waves that are prone to aliasing. Thus, to further attenuate these phenomena (surface waves and aliasing) will be necessary to optimize arrays.
One of the ways to optimize source and receiver arrays is spatial convolution. In this method each complex array consists of several simple arrays which make a large array when joined together. By using simple arrays of 2 and 3 elements and putting the notch points of these arrays in side lobes of the final response of the array, the side lobes can be reduced. Considering that both source and receiver arrays can affect the amplitude of the generated waves, and the total effect of these arrays equates the summation of both effects in decibel, it is also possible to compute the effect of the source array on the final response of source and receiver arrays in order to optimize the final response of array.In seismic data acquisition, source and receiver array is used for attenuation of surface waves caused by the source. On the other hand regarding the spatial sampling, to prevent spatial aliasing, the wavelength should not be shorter than twice the trace interval. Therefore by designing a suitable source and receiver array one may suppress waves that are prone to aliasing. Thus, to further attenuate these phenomena (surface waves and aliasing) will be necessary to optimize arrays.
One of the ways to optimize source and receiver arrays is spatial convolution. In this method each complex array consists of several simple arrays which make a large array when joined together. By using simple arrays of 2 and 3 elements and putting the notch points of these arrays in side lobes of the final response of the array, the side lobes can be reduced. Considering that both source and receiver arrays can affect the amplitude of the generated waves, and the total effect of these arrays equates the summation of both effects in decibel, it is also possible to compute the effect of the source array on the final response of source and receiver arrays in order to optimize the final response of array.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Comparison of numerical integration methods in orbit determination of low earth orbiting satellites■Comparison of numerical integration methods in orbit determination of low earth orbiting satellites■415782217FAMehdiEshaghRoyal Institute of Technology (KTH), SE 100 44 Stockholm, Sweden and Islamic Azad University, Shahr-e-Rey branch,
P.O. Box 18735-334, Tehran, IranMehdiNajafi AlamdariK. N. Toosi University of Technology, P.O. Box 15875-4416, Tehran, IranJournal Article20030831Comparison of some numerical integration methods of solving the differential equation of motion of a satellite is the main subject of this paper. Since the equation of motion of a satellite is a second order differential equation, therefore, six initial values should be introduced to the numerical solution. These six initial values are the components of position and velocity vectors in an inertial frame respectively. Comparing numerically integrated position and velocity vectors with Keplerian orbit; one can obtain the bias of the numerical integration method in a satellite-centered coordinate system. In this research, three methods of Runge-Kutta, Runge-Kutta-Nystrom, and the predictor-corrector method of Adams-Bashforth and Adams-Moulton are investigated for a low earth orbiting satellite. Numerical results show that with integration size of 30 seconds, the Runge-Kutta method, Adams-Bashforth and Adams-Moulton predictor-corrector algorithms, and Runge-Kutta-Nustrom provide closer orbit to the theoretical orbit respectively.Comparison of some numerical integration methods of solving the differential equation of motion of a satellite is the main subject of this paper. Since the equation of motion of a satellite is a second order differential equation, therefore, six initial values should be introduced to the numerical solution. These six initial values are the components of position and velocity vectors in an inertial frame respectively. Comparing numerically integrated position and velocity vectors with Keplerian orbit; one can obtain the bias of the numerical integration method in a satellite-centered coordinate system. In this research, three methods of Runge-Kutta, Runge-Kutta-Nystrom, and the predictor-corrector method of Adams-Bashforth and Adams-Moulton are investigated for a low earth orbiting satellite. Numerical results show that with integration size of 30 seconds, the Runge-Kutta method, Adams-Bashforth and Adams-Moulton predictor-corrector algorithms, and Runge-Kutta-Nustrom provide closer orbit to the theoretical orbit respectively.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023New estimations for the radius of geo-stationary satellite orbits based on spherical and ellipsoidal gravitational models■New estimations for the radius of geo-stationary satellite orbits based on spherical and ellipsoidal gravitational models■596782220FAAli RezaAzmodeh ArdalanDepartment of Surveying and Geomatics Engineering, Center of Excellence of Surveying and Disaster Management, University of TehranTolouSilaviDepartment of Surveying and Geomatics Engineering, Center of Excellence of Surveying and Disaster Management, University of TehranJournal Article20041121Geo-stationary orbit has been in use since 1964 for telecommunication purposes. Our investigation on the orbit of such satellites revealed that the radius of geo-stationary satellites are so far computed based on the point mass W=GM/r gravitational model of the Earth. To investigate the effect of more improved Earth’s gravitational models on the radius of the Geo-stationary satellites we have considered the following models: (i) Bjerhammar gravitational potential field, (ii) first term of ellipsoidal harmonic expansion of the Earth’s gravitational field, (iii) Somigliana-Pizzeti gravitational potential field. According the results of the computations, the radius of the Geo-stationary satellite computed based on the aforementioned gravitational potential models are deviated from the radius computed based on point mass model W = GM/r, by (i) 0.5km, (ii) 2.15km, and (iii) 2.7km, respectively. Since such an error in the radius of geostationary satellites results in small movements of the satellite in orbit that justifies why the current geo-stationary satellites require fuel and racket engine to fine tune the position of the satellite in its orbit at the specific interval of time. These orbital maneuvers are so important that even the life time of those satellites is determined from the time span that their fuel allows for the fine tuning of the orbit. Based on the computations, the radius of 42, 161, 465.71m, computed based on the Somigliana-Pizzeti gravitational potential field is recommended for the placement of geo-stationary satellites.Geo-stationary orbit has been in use since 1964 for telecommunication purposes. Our investigation on the orbit of such satellites revealed that the radius of geo-stationary satellites are so far computed based on the point mass W=GM/r gravitational model of the Earth. To investigate the effect of more improved Earth’s gravitational models on the radius of the Geo-stationary satellites we have considered the following models: (i) Bjerhammar gravitational potential field, (ii) first term of ellipsoidal harmonic expansion of the Earth’s gravitational field, (iii) Somigliana-Pizzeti gravitational potential field. According the results of the computations, the radius of the Geo-stationary satellite computed based on the aforementioned gravitational potential models are deviated from the radius computed based on point mass model W = GM/r, by (i) 0.5km, (ii) 2.15km, and (iii) 2.7km, respectively. Since such an error in the radius of geostationary satellites results in small movements of the satellite in orbit that justifies why the current geo-stationary satellites require fuel and racket engine to fine tune the position of the satellite in its orbit at the specific interval of time. These orbital maneuvers are so important that even the life time of those satellites is determined from the time span that their fuel allows for the fine tuning of the orbit. Based on the computations, the radius of 42, 161, 465.71m, computed based on the Somigliana-Pizzeti gravitational potential field is recommended for the placement of geo-stationary satellites.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Synoptic and dynamic analysis of blocking systems, manner of diagnosis of blocking systems and effects of these systems over the Iranian regionSynoptic and dynamic analysis of blocking systems, manner of diagnosis of blocking systems and effects of these systems over the Iranian region698982221FAFaridehHabibiInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, Iran0000-0002-3831-174XJournal Article20041101This research is an attempt to investigate the blocking systems over the Atlantic Ocean and especially over the Mediterranean sea, which leads to surface cyclogenesis in the east of the Mediterranean Sea, and the influence of their effects over the Iranian region. At the first step of study, Mediterranean topography, local Mediterranean winds and cyclogenesis have studied; then structure and types of existing blocking systems on the synoptic weather maps have been introduced with figures.
In the second step of the study, all of the blocking systems that were formed over the Atlantic Ocean were identified from 1989 to 1997, the results of which are shown in the table included in the paper. Of course this part of the work was done with the help of Prof. Colucci from USA and data are given from NCEP.
Atmospheric blocking systems can lead to a stagnation of weather patterns where the patterns remain for several days or even weeks at the same location during the blocking system. This case can lead to flooding, drought, above normal temperatures, below normal temperatures and other weather extremes. Therefore, it is important to recognize a blocking pattern in its initial development. With this awareness, one can forecast up to several days in advance with a high degree of accuracy.This research is an attempt to investigate the blocking systems over the Atlantic Ocean and especially over the Mediterranean sea, which leads to surface cyclogenesis in the east of the Mediterranean Sea, and the influence of their effects over the Iranian region. At the first step of study, Mediterranean topography, local Mediterranean winds and cyclogenesis have studied; then structure and types of existing blocking systems on the synoptic weather maps have been introduced with figures.
In the second step of the study, all of the blocking systems that were formed over the Atlantic Ocean were identified from 1989 to 1997, the results of which are shown in the table included in the paper. Of course this part of the work was done with the help of Prof. Colucci from USA and data are given from NCEP.
Atmospheric blocking systems can lead to a stagnation of weather patterns where the patterns remain for several days or even weeks at the same location during the blocking system. This case can lead to flooding, drought, above normal temperatures, below normal temperatures and other weather extremes. Therefore, it is important to recognize a blocking pattern in its initial development. With this awareness, one can forecast up to several days in advance with a high degree of accuracy.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Depth and shape factor determination by successive correlation coefficients of least-square of residual anomalies■Depth and shape factor determination by successive correlation coefficients of least-square of residual anomalies■9110582224FARezaAzadmardanInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, IranVahidEbrahimzadeh ArdestaniInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, IranNasrollahKamalianInstitute of Geophysics, University of Tehran, P.O. Box 14155-6466, Tehran, IranJournal Article20050918In this paper with the method considering the correlation coefficients the order of regional field could be determined.
This method can be used for both Bouguer and residual anomalies.
Capability of this method is tested in the case of small models which have not been considered in the main paper.
Depth estimation is one of the most important aspects in the method.In this paper with the method considering the correlation coefficients the order of regional field could be determined.
This method can be used for both Bouguer and residual anomalies.
Capability of this method is tested in the case of small models which have not been considered in the main paper.
Depth estimation is one of the most important aspects in the method.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023A comparative study of the recent geopotential models for synthesizing different gravity functional at the geographical region of Iran■A comparative study of the recent geopotential models for synthesizing different gravity functional at the geographical region of Iran■10711682226FAAli RezaAzmoodeh ArdalanDepartment of Surveying and Geomatics Engineering, Center of Excellence of Surveying and Disaster Management, University of TehranAbdol RezaSafariDepartment of Surveying and Geomatics Engineering, Center of Excellence of Surveying and Disaster Management, University of Tehran0000-0001-5938-5468YahyaJamourNational Cartographic Center of IranJournal Article20060619Recent progress in the field of satellite, airborne, and surface gravimetry has resulted in geopotential models with high resolution. Since the geopotential models are computed based on dense and highly accurate gravity data, they are valuable sources of information for presentation of long and medium wavelength spectrum of the gravity field in the computational algorithm of modern gravity field modeling techniques. Considering the variety of the geopotential models which are currently available, it is necessary to verify their accuracy in synthesizing the gravity functional so that one can select the best geopotential model for a region of interest. In this study 4 geopotential models, namely: EGM96, PGM2000A, Eigen-cgo1c, Eigen-Grace2s, are compared for their accuracy in synthesizing gravity functional of the types: (1) astronomical longitude, (2) astronomical latitude, (3) norm of gravity acceleration, and (4) geoid from GPS/Leveling in the geographical region of Iran. Based on the numerical results EGM96 and PGM2000A perform almost the same in synthesizing the aforementioned gravity functionals in the geographical region of Iran.Recent progress in the field of satellite, airborne, and surface gravimetry has resulted in geopotential models with high resolution. Since the geopotential models are computed based on dense and highly accurate gravity data, they are valuable sources of information for presentation of long and medium wavelength spectrum of the gravity field in the computational algorithm of modern gravity field modeling techniques. Considering the variety of the geopotential models which are currently available, it is necessary to verify their accuracy in synthesizing the gravity functional so that one can select the best geopotential model for a region of interest. In this study 4 geopotential models, namely: EGM96, PGM2000A, Eigen-cgo1c, Eigen-Grace2s, are compared for their accuracy in synthesizing gravity functional of the types: (1) astronomical longitude, (2) astronomical latitude, (3) norm of gravity acceleration, and (4) geoid from GPS/Leveling in the geographical region of Iran. Based on the numerical results EGM96 and PGM2000A perform almost the same in synthesizing the aforementioned gravity functionals in the geographical region of Iran.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023The effect of solid tide in geopotential field of an elastic and inelastic earthThe effect of solid tide in geopotential field of an elastic and inelastic earth1982232FAMehdiEshaghRoyal Institute of Technology (KTH), SE 100 44, Stockholm, Sweden and Islamic Azad University, Shahr-e-Rey branch, P.O.Box 18735-334, Tehran, IranMehdiNajafi AlamdariK.N.Toosi University of Technology, P.O.Box 15875-4416,Tehran, IranJournal Article20030819In this paper the influence of solid tide on the Earth gravity field is considered. In this consideration the Earth can be regarded as either an elastic or inelastic body. Each one of these elastic and inelastic bodies has two main tidal components, frequency dependent and frequency independent components. In this article how to compute the effect of these components in the Earth’s gravity field is presented. In this investigation, an attempt is made to find out whether the Earth should be regarded as an elastic or inelastic body in practical applications. Computations show equivalent effects on the gravity field due to the elastic and inelastic Earth model. The effect of the frequency dependent component of solid tide due to the elastic and inelastic Earth is much smaller than the frequency independent components. It depends on time and tidal constituents and it should be considered in precise applications. Comparisons between the solid tides due to the elastic and inelastic Earth model show that the inelastic Earth is contracted at poles about 3 mm and expanded at equator about 2.5 mm more than the elastic case.In this paper the influence of solid tide on the Earth gravity field is considered. In this consideration the Earth can be regarded as either an elastic or inelastic body. Each one of these elastic and inelastic bodies has two main tidal components, frequency dependent and frequency independent components. In this article how to compute the effect of these components in the Earth’s gravity field is presented. In this investigation, an attempt is made to find out whether the Earth should be regarded as an elastic or inelastic body in practical applications. Computations show equivalent effects on the gravity field due to the elastic and inelastic Earth model. The effect of the frequency dependent component of solid tide due to the elastic and inelastic Earth is much smaller than the frequency independent components. It depends on time and tidal constituents and it should be considered in precise applications. Comparisons between the solid tides due to the elastic and inelastic Earth model show that the inelastic Earth is contracted at poles about 3 mm and expanded at equator about 2.5 mm more than the elastic case.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Least square minimization for depth and shape determination through Micro-gravity dataLeast square minimization for depth and shape determination through Micro-gravity data112182233FAVahidEbrahimzadeh ArdestaniInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, Iran and Center of Excellence in Geomatics Engineering and Disaster Management, Survey Engineering Department, University of TehranJournal Article20060102The least-square minimization approaches for determination of depth, shape and amplitude coefficient expressed by Abdolrahman et al. (2001) for sphere and horizontal cylinder is used for rectangular prisms as synthetic models with and without random noises. The method is also applied for real sources producing micro-gravity data .The capability of the method is tested and discussed in this paper.The least-square minimization approaches for determination of depth, shape and amplitude coefficient expressed by Abdolrahman et al. (2001) for sphere and horizontal cylinder is used for rectangular prisms as synthetic models with and without random noises. The method is also applied for real sources producing micro-gravity data .The capability of the method is tested and discussed in this paper.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023Crustal seismic anisotropy in the south-central Alborz region using Moho Ps converted phasesCrustal seismic anisotropy in the south-central Alborz region using Moho Ps converted phases233282234FAAhmadSadidkhouyInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, IranGholamJavan DoloeiInternational Institute of Earthquake Engineering and Seismology (IIEES), Tehran, IranMohammad RezaGheitanchiInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, IranJournal Article20060131In this paper to find stable anisotropy direction with minimum error, we attempt to automate the choice of analysis of Moho P<sub>s</sub> phase window by performing a grid search over different windows. We have determined the anisotropy direction in the crust of the Alborz region by using 8 events recorded at three IIEES seismic broadband stations (i. e. DAMV, THKV and CHTH).
Analysis of shear wave splitting of the Ps conversion from the crust-mantle boundary indicates a fast azimuth of anisotropy oriented approximately NE-SW. This strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of mountain block ranges.
We estimate the average directions of anisotropy under three stations, DAMV, THKV and CHTH of 42, 55, 55 degree, respectively. Moreover our results show that the magnitude of anisotropy are 0.23, 0.22, 0.27 seconds, respectively, which are in good agreement with the geology and tectonic setting of this region.In this paper to find stable anisotropy direction with minimum error, we attempt to automate the choice of analysis of Moho P<sub>s</sub> phase window by performing a grid search over different windows. We have determined the anisotropy direction in the crust of the Alborz region by using 8 events recorded at three IIEES seismic broadband stations (i. e. DAMV, THKV and CHTH).
Analysis of shear wave splitting of the Ps conversion from the crust-mantle boundary indicates a fast azimuth of anisotropy oriented approximately NE-SW. This strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of mountain block ranges.
We estimate the average directions of anisotropy under three stations, DAMV, THKV and CHTH of 42, 55, 55 degree, respectively. Moreover our results show that the magnitude of anisotropy are 0.23, 0.22, 0.27 seconds, respectively, which are in good agreement with the geology and tectonic setting of this region.Institute of Geophysics, University of TehranJournal of the Earth and Space Physics2538-371X32320061023High porosity anomaly with good reservoir properties in the lower Fahliyan formation (Neocomian) of Darquain Field (SW Iran) by 3D seismicHigh porosity anomaly with good reservoir properties in the lower Fahliyan formation (Neocomian) of Darquain Field (SW Iran) by 3D seismic333982236FAMehrdadMalekiInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, IranAbdolrahimJavaherianInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, IranIrajAbdollahi FardInstitute of Geophysics, University of Tehran, P.O. Box 14155- 6466, Tehran, IranJournal Article20051116Seismic data from the onshore area of Iran over the Darquain structure is used to investigate stratigraphic and structural seismic anomalies inside the Fahliyan formation (Neocomian). The data consists of a 3D grid seismic cube, acquired in the late 1990's. An integrated seismic interpretation and inversion revealed an outstanding feature in the lower part of the Fahliyan formation. This anomaly has N-S trend, parallel to Darquain structure and the results of seismic studies predict high porosity for the observed anomaly in the lower part of Fahliyan as a major oil reservoir. Recent drillings based on the results of seismic interpretation proved its potential.Seismic data from the onshore area of Iran over the Darquain structure is used to investigate stratigraphic and structural seismic anomalies inside the Fahliyan formation (Neocomian). The data consists of a 3D grid seismic cube, acquired in the late 1990's. An integrated seismic interpretation and inversion revealed an outstanding feature in the lower part of the Fahliyan formation. This anomaly has N-S trend, parallel to Darquain structure and the results of seismic studies predict high porosity for the observed anomaly in the lower part of Fahliyan as a major oil reservoir. Recent drillings based on the results of seismic interpretation proved its potential.