/[escript]/trunk/doc/inversion/CookGravity.tex
ViewVC logotype

Diff of /trunk/doc/inversion/CookGravity.tex

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4155 by azadeh, Tue Jan 22 00:23:39 2013 UTC revision 4156 by azadeh, Wed Jan 23 03:45:14 2013 UTC
# Line 2  Line 2 
2    
3  This instrument is designed to prepare practical and obvious information to apply gravity and magnetic inversion packages which have joint application in magnetic and gravity inversion. This codes were written in 2012 at the department of Earth science , University of Queensland.\\  This instrument is designed to prepare practical and obvious information to apply gravity and magnetic inversion packages which have joint application in magnetic and gravity inversion. This codes were written in 2012 at the department of Earth science , University of Queensland.\\
4  It enables geologists and geophysicists, but who is not necessarily versed in the details of inverse theory, to process, visualize and interpret multi-volume geophysical data using attributes and modern visualization techniques.This is not only a functional interpretation system, it is also a research and development environment for geophysics analysis.\\  It enables geologists and geophysicists, but who is not necessarily versed in the details of inverse theory, to process, visualize and interpret multi-volume geophysical data using attributes and modern visualization techniques.This is not only a functional interpretation system, it is also a research and development environment for geophysics analysis.\\
 We present a package for inverting ground or airborn surveys gravity and magnetic data to generate a 2-D or 3-D distribution of density and susceptibility contrast. In this approach, the earth is clearly modeled by using a large number of rectangular cells of constant value such as density or susceptibility, and the final distribution is obtained by minimizing a model objective function subject to fitting the observed data.\\  
5    
6  This package could be used to provide a model which introduce density and susceptibility together to fit a given set of magnetic and gravity anomalies. The given data might contain negative and positive in value for both onshore and offshore region. However the model is a 3-dimentional that any cross or depth sections are extracted from. Being the result of big region as well as small area is one of the advantages of this package.\\  We present a package for inverting ground or airborn surveys gravity and magnetic data to generate a 2-D or 3-D distribution of density and/or susceptibility contrast. In this approach, the earth is clearly modeled by using a large number of rectangular cells of constant value such as density or susceptibility, and the final distribution is obtained by minimizing a model objective function subject to fitting the observed data.\\
7    
8    This package could be used to provide an outstanding model which introduce density and susceptibility together to fit a given set of magnetic and gravity anomalies. The given data might contain negative and positive in value for both onshore and offshore region. However the model is a 3-dimentional that any cross or depth sections are extracted from. Being the result of big region as well as small area is one of the advantages of this package.\\
9  Data type should be prepared attentively to have all corrections and processing. Model is described as a big volume which has changed in density and susceptibility smoothly in all directions. The apparent characteristics of the topography are very sophisticated.\\  Data type should be prepared attentively to have all corrections and processing. Model is described as a big volume which has changed in density and susceptibility smoothly in all directions. The apparent characteristics of the topography are very sophisticated.\\
10  \newpage  \newpage
11    
12  \begin{figure}  \begin{figure}
13  \centering  \centering
14  \includegraphics[width=\textwidth]{pasted1.png}  \includegraphics[width=\textwidth]{density1.png}
15  \caption{Depth image trough a 3 dimensional gravity inversion which presents discrepancy in density. Its surrounded area is before the padding with uncertainty function. The gravity of padding area of the  \caption{Contour image through a 3 dimensional gravity inversion which presents discrepancy in density. Its surrounded area is before the padding with uncertainty function. The gravity of padding area of the
16  model is not defined. Increasing and decreasing in density are indicated  model is not defined. Increasing and decreasing in density are indicated with red color and blue color respectively.}
17  with red color and blue color respectively.}  \end{figure}
18    
19    \begin{figure}
20    \centering
21    \includegraphics[width=\textwidth]{density2.png}
22    \caption{Depth image across previous 3D gravity inversion which presents discrepancy in density. Diversity in density are detected with colors and contours.}
23  \end{figure}  \end{figure}
24    
25  Such inversion calculations invariably depend on assumptions about the characteristic of the forward modeling. From this point of view, we can say that traditional inversion is fundamentally based on a model of the subsurface which have used as a prior knowledge.\\  Such inversion calculations invariably depend on assumptions about the characteristic of the forward modeling. From this point of view, we can say that traditional inversion is fundamentally based on a model of the subsurface which have used as a prior knowledge.\\
# Line 28  The amount of $g$ at sea level is about Line 34  The amount of $g$ at sea level is about
34  The Gravity data are taken from onshore and offshore observation recorded at many gravity stations with high precision of determination in elevation and position (latitude and longitude). All Raw reading gravity observations require to process with many corrections.\\  The Gravity data are taken from onshore and offshore observation recorded at many gravity stations with high precision of determination in elevation and position (latitude and longitude). All Raw reading gravity observations require to process with many corrections.\\
35    
36  The sun and moon gravitational forces make curvature in Earth's shape. These tide effects change figure of oceans, atmosphere and even solid body of the Earth, which impress gravity measurement and it is necessary to compensate it, which vary with location, date and time of the day.\\  The sun and moon gravitational forces make curvature in Earth's shape. These tide effects change figure of oceans, atmosphere and even solid body of the Earth, which impress gravity measurement and it is necessary to compensate it, which vary with location, date and time of the day.\\
37    
38  The surface of the Earth is lumpy on land and water. However for Geophysical and Geological study, a smooth closed surface is assumed. The main one is a spheroid flattened at the poles which is called ellipsoid. The new data are used to defined a best-fitting obtained ellipsoid. The second suggestion is geoid which is really mathematical convenience. There is a uniform mass between gravity stations and ellipsoid, that's effect must be removed with corrections.\\  The surface of the Earth is lumpy on land and water. However for Geophysical and Geological study, a smooth closed surface is assumed. The main one is a spheroid flattened at the poles which is called ellipsoid. The new data are used to defined a best-fitting obtained ellipsoid. The second suggestion is geoid which is really mathematical convenience. There is a uniform mass between gravity stations and ellipsoid, that's effect must be removed with corrections.\\
39    
40  Also the level of topography for hilly and valley measurements is important. The gravity amount which is made up by that equal the mass of hill or valley must be added as a terrian correction to have a measurement on a level surface.\\  Also the level of topography for hilly and valley measurements is important. The gravity amount which is made up by that equal the mass of hill or valley must be added as a terrian correction to have a measurement on a level surface.\\
41  Because gravity descend towards the poles the latitude correction must be added to the observed gravity.\\  Because gravity descend towards the poles the latitude correction must be added to the observed gravity.\\
42  Free-air correction that must be added to observation, ignores the effects of material between the measurement and refrence level which is positive for above sea-level station and is negetive for station below sea-leve.\\  Free-air correction that must be added to observation, ignores the effects of material between the measurement and refrence level which is positive for above sea-level station and is negetive for station below sea-leve.\\
# Line 41  For this adjustment first of all Tidal c Line 49  For this adjustment first of all Tidal c
49  \textbf{Input File} \\  \textbf{Input File} \\
50    
51  For starting up the inversion 2 files are needed. Each of the two files contains a series of parameters which must appear in the correct order, as described in the next paragraph. Each parameter is marked by a keyword, which is followed either on the same line or the next line by one or more parameter values. \\  For starting up the inversion 2 files are needed. Each of the two files contains a series of parameters which must appear in the correct order, as described in the next paragraph. Each parameter is marked by a keyword, which is followed either on the same line or the next line by one or more parameter values. \\
52  The first file is contained the gravity anomaly, including number of points, the accurate location (latitude and longitude) of the observed position and the value of the anomaly after all gravity corrections.\\  The first file is contained the gravity anomaly, including number of points, the accurate location (latitude and longitude) of the observed position and the value of the anomaly after all gravity corrections. The main part in collecting the data is spacing of the observation point which means that all data must be gridded and its spacing affected the simulation resolution wisely.\\
53  The other is run_gravity with py extension which is related to the codes though it needs some constraints to have a good results in inversion.\\  The other is run_gravity with py extension which is related to the codes though it needs some constraints to have a good results in inversion.\\
54    
55  A small part of sample of run_gravity:\\  A small part of sample of run_gravity:\\
56  \begin{verbatim}  \begin{verbatim}
57  mu=100  mu=100
 n_cells_in_data=100  
 depth_offset=0.*U.km  
 l_data = 100 * U.km  
 l_pad=40*U.km  
58  THICKNESS=20.*U.km  THICKNESS=20.*U.km
59  l_air=6*U.km  DATASET='QLD_west.nc'
60    PAD_X = 0.2
61    PAD_Y = 0.2
62    l_air = 6. * U.km
63    n_cells_v = 25
64  \end{verbatim}  \end{verbatim}
65    
66  Run_gravity file consist many options to implement which control how inversion is performed such as padding area, depth , MU factor,\ldots.  Run_gravity file consist many options to implement which control how inversion is performed such as padding area, depth , MU factor,\ldots.
67    
68  \begin{description}      \begin{description}    
69  \item[MAX\_ITER]  \item[MAX\_ITER]
70  Specifying maximum iteration depends on model and the area which have been selected to have an inversion in it also your hard capacity. However the best result were built with 200 iterations. In addition all steps of inversion could be traced and the suited one selected.  Specifying maximum iteration depends on model and the area which have been selected to have an inversion in it and your hard capacity. However the best result were built with 200 iterations. In addition all steps of inversion could be traced and the suited one selected.
71    
72  \item[l_data or PAD\_X, PAD\_Y] To implement the bound constraints in this file, padding area in $x$ and $y$ direction should be determined. In a rectangular area, same padding for both direction is preferable. If directional area will be fixed with elongation in one orientation it does not matter to change the padding area. In addition for 2D inversion padding just add in one direction.  \item[l_pad or PAD\_X, PAD\_Y] To implement the boundry constraints in this file, padding area in $x$ and $y$ direction should be determined. In a rectangular area, same padding for both direction is preferable. If directional area will be fixed with elongation in one orientation it does not matter to change the padding area. In addition for 2D inversion padding just add in one direction.
73    
74  \item[DATASET] In run_gravity file just the location and the name of the data file have to be fixed the source of data and run_gravity must be in a folder.  \item[DATASET] In run_gravity file just the location and the name of the data file have to be fixed. Both source of data and run_gravity fles must be in a folder.
75    
76  \item[THICKNESS] Depth of the model should be assigned to have dipper or shallower inversion also it is assigned to the layer where shows inversion in.The parameter values must be real numbers, and they represent depths in km.  \item[THICKNESS] Depth of the model should be assigned to have dipper or shallower inversion also it is assigned to the layer where shows inversion in.The parameter values must be real numbers, and they represent depths in km.
77    
78  \item[l_air] Length of air is hight of the model above the sea level.  \item[l_air] Length of air is hight of the model above the sea level.
79    
80  \item[n_cells_in_data] The last important part of the inversion property is the number of elements in data of the model which shows the finer or coarser cells in the model so its delimitation have affected on resolution.  \item[n_cells_v] The last important part of the inversion property is the number of elements in data of the model which shows the finer or coarser cells in the model so its delimitation have affected on resolution of the inversion.
81    
82  \item[mu]It is defined in accordance with the noise of data and it has a wide range to select from 0.0001 to 100.  \item[mu]It is defined in accordance with the noise of data and it has a wide range to select from 0.0001 to 100 for each inversion.
83    
84  \end{description}  \end{description}
85    
# Line 79  Specifying maximum iteration depends on Line 87  Specifying maximum iteration depends on
87    
88  At the end of each inversion iteration, package produces a new inversion file with 'gravin.silo' name which is stored in the file. This silo file shows the inversion result which does not have the number of iterations stage.\\  At the end of each inversion iteration, package produces a new inversion file with 'gravin.silo' name which is stored in the file. This silo file shows the inversion result which does not have the number of iterations stage.\\
89   In terminal indicates the specifications of inversion iterations. In this descriptions of the paths of all inversion during stages of modeling are cleared. The format of the file will be described here as it is designed to be used directly for analysis or debugging.\\   In terminal indicates the specifications of inversion iterations. In this descriptions of the paths of all inversion during stages of modeling are cleared. The format of the file will be described here as it is designed to be used directly for analysis or debugging.\\
90  After final iteration the silo file is visible with some software which shows silo file format. It illuminate density distribution in the area (not in the padding) which create the gravity input data.\\  After final iteration the silo file, the result of inversion, is visible with some software which shows silo file format. It illuminate density distribution in the area (not in the padding) which create the gravity input data.\\
91    
92  \textbf{Reference}\\  \textbf{Reference}\\
93    
94  There are three examples for 2D and 3D gravity inversions with artificial input data.\\  There are three examples for 2D and 3D gravity inversions with artificial input data.\\
95  In first step, an area with synthetic density section is suggested. Then based on forward modeling it's gravitational data is collected. Afterwards with generated gravity data, escripts simulate a volume of inverted density. Whilst new density mass could be compared with the synthetic density section to verify the inversion.\\  In first step, an area with synthetic density section was suggested. Then based on forward modeling its gravitational data was collected. Afterwards with generated gravity data, escripts simulated a volume of inverted density. Whilst new density mass could be compared with the synthetic density section to verify the inversion.\\
96    
97  Some of the presumptions are the same for all of the examples to simplify the situation to make a logical comparison between synthetic input and output. which is as followed:\\  Some of the presumptions were the same for all of the examples to simplify the situation to make a logical comparison between synthetic input and output. which is as followed:\\
98  \begin{verbatim}  \begin{verbatim}
99  n_cells_in_data=100  n_cells_in_data=100
100  depth_offset=0.*U.km  depth_offset=0.*U.km
# Line 97  l_air=6*U.km Line 105  l_air=6*U.km
105  \end{verbatim}  \end{verbatim}
106  The others assumptions comes with each examples.\\  The others assumptions comes with each examples.\\
107    
108  1.  A 2D density section with a maximum in centre was assumed. The refrence density and inverted will be shown. The padding area is excluded. \\  1.  A 2D density section with a maximum in centre was assumed. The refrence density and inverted will be shown. The padding area is excluded. (\ref{fig:gravity2D1})\\
109  \begin{verbatim}  \begin{verbatim}
110  n_cells_in_data=100  n_cells_in_data=100
111  n_humbs_h= 1  n_humbs_h= 1
# Line 108  mu=100 Line 116  mu=100
116  \begin{figure}  \begin{figure}
117  \centering  \centering
118  \includegraphics[width=\textwidth]{grav2D1.png}  \includegraphics[width=\textwidth]{grav2D1.png}
119  \caption{2D density model up) reference  down) result}  \caption{2D density model up) reference    down) result}
120    \label{fig:gravity2D1}
121  \end{figure}  \end{figure}
122    
123  2. A 2D density properties with two maximun in corners and one minimum in the centre was inverted. The result have eliminated the effects in padding area. \\  
124    2. A 2D density properties with two maximun in corners and one minimum in the centre was inverted. The result have eliminated the effects in padding area.(\ref{fig:gravity2D3}) \\
125  \begin{verbatim}  \begin{verbatim}
126  n_cells_in_data=100  n_cells_in_data=100
127  n_humbs_h= 3  n_humbs_h= 3
# Line 124  mu=100 Line 133  mu=100
133  \centering  \centering
134  \includegraphics[width=\textwidth]{grav2D3.png}  \includegraphics[width=\textwidth]{grav2D3.png}
135  \caption{2D density model up) reference  down) result}  \caption{2D density model up) reference  down) result}
136    \label{fig:gravity2D3}
137  \end{figure}  \end{figure}
138    
139  3. A 3D model with a maximum in the center was used as input data and the result after simulation in shown in the next image which determined a good distribution of density through the model in the main area.\\  3. A 3D model with a maximum in the center was used as input data and the result after simulation in shown in the next image which determined a good distribution of density through the model in the main area.(\ref{fig:gravity3D} and \ref{fig:gravity3D1})\\
140  \begin{verbatim}  \begin{verbatim}
141  n_cells_in_data=50  n_cells_in_data=50
142  n_humbs_h= 1  n_humbs_h= 1
# Line 139  mu=10 Line 148  mu=10
148  \centering  \centering
149  \includegraphics[width=\textwidth]{density3D-ref.png}  \includegraphics[width=\textwidth]{density3D-ref.png}
150  \caption{3D density model of reference as synthetic data}  \caption{3D density model of reference as synthetic data}
151    \label{fig:gravity3D}
152  \end{figure}  \end{figure}
153    
154  \begin{figure}  \begin{figure}
155  \centering  \centering
156  \includegraphics[width=\textwidth]{gravity3D.png}  \includegraphics[width=\textwidth]{gravity3D.png}
157  \caption{3D density model of result}  \caption{3D density model of result}
158    \label{fig:gravity3D1}
159  \end{figure}  \end{figure}

Legend:
Removed from v.4155  
changed lines
  Added in v.4156

  ViewVC Help
Powered by ViewVC 1.1.26