Skip to content

IFOS2D
IFOS2D
Source
Select Git revision
...
Target
Select Git revision
 
Commits (21)
Expand all Hide whitespace changes
Inline Side-by-side
Showing with 5142 additions and 5050 deletions
README.md
View file @ 04ca7c45
......@@ -9,10 +9,11 @@ IFOS2D is the reverse (inverse) of our 2-D Finite-Difference forward solver [**S
The [**manual**](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/wikis/home) is included in the download archive or can be downloaded [here](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/wikis/home).
# Download and Newsletter
# Download, Newsletter and Contact
Release: [![build status](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/badges/master/build.svg)](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/commits/master) Beta: [![build status](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/badges/develop/build.svg)](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/commits/develop)
You can download the [**latest Release 2.0.3**](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/tags/Release_2.0.3) or the current [**Beta-Version**](https://git.scc.kit.edu/GPIAG-Software/IFOS2D/tree/develop).
To receive news and updates please [register](http://www.gpi.kit.edu/Software-FWI.php) on the email list [IFOS@lists.kit.edu](http://www.gpi.kit.edu/Software-FWI.php).
Please use this list also to ask questions or to report problems or bugs.
To receive news and updates please [register](http://www.gpi.kit.edu/Software-FWI.php) on the e-mail list [IFOS@lists.kit.edu](http://www.gpi.kit.edu/Software-FWI.php).
Please use this list also to ask questions or to report problems or bugs!
\ No newline at end of file
doc/latex/6_Parameter_Definition.tex
View file @ 04ca7c45
......@@ -452,7 +452,7 @@ This section covers some general inversion parameters. The maximum number of ite
If models are read from binary files appropriate file extensions are required for the different models (see section \ref{gen_of_mod}). Depending on the data different components of the seismic sections can be back propagated. For two component data (x- and y-component) set ADJOINT\_TYPE=1, only the y-component (ADJOINT\_TYPE=2) and only the x-component (ADJOINT\_TYPE=3). For the inversion of pressure seismograms ADJOINT\_TYPE=4 has to be used.
During the inversion the misfit values are saved in a log file specified in MISFIT\_LOG\_FILE. The log file consists of eight or nine columns and each line corresponds to one iteration step. The used step length is written in the first column. In the second to fourth column the three test step lengths used for the step length estimation are saved. The corresponding misfit values for these test step lengthes and the test shots are written to column five to seven. Column eight corresponds to the total misfit for all shots and if you use frequency filtering then the ninth column corresponds to the corner frequency of the lowpass filter used in the inversion step.
During the inversion the misfit values are saved in a log file specified in MISFIT\_LOG\_FILE. The log file consists of nine or ten columns and each line corresponds to one iteration step. The used step length is written in the first column. In the second to fourth column the three test step lengths used for the step length estimation are saved. The corresponding misfit values for these test step lengthes and the test shots are written to column five to seven. Column eight corresponds to the total misfit for all shots and if you use frequency filtering then the ninth column corresponds to the corner frequency of the lowpass filter used in the inversion step. If TIMEWIN=1 then the tenth column shows the GAMMA value that was used at the current iteration.
In general IFOS2D tries to minimize the misfit in the particle displacement between the observed data and the synthetic data. If you set the switch VELOCITY to 1 the misfit in the particle velocity seismograms is minimized.
......@@ -489,12 +489,12 @@ Default values are:
With the use of a workflow file, you can define different FWI stages. For instance, one FWI stage could refer to one corner frequency of a low-pass filter or to a specific time window. Every line in the workflow file corresponds to one FWI stage. To use a workflow file the switch USE\_WORKFLOW have to be set to 1. The algorithm will automatically change to the next line of the workflow file if the abort criterium of the current line is reached or if no step length could be found which reduces the misfit. The structure of the variables inside the workflow file is as follow: \\
\noindent
\begin{tabular}{lllllllllllll}
\# & INV\_VS & INV\_VP & INV\_RHO & PRO & TIME\_FILT & F\_HIGH\_PASS & F\_LOW\_PASS & WAVETYPE & 0 & 0 & EPRECOND & EPSILON\_WE\\
\begin{tabular}{llllllllllIlll}
\# & INV\_VS & INV\_VP & INV\_RHO & PRO & TIME\_FILT & F\_HIGH\_PASS & F\_LOW\_PASS & WAVETYPE & 0 & 0 & EPRECOND & EPSILON\_WE & GAMMA \\
\end{tabular}
\ \\
The first column is the number of the line. With INV\_* etc. you can activate the inversion for VS, VP or RHO, respectively. The abort criterium in percent for this FWI stage will be the declared in the variable PRO. With TIME\_FILT you can activate the frequency filtering with the corner frequency F\_HIGH\_PASS of the high-pass and F\_LOW\_PASS of the low-pass filter. WAVETYPE can be used to switch between PSV and SH modeling, however it is recommended to use PSV modeling (WAVETYPE==1) due to the current development on SH modeling. The following to zeros are placeholders for an upcoming update. With EPRECOND and EPSILON\_WE you can control the approx. Hessian. Please note, that all features which are used eg. TIME\_FILT (see section \ref{sec:filtering}) within the workflow have to be activated in the .JSON file.
The first column is the number of the line. With INV\_* etc. you can activate the inversion for VS, VP or RHO, respectively. The abort criterium in percent for this FWI stage will be the declared in the variable PRO. With TIME\_FILT you can activate the frequency filtering with the corner frequency F\_HIGH\_PASS of the high-pass and F\_LOW\_PASS of the low-pass filter. WAVETYPE can be used to switch between PSV and SH modeling, however it is recommended to use PSV modeling (WAVETYPE==1) due to the current development on SH modeling. The following two zeros are placeholders for an upcoming update. With EPRECOND and EPSILON\_WE you can control the approx. Hessian. Please note, that all features which are used eg. With GAMMA you can define the exponential taper applied to both synthetic and observed data. TIME\_FILT (see section \ref{sec:filtering}) within the workflow have to be activated in the .JSON file.
For an example of a workflow file, have a look in the par/ folder.
......@@ -731,7 +731,7 @@ If you are using frequeny filtering (TIME\_FILT==1) during the inversion, you ca
"PICKS_FILE" : "./picked_times/picks"
"TWLENGTH_PLUS" : "0.01",
"TWLENGTH_MINUS" : "0.01",
"GAMMA" : "100000",
"GAMMA" : "30",
\end{verbatim}}}
{\color{red}{\begin{verbatim}
......@@ -741,7 +741,7 @@ Default values are:
To apply time windowing in a time series the parameter TIMEWIN must set to 1. A automatic picker routine is not integrated. The point in time (picked time) for each source must be specified in separate files. The folder and file name can be set with the parameter PICKS\_FILE. The files must be named like this PICKS\_FILE\_<sourcenumber>.dat. So the number of sources must be equal to the number of files. Each file must contain the picked times for every receiver.
The parameters TWLENGTH\_PLUS and TWLENGTH\_MINUS specify the length of the time window after (PLUS) and before (MINUS) the picked time. The unit is seconds. The damping factor GAMMA must be set individually.
The parameters TWLENGTH\_PLUS and TWLENGTH\_MINUS specify the length of the time window after (PLUS) and before (MINUS) the picked time. The unit is seconds. The damping factor GAMMA must be set individually. In case where the workflow is used you can specify different values for GAMMA per stage.
When using TW\_IND = 1 three columns are expected in PICK\_FILE for each trace with the picked time (PT), time window length plus (TWLP) and time window length minus (TWLM). In this case TWLENGTH\_PLUS and TWLENGTH\_MINUS are ignored. It is also possible to use two time windows for each trace which can be utilized with TW\_IND = 2. PICK\_FILE must then contain six columns with values for PT\_1, TWLP\_1, TWLM\_1, PT\_2, TWLP\_2 and TWLM\_2 for each trace.
......
par/workflow.txt
View file @ 04ca7c45
# VS VP RHO PRO F_FILT F_HIGH F_LOW WT J J EPRE EPSI
1 1 0 0 0.01 0 0 0 2 0 0 0 0.005
2 1 0 0 0.01 0 0 0 2 0 0 0 0.005
3 1 0 0 0.01 0 0 0 2 0 0 0 0.005
4 1 0 0 0.01 0 0 0 2 0 0 0 0.005
# VS VP RHO PRO F_FILT F_HIGH F_LOW WT J J EPRE EPSI GAMMA
1 1 0 0 0.01 0 0 0 2 0 0 0 0.005 20
2 1 0 0 0.01 0 0 0 2 0 0 0 0.005 10
3 1 0 0 0.01 0 0 0 2 0 0 0 0.005 5
4 1 0 0 0.01 0 0 0 2 0 0 0 0.005 0
src/IFOS2D.c
View file @ 04ca7c45
This diff is collapsed.
src/apply_workflow.c
View file @ 04ca7c45
......@@ -38,6 +38,7 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST
extern float JOINT_INVERSION_PSV_SH_ALPHA_RHO;
extern int EPRECOND;
extern float EPSILON_WE;
extern float GAMMA;
extern int GRAD_METHOD;
extern int WORKFLOW_STAGE;
......@@ -126,10 +127,12 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST
if(EPRECOND==0 && workflow[WORKFLOW_STAGE][12]!=0){
if(MYID==0) printf(" WARNING: EPRECOND have to be set >0 in JSON (if so, ignore this message)");
}
EPRECOND=workflow[WORKFLOW_STAGE][12];
EPSILON_WE=workflow[WORKFLOW_STAGE][13];
GAMMA=workflow[WORKFLOW_STAGE][14];
if(*LBFGS_iter_start==*iter && GRAD_METHOD==2){
if(MYID==0)printf("\n L-BFGS will be used from iteration %d on.",*LBFGS_iter_start+1);
}
}
\ No newline at end of file
}
src/fd.h
View file @ 04ca7c45
This diff is collapsed.
src/read_par_json.c
View file @ 04ca7c45
......@@ -463,7 +463,9 @@ void read_par_json(FILE *fp, char *fileinp){
declare_error("Variable FORWARD_ONLY could not be retrieved from the json input file!");
else {
if (FORWARD_ONLY==0) { /* FWI is calculated */
/* Overwrite IDX/IDY option from forward modeling (used for snapshots), interpolation for FWI not yet implemented*/
IDX=1;
IDY=1;
/* General inversion parameters */
if (get_int_from_objectlist("ITERMAX",number_readobjects,&ITERMAX,varname_list, value_list))
declare_error("Variable ITERMAX could not be retrieved from the json input file!");
......
src/readmod_elastic.c
View file @ 04ca7c45
......@@ -104,7 +104,7 @@ void readmod_elastic(float ** rho, float ** pi, float ** u){
}
if(feof(fp_vs) && feof(fp_rho)){
if(feof(fp_vs) || feof(fp_rho)){
declare_error("Model file VS or RHO is to small. Check dimensions NX*NY of file.");
}
......@@ -141,7 +141,7 @@ void readmod_elastic(float ** rho, float ** pi, float ** u){
fread(&vs, sizeof(float), 1, fp_vs);
fread(&rho, sizeof(float), 1, fp_rho);
if(!feof(fp_vs) && !feof(fp_rho)){
if(!feof(fp_vs) || !feof(fp_rho)){
declare_error("Model file VS or RHO is to big. Check dimensions NX*NY of file.");
}
fclose(fp_vs);
......
src/saveseis_glob.c
View file @ 04ca7c45
This diff is collapsed.
src/time_window.c
View file @ 04ca7c45
......@@ -139,10 +139,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times_m[1][i]+picked_times_m[2][i]){
sectiondata[i][j] = sectiondata[i][j] * taper;}
......@@ -160,10 +160,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
dummysection[i][j] = sectiondata[i][j];
......@@ -174,10 +174,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
dummysection[i][j] = sectiondata[i][j] * taper1;}
dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]);
taper2 = exp(-GAMMA*dumpc*dumpc);
taper2 = exp(-GAMMA*dumpc);
dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]);
taper3 = exp(-GAMMA*dumpd*dumpd);
taper3 = exp(-GAMMA*dumpd);
if(time>=picked_times_m[4][i]+picked_times_m[5][i]){
sectiondata[i][j] = sectiondata[i][j] * taper2;}
......@@ -196,10 +196,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT);
dumpa = (time-picked_times[i]-TWLENGTH_PLUS);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times[i]+TWLENGTH_MINUS);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times[i]+TWLENGTH_PLUS){
sectiondata[i][j] = sectiondata[i][j] * taper;}
......@@ -221,4 +221,4 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
free_vector(picked_times,1,ntr_glob);
}
} /* end of function time_window.c */
\ No newline at end of file
} /* end of function time_window.c */
src/time_window_glob.c
View file @ 04ca7c45
......@@ -101,10 +101,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times_m[1][i]+picked_times_m[2][i]){
sectiondata[i][j] = sectiondata[i][j] * taper;}
......@@ -123,10 +123,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
dummysection[i][j] = sectiondata[i][j];
......@@ -137,10 +137,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
dummysection[i][j] = sectiondata[i][j] * taper1;}
dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]);
taper2 = exp(-GAMMA*dumpc*dumpc);
taper2 = exp(-GAMMA*dumpc);
dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]);
taper3 = exp(-GAMMA*dumpd*dumpd);
taper3 = exp(-GAMMA*dumpd);
if(time>=picked_times_m[4][i]+picked_times_m[5][i]){
sectiondata[i][j] = sectiondata[i][j] * taper2;}
......@@ -159,10 +159,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT);
dumpa = (time-picked_times[i]-TWLENGTH_PLUS);
taper = exp(-GAMMA*dumpa*dumpa);
taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times[i]+TWLENGTH_MINUS);
taper1 = exp(-GAMMA*dumpb*dumpb);
taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times[i]+TWLENGTH_PLUS){
sectiondata[i][j] = sectiondata[i][j] * taper;}
......@@ -184,4 +184,4 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
free_vector(picked_times,1,ntr_glob);
}
} /* end of function time_window_glob.c */
\ No newline at end of file
} /* end of function time_window_glob.c */
src/update_s_elastic_PML_SH.c
View file @ 04ca7c45
......@@ -103,8 +103,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......@@ -159,8 +159,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......@@ -216,8 +216,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......@@ -275,8 +275,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......@@ -336,8 +336,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......@@ -399,8 +399,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT;
uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]=(f*vzy)/DT;
}
}
......