Hasil Pencarian  ::  Simpan CSV :: Kembali

"Metode inversi seismik deterministik merupakan salah satu metode yang digunakan dalam proses karakterisasi reservoar. Metode ini digunakan karena ia dapat menggambarkan kondisi bawah permukaan melalui nilai impedansi akustik yang dihasilkannya. Akan tetapi, nilai impedansi akustik yang dihasilkan oleh inversi seismik deterministik dinilai kurang representatif dalam proses karakterisasi reservoar karena dalam proses inversinya hanya dihasilkan satu realisasi saja. Oleh karena itu, dikembangkan metode inversi stokastik yang menggunakan konsep geostatistik. Di mana ia akan menghasilkan beberapa realisasi nilai impedansi akustik untuk meminimalisir ketidakpastian. Penelitian ini memanfaatkan metode inversi geostatistik untuk proses karakterisasi reservoar pada lapangan “X” di Laut Utara. Pada penelitian ini digunakan data seismik 3D post-stack migration dan dua sumur. Melalui penerapan inversi geostatistik di zona target pada lapangan “X” diketahui bahwa pada zona target, litologinya merupakan batu pasir dan shale yang posisinya berselingan dan tidak terlihat indikasi keberadaan fluida hidrokarbon.

---

Deterministic seismic inversion method is one of the methods used in the reservoir characterization process. This method is used because it can describe the subsurface conditions through acoustic impedance value. However, the value of acoustic impedance produced by deterministic seismic inversion is considered less representative in the reservoir characterization process because the inversion only produces one realization of acoustic impedance. Therefore, a stochastic inversion method was developed. The method used a geostatistical concept where it will result in several realizations of acoustic impedance values to minimize uncertainty. This study utilizes the geostatistical inversion method for the reservoir characterization process in the "X" field in the North Sea. In this study, 3D post-stack migration seismic data and two wells were used. Through the application of geostatistical inversion in the target zone on the "X" field, it is known that in the target zone, the lithology is sandstone and shale which are alternately positioned and there is no indication of the presence of hydrocarbon fluids."

"Cekungan Jawa Timur Utara tersebut merupakan cekungan Tersier yang produktif akan hidrokarbon pada setiap interval stratigrafi (Marianto et al., 2017). Penelitian ini dilakukan pada lapangan B dan  terfokus pada formasi Kujung. Formasi Kujung adalah formasi yang litologinya di dominasi oleh karbonat pada reservoirnya dan memiliki batuan penutup berupa shale. Karakterisasi reservoir karbonat merupakan suatu tantangan sebab batuannya memiliki struktur pori yang kompleks akibat proses diagenesis. Sehingga mempengaruhi kecepatan gelombang seismik yang merambat pada karbonat. Data yang  digunakan dalam penelitian ini adalah data Seismik 3D PSTM dengan kontrol 3 sumur berbeda. Metode yang digunakan untuk mengkarakterisasi reservoir karbonat ini adalah dengan menggunakan metode seismik inversi simultan. Metode Seismic Inversi Simultan tersebut akan menghasilkan model Impedansi-P (AI), Impedansi-S (SI) dan densitas yang kemudian akan di transformasikan menjadi Lambda-Rho dan Mu-Rho. Model dari parameter hasil inversi simultan dan hasil transformasi tersebut akan digunakan untuk mengetahui sebaran litologi, porositas dan konten fluida pada pori batuan. Berdasarkan hasil analisis, daerah prospek pada penelitian ini terletak pada build up karbonat serta memiliki porositas yang tinggi dengan nilai AI sebesar 16000-31000 ((ft/s)(g/cc)) dan Mu-rho sebesar 10-27 (GPA g/cc) dan fluida minyak dengan nilai λÏ? 10-36 (GPA g/cc.

---

North East Java Basin is a tertiary basin which prolific of hydrocarbon on almost all stratigraphic intervals (Marianto et al., 2017). The study area is located on B field and focused on Kujung formation. Kujung formation is dominated by carbonate on its reservoir and has shale as its caprock or seal. Reservoir characterization considered as a huge challenge because of its complex structural pore caused by diagenesis processes. Therefore, affect the seismic wave propagation which travels through carbonate body. This study use 3D Seismic PSTM with 3 different wells as a control. The method used in this study is Simultaneous Seismic Inversion. The Seismic Simultaneous Inversion will extract P-Impedance (AI), S-Impedance (SI) and density as output models. Those three models will be transformed into Lambda-Rho dan Mu-Rho. The Output models of inversion and transformation will be used for characterizing distribution of the litology, porosity and fluid content that fills rocks pore. Based on analysis result, the prospect area in this study located on carbonate build up, high porosity indicated by AI value between 16000-31000 ((ft/s)(g/cc)) and Mu-rho value between 10-27 (GPA g/cc) and oil by value between 10-36 (GPA g/cc)."

"[ABSTRAKLapangan Athar mulai berproduksi sejak tahun 1975 dengan produksi kumulatif mencapai 900 MMbbls dan RF 50%. Pada tahun 2011 dilakukan akuisisi seismik 3D untuk melihat remaining potential yang ada di lapangan ini. Analisa seismik 3D dengan atribut seismik dan inversi simultan memberikan hasil yang cukup baik untuk menentukan penyebaran batupasir dan hidrokarbon yang ada di dalam nya. Volume densitas hasil inversi, atribut minimum amplitude, dan atribut arc length membantu dalam mengidentifikasi penyebaran reservoar. Lambda-Rho dan AI membantu dalam mengidentifikasi area-area yang mengandung hidrokarbon. Interpretasi struktur di seismik menunjukkan adanya sesar minor, yang sebelumnya tidak teridentifikasi pada seismik 2D. Hasil analisa seismik digunakan dalam pembuatan 3D geomodel. Penentuan batas channel, dan area hidrokarbon di dipandu hasil dari analisa seismik dan data sumur. Empat tubuh channel teridentifikasi dari analisis tersebut, sementara overbank deposit disebarkan secara statistik. Data interpreasi struktur digunakan untuk menentukan jumlah segmen atau kompartemen, dan hasilnya adalah lima segmen selatan dan satu segmen utara terdapat di zona dangkal Lapangan Athar. Hasil pemodelan geologi menunjukkan bahwa segmen 5 menjadi area yang masih memiliki prosepek untuk produksi minyak, sementara segmen 2 memiliki prospek untuk produksi gas.

---

ABSTRACTAthar Field start producing since 1975 up to now with cummulative production reached 900 MMbbls and RF 50%. 3D seismic was acquired in 2011 to identify the remaining potential in this field. Analysis of 3D seismic with attribute and simultaneous inversion provide a good result to determine reservoar distribution and hydrocarbon contained. Density cube from inversion, minimum amplutide attribute, anda arc length attribute are used to map the sand distribution. Lambda Rho and AI (Ip) are usefull to identify the remaining hyrocarbon area. Structural interpretation from seismic shows there are minor faults which is not identified before with 2D seismic. The result of the analysis was used to create 3D Geomodel. The channel limit determination, and remaining hydrocarbon area guided by analysis from seismic and well data. As a result, four channel bodies were identified, whilst overbank deposit distributed statistically. Structural interpretation data used to determine the number of segment or compartment, and the result are five segments in the south and one segment in the north identified in the shallow zone Athar Field. The result of geological modeling shows that segment 5 still has prospect in oil production, whilst segment 2 has prospect in gas production.;Athar Field start producing since 1975 up to now with cummulativeproduction reached 900 MMbbls and RF 50%. 3D seismic was acquired in 2011 toidentify the remaining potential in this field.Analysis of 3D seismic with attribute and simultaneous inversion provide agood result to determine reservoar distribution and hydrocarbon contained. Densitycube from inversion, minimum amplutide attribute, anda arc length attribute are usedto map the sand distribution. Lambda Rho and AI (Ip) are usefull to identify theremaining hyrocarbon area. Structural interpretation from seismic shows there areminor faults which is not identified before with 2D seismic.The result of the analysis was used to create 3D Geomodel. The channel limitdetermination, and remaining hydrocarbon area guided by analysis from seismic andwell data. As a result, four channel bodies were identified, whilst overbank depositdistributed statistically. Structural interpretation data used to determine the number ofsegment or compartment, and the result are five segments in the south and onesegment in the north identified in the shallow zone Athar Field.The result of geological modeling shows that segment 5 still has prospect inoil production, whilst segment 2 has prospec in gas production, Athar Field start producing since 1975 up to now with cummulativeproduction reached 900 MMbbls and RF 50%. 3D seismic was acquired in 2011 toidentify the remaining potential in this field.Analysis of 3D seismic with attribute and simultaneous inversion provide agood result to determine reservoar distribution and hydrocarbon contained. Densitycube from inversion, minimum amplutide attribute, anda arc length attribute are usedto map the sand distribution. Lambda Rho and AI (Ip) are usefull to identify theremaining hyrocarbon area. Structural interpretation from seismic shows there areminor faults which is not identified before with 2D seismic.The result of the analysis was used to create 3D Geomodel. The channel limitdetermination, and remaining hydrocarbon area guided by analysis from seismic andwell data. As a result, four channel bodies were identified, whilst overbank depositdistributed statistically. Structural interpretation data used to determine the number ofsegment or compartment, and the result are five segments in the south and onesegment in the north identified in the shallow zone Athar Field.The result of geological modeling shows that segment 5 still has prospect inoil production, whilst segment 2 has prospec in gas production]"

"Karakterisasi reservoar pada Lapangan "X" telah dilakukan menggunakan inversi impedansi akustik dan geostatistik. Kedua teknik ini dapat menghasilkan peta persebaran reservoar yang didapatkan dari impedansi akustik hasil inversi dan porositas. Selain itu, peta persebaran reservoar dapat digunakan untuk melihat litologi batuan reservoar. Daerah penelitian terletak di Lapangan "X" pada Cekungan Sumatera Tengah, Provinsi Riau. Teknik inversi impedansi akustik yang digunakan adalah inversi model based, sedangkan teknik geostatistik yang digunakan adalah teknik cokriging. Hasil analisa terintegrasi dari kedua teknik menunjukkan bahwa daerah target mempunyai litologi batupasir ber-porous disekitar well A-1 dengan porositas tinggi pada nilai impedansi akustik 15.000-19592 ((ft/s)*(g/cc)) dan mempunyai litologi batupasir tight disekitar well A-12 dan well A-16 dengan porositas cukup rendah pada nilai impedansi akustik 15306-20204 ((ft/s)*(g/cc)) dan 15306-21429 ((ft/s)*(g/cc)).

---

ABSTRAKReservoir characterization in 'X' field has been carried out using acoustic impedance inversion and geostatistics. Both of these techniques can produce a map of the distribution reservoir that obtained from acoustic impedance inversion and porosity results. Beside that, a map of the distribution reservoir can be used to view the lithology of rocks. The research area is in "X" field at Central Sumatera Basin, Riau, Indonesia. The acoustic impedance inversion is performed by using "model based" inversion, while the cokriging technique had been used as part of geostatistic technique. The integrated analysis results from these methods show that the reservoir in research field has porous sandstones around the well A-1 with high porosity in acoustic impedance values 15000-19592 ((ft/s) *(g/cc)) and has tight sandstone around well A-12 and well A-16 with low porosity in acoustic impedance values 15306-20204 ((ft/s) *(g/cc)) and 15306-21429 ((ft/s)*(g/cc))."

"Lapangan A merupakan lapangan offshore yang berada di Laut mid-Norway dan telah memproduksi minyak dan gas bumi cukup lama. Tiga buah sumur digunakan pada penelitian ini dengan kelengkapan data Gamma Ray, NPHI, RHOB, dan log Sonik. Untuk mengidentifikasi reservoir hidrokarbon digunakan metode seismik inversi untuk menghasilkan nilai impedansi akustik dan metode well log berupa gamma ray dan porositas guna menggambarkan penyebaran batu pasir pada reservoir hidrokarbon. Dalam proses inversi dibuat dua buah horison yaitu top dan bottom dengan kontrol tiga buah sumur.

---

Field A is an offshore field, located in the mid-Norway Sea and has been producing oil and gas for a while. Three wells were used in this study that have Gamma Ray, NPHI, RHOB, and sonic logs. Seismic inversion method is used to identify the hydrocarbon reservoir then to deliver an acoustic impedance. Well log analysis such as gamma ray, porosity, and density are used to illustrate lithology of hydrocarbon reservoir. Inversion method created two horizons ,the top and bottom, with three wells controls."

"ABSTRAKDaerah penelitian merupakan wilayah kerja PT. Pilona Tanjung Lontar yangterletak di barat daya, sub cekungan palembang selatan. Terdapat lapangan minyakA, B, C dan D yang beroperasi hampir 20 tahun. Keempat lapangan ini berproduksiminyak dari Batu Pasir Formasi Muara Enim dan menjadi fokus reservoar padapenelitian ini. Lokasi penelitian yang dinamakan Lapangan B-1, merupakan areaperluasan (?step out?) dari lapangan B. Hasil interpretasi seismik menyatakanbahwa Lapangan B-1 merupakan sayap antiklin dari struktur lapangan B sehinggakecil kemungkinan akumulasi hidrokarbon pada daerah rendahan dan tiada aktivitaslanjut untuk pengembangan. Akan tetapi pada daerah sayap antiklin berpeluangterdapatnya akumulasi hidrokarbon berupa perangkap stratigrafi yang disebabkanperubahan fasies secara bertahap ke arah vertikal yang ditunjukkan oleh indikasiisolated strong amplitude dibeberapa penampang seismik. Kombinasi analisis datasumur, analisis inversi seismik dan analisis AVO ini dilakukan agar dapatmengevaluasi indikasi perangkap stratigrafi dan memetakan prospek hidrokarbonpada area ini. Berdasarkan hasil inversi bandlimited, Lapangan B-1 terdapatanomali lokal yang memiliki nilai impedansi akustik dan densitas yang relatif lebihbesar daripada daerah disekitarnya yang juga berkorelasi dengan hasil analisis AVOyang dicurigai sebagai akumulasi hidrokarbon pada perangkap stratigrafi. Hasilanalisis atribut AVO intercept*gradient didapati anomali nya berupa kelas 2p yangberarti batu pasir tersaturasi dengan minyak.

---

ABSTRACTThe research area is the working area of PT. Pilona Tanjung Lontar is located in thesouthwest, sub basin palembang south. There are oil fields A, B, C and D whichoperates nearly 20 years. The fourth field is producing oil from Muara EnimFormation Sandstone and the focus of the reservoir in this study. Location of thestudy, called Fields B-1, an expansion area ( "step out") from the field B. The resultsof the seismic interpretation states that Fields B-1 is a wing of the anticline structureof the field and making it less likely the accumulation of hydrocarbons in the lowerarea and no further activity for development. On the other side, the presence ofhydrocarbon accumulation at wing area anticline in the form of stratigraphic trapscaused gradual facies changes in the vertical direction indicated by the indicationof some isolated strong seismic amplitude. The combination of well data analysis,analysis of seismic inversion and AVO analysis is performed in order to evaluatethe indication of stratigraphic traps and map hydrocarbon prospects in this area.Based on the results bandlimited inversion, Field B-1 there is a local anomaly thathas an acoustic impedance value and relative density larger than the area around itis also correlated with the results of AVO analysis of suspected hydrocarbonaccumulation in stratigraphic traps. The results of the AVO attribute analysisintercept * gradient anomalies found his form of class 2p means sandstone saturatedwith oil."

"ABSTRAKStudi inversi AVO melalui pendekatan simultaneous inversion di lapangan Alphadidesain untuk mendapatkan model geologi dan geofisika di lapangan tersebut,termasuk didalamnya distribusi dan geometri reservoar serta identifikasi fluidahidrokarbon. Simultaneous inversion dilakukan pada data seismik angle stack danditujukan untuk menghitung parameter elastis seperti densitas, Vp, Vs, impedansiakustik (Zp) dan impedansi shear (Zs). Parameter tersebut selanjutnya digunakanuntuk karakterisasi reservoar deltaik dan mendelineasi hidrokarbon melaluipendekatan lambda mu rho. Sebagai komparasi, dilakukan pula analisis atributberbasis amplitude pada data post stack dan angle-stackLapangan Alpha merupakan bagian dari Sanga-Sanga PSC dan terletak dibagianselatan cekungan Kutai. Tebal reservoar batupasir umumnya kurang dari 10 meterdan interbedded dengan coal. Seismically, batupasir dibawah resolusi tuningthickness serta refleksi coal yang kuat berinterferensi dengan data seismik post-stackkonvensional.Penampang impedansi akustik (Zp) turunan dari metoda ini telah memberikanpemahaman yang lebih baik mengenai konektivitas, kontinyuitas, dan kompleksitasreservoar deltaik secara vertikal, terutama pada area-area disekitar kontrol sumur.Namun demikian, hasil simultaneous inversion secara umum menunjukkan bahwametoda LMR tidak memberikan separasi litologi dan identifikasi kandungan fluidasesuai yang diharapkan. Sebaliknya, analisis atribut berbasis amplitude pada datapost-stack dan far-angle stack memberikan hasil interpretasi yang lebih memuaskan.

---

AbstractThe Alpha field simultaneous AVO inversion study was designed to providegeological and geophysical model in the field including reservoir distribution andgeometry, and fluid potential identification. Simultaneous inversion was performed inseismic angle stack volumes to obtain elastic parameters such as density, Vp, Vs,acoustic impedance (Zp) and shear impedance (Zs). Those parameters are used tocharacterize reservoir and delineate hydrocarbon throughout lambda mu rhoapproaches. Moreover, amplitude-based attribute analysis also done to compare andcomplete the analysis.The study area is located in the sanga-sanga PSC to the south of Kutai basin. Thereservoir sands, statistically are less than 10 meters and inter-bedded with coals.Seismically, the reservoir sands are below seismic tuning thickness resolution andstrong coal reflections interfere with the conventional post-stack seimic data.The acoustic impedance section derived from this study has been giving betterunderstanding relates to reservoir connectivity, continuity, and complexity of thedeltaic reservoir, mainly in the area which closes to the well control. However, resultsof inversion in general cannot meet the expectation of giving good lithologyseparation and fluid content identification. On the other hand, amplitude-basedattribute analysis on post-stack and far-angle stack volume gives more specific resultin terms of reservoir geometry and fluid identification. "

"Metode seismik merupakan suatu metode geofisika yang hingga saat ini merupakan metode dipercaya dapat memberikan gambaran bawah permukaan dari suatu lapangan. Seismik inversi adalah salah satu dari banyak metode yang digunakan untuk karakterisasi reservoar. Dengan menggunakan inversi seismik simultan mampu mendapatkan jawaban yang lebih pasti dengan cara menganalisa impedansi P, impedansi S dan densitas. Lapangan X, pada zona dangkal merupakan zona produksi yang cukup menjanjikan. Sehingga, pengembangan pada zona ini cukup gencar. Namun, beberapa kali prediksi reservoar meleset dikarenakan lapisan coal yang melimpah. Sebelum melakukan analisa inversi seismik simultan, perlu diketahui terlebih dahulu hubungan sifat fisik batuan (Zp, Zs, VpVs, lamda-rho dan mhu-rho) dengan properti batuannya (densitas dan porositas). Studi kelayakan dilakukan dengan cara melakukan cross plot parameter fisik batuan untuk dapat mendefinisikan litologi dan fluida yang ada. Selanjutnya well-seismic tie dilakukan untuk mendapatkan korelasi dan koherensi antara data dari log sumur dengan data seismik. Interpretasi horison, struktur dan pembuatan model frekuensi rendah dilakukan untuk memahami keadaan geologi dari daerah penelitian. Pada akhirnya, inversi seismik simultan dapat dilakukan dengan tepat. Hasil dari analisa inversi seismik simultan adalah kisaran nilai Zp, Zs dan Dn untuk mengkarakterisasi reservoar. Selain itu, analisa mengenai LMR (lamda-mhurho) juga dilakukan untuk mendapatkan hasil yang lebih baik. Sehingga, persebaran reservoar pada daerah penelitian dapat diketahui.

---

Seismic is one of the methods in geophysics that until now still reliable for sub-surface imaging and interpretation of a field. Seismic inversion is one of the methods to characterize the reservoir. Using simultaneous seismic inversion, analyzing the P-impedance, S-impedance and density can be more convenient and certain.

X Field, especially in the shallow zone is a promising production zone. Hence, development in this zone is very incentive. Unfortunately, miss prediction sometimes happens due to coal layer that abundant. Relationship between rock physic (Zp, Zs, Vp/Vs, lamda-rho and mhu-rho) and rock property (density and porosity) must be understood prior to simultaneous seismic inversion analysis. Feasibility study was conducted by cross plotting among some parameters to define the lithology and fluids.

Well-seismic tie was conducted to have a good correlation and coherency between well-log data and seismic data. Horizon, structural interpretation and low frequency model were performed to have geological understanding of research area. The result of series steps previously then analyzed to have a good quality data. Eventually, the simultaneous seismic inversion can be performed in a proper way.

Result of the simultaneous seismic inversion analysis is a value of Zp, Zs and Dn for reservoir characterization. Moreover, LMR (lamda-mhu-rho) analysis can be performed to give more preferable result. Eventually, the distribution of gas-sand reservoir can be understood."

"ABSTRAKStudi menggunakan metoda seismik telah diterapkan dalam penelitian ini untuk menginterpretasi dan mendelineasi penyebaran reservoar Gita sand di Lapangan S. Dengan menggunakan inversi seismik, kita bisa mendapatkan model geologi penyebaran reservoar melalui penampang dan slicing P-impedance. Lapangan S terletak di sebelah utara Cekungan Asri dan merupakan lapangan eksplorasi yang belum pernah ditemukan potensi hidrokarbon. Hipotesis awal dari analisa sensifitas menunjukkan bahwa data sumur memiliki sensifitas terhadap P-impedance untuk memisahkan litologi. Analisa hasil inversi membantu interpretasi model geologi menjadi lebih akurat dan juga menunjukkan bahwa reservoar tersebar di sekitar prospek eksplorasi Undip B-1. Hasil penelitian membuktikan bahwa prospek Undip B-1 dapat menjadi prospek hidrokarbon baru di lapangan S.

---

AbstractSeismic inversion method has been applied to interpret and delineate reservoir distribution of Gita sand in S field. Using seismic inversion, it can get reservoir distribution model through P-impedance section and slicing. S field is exploration field located in the northern part of Asri Basin and there was no found hydrocarbon potential previously in this field. Initial review from sensitivity analysis shows that well data has sensitivity to p=impedance to classify lithology. The inversion result helps geological model interpretation to be more accurate and also shows resevoir distribution around exploration prospect Updip B-1. The result of this study prove that exploration prospect Updip B-1 can be new hydrocarbon prospect in S field"

"ABSTRAKPenelitian yang dilakukan di cekungan Sumatera Utara bertujuan untuk mengetahhui bagaimana kondisi batuan sedimen pra-tersier di bawah permukaan. Hal ini dilakukan karena adanya isu dimana batuan sedimen yang telah berumur pra-tersier diisukan belum semua nya termertamorfkan. Hal ini menjadi menarik karena ada kemungkinan cadangan minyak dan gas bumi yang tersimpan pada batuan sedimen pra-tersier. Metode gravitasi merupakan metode geofisika yang digunakan untuk memetakan struktur bawah permukaan. Karena nilai yang terukur oleh metode ini adalah nilai anomali gravitasi dari batuan yang ada di bawah permukaan. Metode inversi menggunakan persamaan Fast Fourier Transform pada aplikasi Matlab yang mana bertujuan untuk dapat membuat pemodelan dari bawah permukaan. Data masukan yang digunakan pada pemograman Matlab berupa nilai data anomali graitasi, jumlah kolom dan baris atau koordinat titik pengambilan data, nilai kedalaman rata-rata batuan atau lapisan targer, nilai kontras densitas batuan dan nilai kriteria konvergensi. Hasil pemograman menampilkan horizon bawah permukaan pada nilai kedalaman dan horizon hasil perhitungan nilai kedalaman dengan nilai anomali gravitasi. Proses inversi dilakukan pada saat melakukan proses iterasi dan akan berhenti jika nilai kesalahan RMS yang didapat dari proses iterasi telah mendekati nilai yang ditentukan sebagai kriteria konvergensi atau sampai iterasi maksimum tercapai. Hasil iterasi berhenti saat nilai RMS error sekitar 0.00016879 % pada iterasi kedua. Horizon yang terbentuk sesuai dengan keadaan geologi yang terbentuk di bawah permukaan dimana batuan sedimen pra-tersier berada di kedalaman rata-rata 5.8 km di bawah permukaan dengan ketebalan dari 5 km hingga 6 km.

---

ABSTRACTThe research conducted in the North Sumatra basin aims to find out how the conditions of pre-tertiary sedimentary rocks below the surface. This is done because there are issues where sedimentary rocks that have been pre-tertiary are not yet included. This is interesting because it is possible that oil and gas reserves are stored in pre-tertiary sedimentary rocks. The gravity method is a geophysical method used to map the subsurface structure. Because the value measured by this method is the value of the gravitational anomaly of the rock below the surface. The inversion method uses the Fast Fourier Transform equation in the Matlab application, which aims to be able to make modeling from the subsurface. Input data used in Matlab programming is in the form of gravity anomalous data values, column number and row or coordinates of data retrieval points, average rock depth or targer layer values, rock density contrast values and convergence criterion values. The programming results show the subsurface horizon at the depth and horizon values resulting from the calculation of the depth value with the value of gravity anomalies. The inversion process is carried out when carrying out the iteration process and will stop if the RMS error value obtained from the iteration process has approached the value specified as the convergence criterion or until the maximum iteration is reached. The iteration results stop when the RMS error value is around 0.00016879% in the second iteration. Horizon is formed in accordance with the geological conditions formed below the surface where pre-tertiary sedimentary rocks are at an average depth of 5.8 km below the surface with a thickness of 5 km to 6 km."