Submit manuscript...
MOJ
eISSN: 2573-2919

Ecology & Environmental Sciences

Research Article Volume 4 Issue 2

Diffusion coefficient fractal dimension for characterizing Shajara reservoirs of the Permo – carboniferous Shajara formation, Saudi Arabia

Khalid Elyas Mohamed Elameen Alkhidir

Department of Petroleum and Natural Gas Engineering, College of Engineering, King Saud University, Saudi Arabia

Correspondence: Prof. Khalid Elyas Mohamed Elameen Alkhidir, Ph.D, Department of Petroleum and Natural Gas Engineering, College of Engineering, King Saud University, Riyadh, Saudi Arabia, Tel 966114679118

Received: March 21, 2019 | Published: April 5, 2019

Citation: Alkhidir KEME. Diffusion coefficient fractal dimension for characterizing Shajara reservoirs of the Permo – carboniferous Shajara formation, Saudi Arabia. MOJ Eco Environ Sci. 2019;4(2):85-90. DOI: 10.15406/mojes.2019.04.00138

Download PDF

Abstract

The quality and assessment of a reservoir can be documented in details by the application of diffusion coefficient. This research aims to calculate fractal dimension from the relationship among diffusion coefficient, maximum diffusion coefficient and wetting phase saturation and to approve it by the fractal dimension derived from the relationship among capillary pressure and wetting phase saturation. In this research, porosity was measured on real collected sandstone samples and permeability was calculated theoretically from capillary pressure profile measured by mercury intrusion contaminating the pores of sandstone samples in consideration. Two equations for calculating the fractal dimensions have been employed. The first one describes the functional relationship between wetting phase saturation, diffusion coefficient, maximum diffusion coefficient and fractal dimension. The second equation implies to the wetting phase saturation as a function of capillary pressure and the fractal dimension. Two procedures for obtaining the fractal dimension have been utilized. The first procedure was done by plotting the logarithm of the ratio between diffusion coefficient and maximum diffusion coefficient versus logarithm wetting phase saturation. The slope of the first procedure=3- Df (fractal dimension). The second procedure for obtaining the fractal dimension was determined by plotting the logarithm of capillary pressure versus the logarithm of wetting phase saturation. The slope of the second procedure=Df -3. On the basis of the obtained results of the fabricated stratigraphic column and the attained values of the fractal dimension, the sandstones of the Shajara reservoirs of the Shajara Formation were divided here into three units. The obtained units from bottom to top are: Lower, Middle and Upper Shajara Diffusion Coefficient Fractal Dimension Units. It was found that fractal dimension increases with increasing grain size and permeability.

Keywords: Shajara reservoirs, Shajara formation, diffusion coefficient fractal dimension, capillary pressure fractal dimension

Introduction

Seismo electric effects related to electro kinetic potential, dielectric permittivity, pressure gradient, fluid viscosity, and electric conductivity was first reported.1 Capillary pressure follows the scaling law at low wetting phase saturation was reported.2 Seismo electric phenomenon by considering electro kinetic coupling coefficient as a function of effective charge density, permeability, fluid viscosity and electric conductivity was reported.3 The magnitude of seismo electric current depends porosity, pore size, zeta potential of the pore surfaces, and elastic properties of the matrix was investigated.4 The tangent of the ratio of converted electric field to pressure is approximately in inverse proportion to permeability was studied.5 Permeability inversion from seism electric log at low frequency was studied.6 They reported that, the tangent of the ratio among electric excitation intensity and pressure field is a function of porosity, fluid viscosity, frequency, tortuosity, fluid density and Dracy permeability. A decrease of seismo electric frequencies with increasing water content was reported.7 An increase of seismo electric transfer function with increasing water saturation was studied.8 An increase of dynamic seismo electric transfer function with decreasing fluid conductivity was described.9 The amplitude of seismo electric signal increases with increasing permeability which means that the seismo electric effects are directly related to the permeability and can be used to study the permeability of the reservoir was illustrated.10 Seismo electric coupling is frequency dependent and decreases expontialy when frequency increases was demonstrated.11 An increase of permeability with increasing pressure head and bubble pressure fractal dimension was reported.12 An increase of geometric and arithmetic relaxation time of induced polarization fractal dimension with permeability increasing was described.13–15

Material and methods

Sandstone samples were collected from the surface type section of the Permo-Carboniferous Shajara Formation, latitude 26˚ 52' 17.4", longitude 43˚ 36' 18" (Figure 1). Porosity was measured on collected samples using mercury intrusion Porosimetry and permeability was derived from capillary pressure data. The purpose of this paper is to obtain diffusion coefficient fractal dimension and to confirm it by capillary pressure fractal dimension. The fractal dimension of the first procedure is determined from the positive slope of the plot of logarithm of the ratio of diffusion coefficient to maximum diffusion coefficient log (D1/2/Dmax1/2) versus log wetting phase saturation (log Sw). Whereas the fractal dimension of the second procedure is determined from the negative slope of the plot of logarithm of log capillary pressure (log Pc) versus logarithm of wetting phase saturation (log Sw). 

The diffusion coefficient can be scaled as

Sw= [ D 1 2 D max 1 2 ] 3Df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadofacaWG3b Gaeyypa0ZaamWaaeaadaWcaaqaaiaadseadaahaaWcbeqaamaalaaa baGaaGymaaqaaiaaikdaaaaaaaGcbaGaamiraiGac2gacaGGHbGaai iEamaaCaaaleqabaWaaSaaaeaacaaIXaaabaGaaGOmaaaaaaaaaaGc caGLBbGaayzxaaWaaWbaaSqabeaacaaIZaGaeyOeI0IaamiraiaadA gaaaaaaa@4873@     (1)

Equation 1 can be proofed from

De= [ D*ϕ*connectivity ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadseacaWGLb Gaeyypa0ZaaSaaaeaadaWadaqaaiaadseacaGGQaGaeqy1dyMaaiOk aiaadogacaWGVbGaamOBaiaad6gacaWGLbGaam4yaiaadshacaWGPb GaamODaiaadMgacaWG0bGaamyEaaGaay5waiaaw2faaaqaaiaadsfa aaaaaa@4D06@     (2)

Where De the effective diffusion coefficient in square micrometer per second; D diffusion coefficient in square micrometer per second; ϕ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiabew9aMbaa@39F2@  porosity (dimensionless); connectivity (dimensionless); and Τ tortuosity (dimensionless).

The diffusion coefficient can be scaled as

D= A 2 [ 2*Tcc ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadseacqGH9a qpdaWcaaqaaiaadgeadaahaaWcbeqaaiaaikdaaaaakeaadaWadaqa aiaaikdacaGGQaGaamivaiaadogacaWGJbaacaGLBbGaayzxaaaaaa aa@41C7@     (3)

Where D the diffusion coefficient in square micrometer per second; A pore throat radius in micrometer; Tcc cole time in seconds. The effective diffusion coefficient can be scaled as

De= A e 2 [ 2*Tcceff ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadseacaWGLb Gaeyypa0ZaaSaaaeaacaWGbbGaamyzamaaCaaaleqabaGaaGOmaaaa aOqaamaadmaabaGaaGOmaiaacQcacaWGubGaam4yaiaadogacaWGLb GaamOzaiaadAgaaiaawUfacaGLDbaaaaaaaa@465B@     (4)

Where De the effective diffusion coefficient in square micrometer per second; Ae effective pore throat radius in micrometer; and Tcceff the effective cole time in seconds. Insert equation 4 into equation 2.

A e 2 [ 2*Tcceff ] = [ D*ϕ*connectivity ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaGaam yqaiaadwgadaahaaWcbeqaaiaaikdaaaaakeaadaWadaqaaiaaikda caGGQaGaamivaiaadogacaWGJbGaamyzaiaadAgacaWGMbaacaGLBb Gaayzxaaaaaiabg2da9maalaaabaWaamWaaeaacaWGebGaaiOkaiab ew9aMjaacQcacaWGJbGaam4Baiaad6gacaWGUbGaamyzaiaadogaca WG0bGaamyAaiaadAhacaWGPbGaamiDaiaadMhaaiaawUfacaGLDbaa aeaacaWGubaaaaaa@56CB@     (5)

The effective pore throat radius can be scaled as

A e 2 =8*k*F MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadgeacaWGLb WaaWbaaSqabeaacaaIYaaaaOGaeyypa0JaaGioaiaacQcacaWGRbGa aiOkaiaadAeaaaa@3FAC@     (6)

Where Ae the effective pore throat radius in micrometer; k permeability in square micrometer; and F the formation resistivity factor (dimensionless). Insert equation 6 into equation 5

[ 8*k*F ] [ 2*Tcceff ] = [ D*ϕ*connectivity ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaWaam WaaeaacaaI4aGaaiOkaiaadUgacaGGQaGaamOraaGaay5waiaaw2fa aaqaamaadmaabaGaaGOmaiaacQcacaWGubGaam4yaiaadogacaWGLb GaamOzaiaadAgaaiaawUfacaGLDbaaaaGaeyypa0ZaaSaaaeaadaWa daqaaiaadseacaGGQaGaeqy1dyMaaiOkaiaadogacaWGVbGaamOBai aad6gacaWGLbGaam4yaiaadshacaWGPbGaamODaiaadMgacaWG0bGa amyEaaGaay5waiaaw2faaaqaaiaadsfaaaaaaa@59F3@     (7)

Equation 7 after simplification and division both sides by will become

[ 4*k*F ] [ Tcceff*ϕ ] = [ D*ϕ*connectivity ] T*ϕ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaWaam WaaeaacaaI0aGaaiOkaiaadUgacaGGQaGaamOraaGaay5waiaaw2fa aaqaamaadmaabaGaamivaiaadogacaWGJbGaamyzaiaadAgacaWGMb GaaiOkaiabew9aMbGaay5waiaaw2faaaaacqGH9aqpdaWcaaqaamaa dmaabaGaamiraiaacQcacqaHvpGzcaGGQaGaam4yaiaad+gacaWGUb GaamOBaiaadwgacaWGJbGaamiDaiaadMgacaWG2bGaamyAaiaadsha caWG5baacaGLBbGaayzxaaaabaGaamivaiaacQcacqaHvpGzaaaaaa@5D71@     (8)

Equation 8 after simplification will become

[ 4*k*F ] [ Tcceff*ϕ ] = [ D*connectivity ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaWaam WaaeaacaaI0aGaaiOkaiaadUgacaGGQaGaamOraaGaay5waiaaw2fa aaqaamaadmaabaGaamivaiaadogacaWGJbGaamyzaiaadAgacaWGMb GaaiOkaiabew9aMbGaay5waiaaw2faaaaacqGH9aqpdaWcaaqaamaa dmaabaGaamiraiaacQcacaWGJbGaam4Baiaad6gacaWGUbGaamyzai aadogacaWG0bGaamyAaiaadAhacaWGPbGaamiDaiaadMhaaiaawUfa caGLDbaaaeaacaWGubaaaaaa@5885@     (9)

Equation 9 after rearrange will become

k ϕ = [ D*connectivity*Tcceff ] [ 4*k*F ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaalaaabaGaam 4Aaaqaaiabew9aMbaacqGH9aqpdaWcaaqaamaadmaabaGaamiraiaa cQcacaWGJbGaam4Baiaad6gacaWGUbGaamyzaiaadogacaWG0bGaam yAaiaadAhacaWGPbGaamiDaiaadMhacaGGQaGaamivaiaadogacaWG JbGaamyzaiaadAgacaWGMbaacaGLBbGaayzxaaaabaWaamWaaeaaca aI0aGaaiOkaiaadUgacaGGQaGaamOraaGaay5waiaaw2faaaaaaaa@56AA@     (10)

Take the square root of equation 10

k ϕ = [ D*connectivity*Tcceff ] [ 4*k*F ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaakaaabaWaaS aaaeaacaWGRbaabaGaeqy1dygaaaWcbeaakiabg2da9maakaaabaWa aSaaaeaadaWadaqaaiaadseacaGGQaGaam4yaiaad+gacaWGUbGaam OBaiaadwgacaWGJbGaamiDaiaadMgacaWG2bGaamyAaiaadshacaWG 5bGaaiOkaiaadsfacaWGJbGaam4yaiaadwgacaWGMbGaamOzaaGaay 5waiaaw2faaaqaamaadmaabaGaaGinaiaacQcacaWGRbGaaiOkaiaa dAeaaiaawUfacaGLDbaaaaaaleqaaaaa@56EA@     (11)

But; k ϕ =RQI MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadkeacaWG1b GaamiDaiaacUdadaGcaaqaamaalaaabaGaam4Aaaqaaiabew9aMbaa aSqabaGccqGH9aqpcaWGsbGaamyuaiaadMeaaaa@4211@     (12)

Equation 12 will become

RQI= [ D*connectivity*Tcceff ] [ 4*k*F ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadkfacaWGrb Gaamysaiabg2da9maakaaabaWaaSaaaeaadaWadaqaaiaadseacaGG QaGaam4yaiaad+gacaWGUbGaamOBaiaadwgacaWGJbGaamiDaiaadM gacaWG2bGaamyAaiaadshacaWG5bGaaiOkaiaadsfacaWGJbGaam4y aiaadwgacaWGMbGaamOzaaGaay5waiaaw2faaaqaamaadmaabaGaaG inaiaacQcacaWGRbGaaiOkaiaadAeaaiaawUfacaGLDbaaaaaaleqa aaaa@5678@     (13)

Where RQI the reservoir quality index in micrometer. The maximum reservoir quality index can be scaled as

RQImax= [ Dmax*connectivity*Tcceff ] [ 4*k*F ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadkfacaWGrb Gaciysaiaac2gacaWGHbGaamiEaiabg2da9maakaaabaWaaSaaaeaa daWadaqaaiaadseaciGGTbGaaiyyaiaacIhacaGGQaGaam4yaiaad+ gacaWGUbGaamOBaiaadwgacaWGJbGaamiDaiaadMgacaWG2bGaamyA aiaadshacaWG5bGaaiOkaiaadsfacaWGJbGaam4yaiaadwgacaWGMb GaamOzaaGaay5waiaaw2faaaqaamaadmaabaGaaGinaiaacQcacaWG RbGaaiOkaiaadAeaaiaawUfacaGLDbaaaaaaleqaaaaa@5C21@     (14)

The reservoir quality index can be scaled as

Sw= [ RQI RQImax ] 3Df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadofacaWG3b Gaeyypa0ZaamWaaeaadaWcaaqaaiaadkfacaWGrbGaciysaaqaaiaa dkfacaWGrbGaciysaiaac2gacaWGHbGaamiEaaaaaiaawUfacaGLDb aadaahaaWcbeqaaiaaiodacqGHsislcaWGebGaamOzaaaaaaa@485D@     (15)

Where Sw the water saturation; RQI reservoir quality index in micrometer; RQImax maximum reservoir quality index in micrometer; and Df fractal dimension (dimensionless). Insert equation 13 into equation 15

Sw= [ [ D*connectivity*Tcceff ] [ 4*k*F ] RQImax ] 3Df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadofacaWG3b Gaeyypa0ZaamWaaeaadaWcaaqaamaakaaabaWaaSaaaeaadaWadaqa aiaadseacaGGQaGaam4yaiaad+gacaWGUbGaamOBaiaadwgacaWGJb GaamiDaiaadMgacaWG2bGaamyAaiaadshacaWG5bGaaiOkaiaadsfa caWGJbGaam4yaiaadwgacaWGMbGaamOzaaGaay5waiaaw2faaaqaam aadmaabaGaaGinaiaacQcacaWGRbGaaiOkaiaadAeaaiaawUfacaGL DbaaaaaaleqaaaGcbaGaamOuaiaadgfaciGGjbGaaiyBaiaadggaca WG4baaaaGaay5waiaaw2faamaaCaaaleqabaGaaG4maiabgkHiTiaa dseacaWGMbaaaaaa@60B8@     (16)

Insert equation 14 into equation 16

Sw= [ [ D*connectivity*Tcceff ] [ 4*k*F ] [ Dmax*connectivity*Tcceff ] [ 4*k*F ] ] 3Df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadofacaWG3b Gaeyypa0ZaamWaaeaadaWcaaqaamaakaaabaWaaSaaaeaadaWadaqa aiaadseacaGGQaGaam4yaiaad+gacaWGUbGaamOBaiaadwgacaWGJb GaamiDaiaadMgacaWG2bGaamyAaiaadshacaWG5bGaaiOkaiaadsfa caWGJbGaam4yaiaadwgacaWGMbGaamOzaaGaay5waiaaw2faaaqaam aadmaabaGaaGinaiaacQcacaWGRbGaaiOkaiaadAeaaiaawUfacaGL DbaaaaaaleqaaaGcbaWaaOaaaeaadaWcaaqaamaadmaabaGaamirai Gac2gacaGGHbGaaiiEaiaacQcacaWGJbGaam4Baiaad6gacaWGUbGa amyzaiaadogacaWG0bGaamyAaiaadAhacaWGPbGaamiDaiaadMhaca GGQaGaamivaiaadogacaWGJbGaamyzaiaadAgacaWGMbaacaGLBbGa ayzxaaaabaWaamWaaeaacaaI0aGaaiOkaiaadUgacaGGQaGaamOraa Gaay5waiaaw2faaaaaaSqabaaaaaGccaGLBbGaayzxaaWaaWbaaSqa beaacaaIZaGaeyOeI0IaamiraiaadAgaaaaaaa@7913@     (17)

Equation 17 after simplification will become

Sw= [ D Dmax ] 3D = [ D 1 2 D max 1 2 ] 3Df MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadofacaWG3b Gaeyypa0ZaamWaaeaadaWcaaqaamaakaaabaGaamiraaWcbeaaaOqa amaakaaabaGaamiraiGac2gacaGGHbGaaiiEaaWcbeaaaaaakiaawU facaGLDbaadaahaaWcbeqaaiaaiodacqGHsislcaWGebaaaOGaeyyp a0ZaamWaaeaadaWcaaqaaiaadseadaahaaWcbeqaamaalaaabaGaaG ymaaqaaiaaikdaaaaaaaGcbaGaamiraiGac2gacaGGHbGaaiiEamaa CaaaleqabaWaaSaaaeaacaaIXaaabaGaaGOmaaaaaaaaaaGccaGLBb GaayzxaaWaaWbaaSqabeaacaaIZaGaeyOeI0IaamiraiaadAgaaaaa aa@52D5@     (18)

Equation 18 is the proof of equation 1 which relates water saturation; diffusion coefficient; maximum diffusion coefficient; and fractal dimension. The capillary pressure can be scaled as

Log Sw=[ Df3 ]*logPc+constant MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadYeacaWGVb Gaam4zaabaaaaaaaaapeGaaiiOa8aacaWGtbGaam4Daiabg2da9maa dmaabaGaamiraiaadAgacqGHsislcaaIZaaacaGLBbGaayzxaaGaai OkaiaacYgacaGGVbGaai4zaiaaccfacaGGJbGaey4kaSIaai4yaiaa c+gacaGGUbGaai4CaiaacshacaGGHbGaaiOBaiaacshaaaa@51FB@     (19)

Where Sw the water saturation, Pc the capillary pressure, and Df the fractal dimension

Result and discussion

Based on field observation the Shajara Reservoirs of the Shajara Formation of the Permo-Carboniferous Unayzah Group were divided here into three units as described in Figure 1. These units from bottom to top are: Lower Shajara Reservoir, Middle Shajara reservoir, and Upper Shajara Reservoir.

Figure 1 Surface type section of the Shajara Reservoirs of the Permo – Carboniferous Shajara Formation at latitude 26˚ 52' 17.4", longitude 43˚ 36' 18".

Their acquired results of the diffusion coefficient fractal dimension and capillary pressure fractal dimension are displayed in Table 1. Based on the attained results it was found that the diffusion coefficient fractal dimension is equal to the capillary pressure fractal dimension. The maximum value of the fractal dimension was found to be 2.7872 assigned to sample SJ13 from the Upper Shajara Reservoir as verified in Table 1. Whereas the minimum value of the fractal dimension 2.4379 was reported from sample SJ3 from the Lower Shajara reservoir as displayed in Table 1. The diffusion coefficient fractal dimension and capillary pressure fractal dimension were observed to increase with increasing permeability as proofed in Table 1 owing to the possibility of having interconnected channels.

Formation

Reservoir

Sample

Porosity %

k md

Positive slope of the first procedure Slope=3-Df

Negative slope of the second procedure Slope=Df-3

Diffusion coefficient fractal dimension

Capillary pressure fractal dimension

Permo-Carboniferous Shajara Formation

Upper Shajara Reservoir

SJ13

25

973

0.2128

-0.2128

2.7872

2.7872

SJ12

28

1440

0.2141

-0.2141

2.7859

2.7859

SJ11

36

1197

0.2414

-0.2414

2.7586

2.7586

Middle Shajara Reservoir

SJ9

31

1394

0.2214

-0.2214

2.7786

2.7786

SJ8

32

1344

0.2248

-0.2248

2.7752

2.7752

SJ7

35

1472

0.2317

-0.2317

2.7683

2.7683

Lower Shajara Reservoir

SJ4

30

176

0.3157

-0.3157

2.6843

2.6843

SJ3

34

56

0.5621

-0.5621

2.4379

2.4379

SJ2

35

1955

0.2252

-0.2252

2.7748

2.7748

SJ1

29

1680

0.2141

-0.2141

2.7859

2.7859

Table 1 Petrophysical model showing the three Shajara Reservoir Units with their corresponding values of diffusion coefficient fractal dimension and capillary pressure fractal dimension

The Lower Shajara reservoir was denoted by six sandstone samples (Figure 1), four of which label as SJ1, SJ2, SJ3 and SJ4 were selected for capillary pressure measurement as confirmed in Table 1. Their positive slopes of the first procedure (log of the ratio of diffusion coefficient to maximum diffusion coefficient versus log wetting phase saturation) and negative slopes of the second procedure (log capillary pressure versus log wetting phase saturation (log Sw) is delineated in Figures 2–5. Their diffusion coefficient fractal dimension and capillary pressure fractal dimension values are shown in Table 1. As we proceed from sample SJ2 to SJ3 a pronounced reduction in permeability due to compaction was reported from 1955md to 56md which reflects decrease in diffusion coefficient fractal dimension from 2.7748 to 2.4379 as specified in Table 1. Again, an increase in grain size and permeability was verified from sample SJ4 whose diffusion coefficient fractal dimension and capillary pressure fractal dimension was found to be 2.6843 as described in Table 1.

Figure 2 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ1.

Figure 3 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ2.

Figure 4 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ3.

Figure 5 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ4.

In contrast, the Middle Shajara reservoir which is separated from the Lower Shajara reservoir by an unconformity surface as shown in Figure 1. It was designated by four samples (Figure 1), three of which named as SJ7, SJ8, and SJ9 as illustrated in Table 1 were selected for capillary measurements as described in Table 1. Their positive slopes of the first procedure and negative slopes of the second procedure are shown in Figures 6–8. Additionally, their diffusion coefficient fractal dimensions and capillary pressure fractal dimensions show similarities as delineated in Table 1.Their fractal dimensions are higher than those of samples SJ3 and SJ4 from the Lower Shajara Reservoir due to an increase in their permeability as explained in Table 1.

Figure 6 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ7.

Figure 7 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ8.

Figure 8 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ9.

On the other hand, the Upper Shajara reservoir separated from the Middle Shajara reservoir by yellow green mudstone as revealed in Figure 1. It is defined by three samples so called SJ11, SJ12, SJ13 as explained in Table 1. Their positive slopes of the first procedure and negative slopes of the second procedure are displayed in Figures 9–11. Moreover, their diffusion coefficient fractal dimension and capillary pressure fractal dimension are also higher than those of sample SJ3 and SJ4 from the Lower Shajara Reservoir due to an increase in their permeability as clarified in Table 1.

Figure 9 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ11.

Figure 10 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ12.

Figure 11 Log (D1/2/Dmax1/2) & log Pc versus log Sw for sample SJ13.

Overall a plot of diffusion coefficient fractal dimension versus capillary pressure fractal dimension as shown in Figure 12 reveals three permeable zones of varying Petrophysical properties. Such variation in fractal dimension can account for heterogeneity which is a key parameter in reservoir quality assessment. This reservoir heterogeneity was also confirmed by plotting positive slope of the first procedure versus negative slope of the second procedure as described in Figure 13.

Figure 12 Diffusion coefficient fractal dimension versus capillary pressure fractal dimension.

Figure 13 Slope of the first procedure versus slope of the second procedure.

Conclusion

The sandstones of the Shajara Reservoirs of the Permo-Carboniferous Shajara Formation were divided here into three units based on diffusion coefficient fractal. The Units from bottom to top are Lower Shajara Diffusion Coefficient Fractal dimension unit, Middle Shajara Diffusion Coefficient Fractal dimension unit, and Upper Shajara Diffusion Coefficient Fractal dimension unit. These units were also confirmed by capillary pressure fractal dimension. The heterogeneity increases with increasing permeability, increasing fractal dimension, decreasing compaction owing to possibility of having interconnected channels.

Acknowledgments

The author would like to thank College of Engineering, King Saud University, Department of Petroleum and Natural Gas Engineering, Department of Chemical Engineering, Research Centre at College of Engineering and King Abdulla Institute for Research and Consulting Studies for their supports.

Conflicts of interest

The authors declared there is no conflict of interest.

References

  1. Frenkel J. On the theory of seismic and seismoelectric phenomena in a moist soil. Journal of Physics. 1944;3:230–241.
  2. Li K, Williams W. Determination of capillary pressure function from resistivity data. Transp Porous Media. 2007;67(1):1–15.
  3. Revil A, Jardani A. Seismoelectric response of heavy oil reservoirs: theory and numerical modelling. Geophys J Int. 2010;180(2):781–797.
  4. Dukhin A, Goetz P, Thommes M. Seismoelectric effect: a non-isochoric streaming current. 1 experiment. J Colloid Interface Sci. 2010;345(2):547–553.
  5. Guan W, Hu J, Wang Z. Permeability inversion from low-frequency seismoelectric logs in fluid-saturated porous formations. Geophys Prospect. 2012;61(1):120–133.
  6. Hu H, Guan W, Zhao W. Theoretical studies of permeability inversion from seismoelectric logs. Geophysical Research Abstracts. 2012.
  7. Borde C, S´en´echal P, Barri`ere J, et al. Impact of water saturation on seismoel+ectric transfer functions: a laboratory study of co-seismic phenomenon. Geophys J Int. 2015;200(3):1317–1335.
  8. Jardani A, Revi A. Seismoelectric couplings in a poroelastic material containing two immiscible fluid phases. Geophys J Int. 2015;202(2):850–870.
  9. Holzhauer J, Brito D, Bordes C, et al. Experimental quantification of the seismoelectric transfer function and its dependence on conductivity and saturation in loose sand. Geophysics Prospect. 2016;65(4):1097–1120.
  10. Rong P, Wei J-X, Di B-R, et al. Experimental research on seismoelectric effects in sandstone. Applied Geophysics. 2016;13(3):425–436.
  11. Djuraev U, Jufar S R, Vasant P. Numerical Study of frequency-dependent seismoelectric coupling in partially-saturated porous media. MATEC Web of Conferences. 2017.
  12. Alkhidir KEME. Pressure head fractal dimension for characterizing Shajara Reservoirs of the Shajara Formation of the Permo-Casrboniferous Unayzah Group, Saudi Arabia. Arch Pet Environ Biotech. 2017;(2):1–7.
  13. Alkhidir KEME. Geometric relaxation time of induced polarization fractal dimension for characterizing Shajara Reservoirs of the Shajara Formation of the Permo-Carboniferous Unayzah Group, Saudi Arabia. Scifed J Petroleum. 2018;2(1):1–6.
  14. Alkhidir KEME. Geometric relaxation time of induced polarization fractal dimensionfor characterizing Shajara Reservoirs of the Shajara formation of the Permo-Carboniferous Unayzah Group-Permo. Int J Pet Res. 2018;2(1):105–108.
  15. Alkhidir KEME. Arithmetic relaxation time of induced polarization fractal dimension for characterizing Shajara Reservoirs of the Shajara Formation. Nanosci and Nanotechnol. 2018;2(1):1–8.
Creative Commons Attribution License

©2019 Alkhidir. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.