leileely / fdwave3d

FDwave3D

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

FDWAVE3D IS DEVELOPED BASED ON A 2D VERSION PACKAGE OF VECTORIZED FD OPERATOR THE PACKAGE AND THE TESTS ARE CONDUCTED WITH MATLAB 2016b UNDER BOTH LINUX and WINDOWS SYSTEMS

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

CHARACTERSTICS OF THIS 3D FINITE DIFFERENCE MODELLING CODE PACKAGE ARE

1) In- time domain 
2) For- anisotropic elastic media and moment tensor sources
3) Over- staggered grid 
4) Uses- vectorized finite-difference operator
5) With- Matlab

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

CONTENT IN PACKAGE

manual.pdf:
Instruction and documentation of the code package.

FDwave: 
This directory contains all the programs & functions related to seismic modeling.

Matlab_vs_CPP:
The codes of the test for the efficiency comparison between Matlab and C++ codes.

The scripts used to reproduce the records and figures in the related paper are also contained.

SEG/EAGE overthrust model 
https://wiki.seg.org/wiki/SEG/EAGE_Salt_and_Overthrust_Models
or download via the link below (about 140 M):
https://drive.google.com/file/d/1vTemFS0poXAUMhHfea-nVGRSvuLPNa5K/view?usp=sharing

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

HOW TO RUN THE PROGRAM

More detailed instructions are presented in the manual file
The figures can be found in the related manuscript

e.g.,
Figure 2: 
Simulation3_homo_downhole.m: 		for the 3D homogeneous and isotropic model
compare_homo.m

Figure 3:
Simulation3_layer3_TI.m: 			for the layered anisotropic model
showwavefield_layer.m:			show original wavefields of the layered model (Figures 3a-c)
showrecord_layer.m:			show original seismograms of the layered model (Figures 3d-f)

Figure 4: show_overthrust.m

Figure 5:
Simulation3_overthrust_TI.m: 			for the anisotropic overthrust model 
showwavefield_overthrust.m:			show original wavefields of the overthrust model (Figures 5a-c)
showrecord_overthrust.m:			show original seismograms of the overthrust model (Figures 5d-f)

Figure 6: compare_cost.m


All the results can be directly downloaded via the link below (about 150 MB) for a quick validation.
https://drive.google.com/file/d/1mCKpKfma-oWOW9pfu3bvknVzmX1hsEFk/view?usp=sharing

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

GETTING HELP

Details about a function in the FDwave folder can be obtained by typing "help fun_name" in the command window.
More details of the package can be found in above references and the manual file.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

REFERENCE

Malkoti, A., Vedanti, N., Tiwari, R.K.: An algorithm for fast elastic wave simulation using a vectorized finite difference operator. Computers & Geosciences. 116, 23-31 (2018). https://doi.org/10.1016/j.cageo.2018.04.002.

Li, L., Tan, J., Zhang, D., Malkoti, A., Abakumov, I., Xie, Y.: FDwave3D: A MATLAB solver for the 3D anisotropic wave equation using the finite-difference method. Computational Geosciences. 2021. https://doi.org/10.1007/s10596-021-10060-3.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

CONTACT

Lei Li: [email protected] Ajay Malkoti: [email protected]

fdwave3d's Issues

DLL file missing

Hello,
It look like to dll file is missing for c-mex and I can't run it.
What can I do?

Best,
Maayan

在上边界为自由表面条件下，模型的网格点数貌似不能正确生成

李老师，你好！

如果是自由表面，nz貌似加了两次nab。

如果修改的话，是不是将“tl(1)=tl(1)+nab*dx”改为“tl(1)=tl(1)”就可以了。