Studi ini dilakukan dengan maksud dan tujuan untuk melakukan analisis Natural Coke (NC) akibat Intrusi Batuan Beku. Metode pemodelan yang dilakukan dalam studi ini adalah menggunakan metode geostatistika indicator kriging dengan bantuan software Isatis.neo sebagai tools. Kriging adalah metode geostatistik yang digunakan untuk memperkirakan nilai pada titik yang belum terukur berdasarkan hubungan spasial antara titik sampel. Metode ini menggunakan variogram untuk mengukur korelasi spasial. Indicator kriging yang tidak memerlukan asumsi normalitas sangat cocok untuk data biner. Studi ini dilakukan dengan studi kasus pada Site B di PT XYZ. Probabilitas persebarannya ditinjau korelasinya dengan formasi geologi sehingga dapat diketahui penyebab terbentuknya NC di wilayah studi. Probabilitas persebaran NC dan IGN terjadi di Formasi Batupasir Haloq (Teh). Formasi geologi ini dianalisis dan dikorelasikan dengan data litologi lubang bor hasil eksplorasi. Dengan demikian, dapat diketahui bahwa intrusi batuan beku memengaruhi terbentukmya NC yang dibuktikan dengan adanya intrusi magma dan formasi vulkanik di sekitar formasi dengan probabilitas keterdapatan NC > 0,5.