Scribd
Upload
141 views3 pages

Gravity Data Processing

The document describes the process of gravity data processing and inversion that the author performed on raw gravity data. This included correcting the data for drift, applying gravity corrections to calculate the Bouguer anomaly using Matlab software, constructing a contour map using Surfer, extracting 2D profile data, performing regional-residual separation to filter near surface anomalies, and applying a 2D gravity inversion algorithm based on a research paper to obtain subsurface mass distribution.

Uploaded by

Santosh Kumar
Copyright
© Public Domain
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
141 views3 pages

Gravity Data Processing

The document describes the process of gravity data processing and inversion that the author performed on raw gravity data. This included correcting the data for drift, applying gravity corrections to calculate the Bouguer anomaly using Matlab software, constructing a contour map using Surfer, extracting 2D profile data, performing regional-residual separation to filter near surface anomalies, and applying a 2D gravity inversion algorithm based on a research paper to obtain subsurface mass distribution.

Uploaded by

Santosh Kumar
Copyright
© Public Domain
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Work Done!

I was given raw gravity data and using that I did gravity data processing and inversion. First of all I corrected the data
for ​Drift changes​ due to instrument error then applied all the gravity corrections to get ​Bouguer anomaly​ caused by
anomalous density distribution using Matlab software.

%%===============Instrument Constant==============>
IC = 0.1008;

%%==========Base station absolute reading===========>


absolute_base_station_reading = 978796.93;

%=============Load Data In Matlab==============>%


data = xlsread('data.xlsx');

%%=============Separate all the datas-============>


latitude = data(:,1);
longitude = data(:,2);
elevation = data(:,3);
time = data(:,4);
reading = data(:,5);
del_time = time - time(1);

%%============Drift Correction==============>
len = length(reading);
format long;
drift_rate = (reading(len) - reading(1))/del_time(len);
total_drift = drift_rate * del_time;
corrected_reading = reading - total_drift;
division = corrected_reading - corrected_reading(1);
difference_mGal = division * IC;
G_observed = difference_mGal + absolute_base_station_reading;

%%===========Theoretical gravity calculation ============>


phi = (latitude * pi)/180;
alpha = 0.0053204;
bita = 0.0000058;
G_theoretical = 978031.8*( (1) + (alpha*(sin(phi).^2)) + (bita*(sin(2*phi).^2)) );

%%==========Latitude correction calculation================>


del_angle = (latitude - latitude(1));
surface_distance = 111 * del_angle ;
G_lattitude = 0.812 * sin(2*phi(1)) * surface_distance;

%%===========Free air correction=================>


G_air = (-1) * 0.308595801 * elevation;

%%========Bouguer plate correction=================>


density = 2.67; %in gm/cc
G_bouguer = 0.04188 * density * elevation;

%=========Bouguer anomaly calculation ==============>


format short;
bouguer_anomaly = G_observed + G_lattitude - G_air - G_bouguer - G_theoretical;
After getting ​Bouguer anomaly​ I constructed contour map using Surfer software. Then extracted the data from
contour map for 2D profile which I used to do 2D gravity inversion.
Regional-residual separation​ is done to remove the effect of ​deep anomalous body​ and filter out the ​near surface
anomaly. ​After that I used ​GRAVMAG​ software to visualize the the possible density distribution in the subsurface

After this I applied 2D gravity data inversion algorithm based on tomography problem to get actual subsurface mass
distribution.

** For 2D inversion. I took reference from of research paper published by Animesh Mandal, research scholar at IIT
Kharagpur.

Link: ​Integrated Geophysical Survey for the Delineation of Sub-surface Structure Associated with Uranium 
Mineralization, South Purulia Shear Zone, India 

You might also like