WaveSimulation_Matlab
Author: zhaoqingwei
Date: April 24, 2023
email: zhaoqwei001@163.com
SIMPLE SIMPLE SIMPLE!!!
Main Purpose:
I majored in geophysics with a direction in seismic exploration. Wave Simulation is the key to professional entry, so summarize your own code and share it with you. So the codes is just simple simple simple! The kernel code is just a few lines. Industrialization code, I recommend C, CUDA, MPI.
Kernel Code:
- 2D acoustic wave (three lines)
const1=v.*v*DT*DT/DH/DH;
UU=imfilter(u2,dd);
u3=2*u2-u1+const1.*UU+s(k)*f;
UU=imfilter(u2,dd);
u3=2*u2-u1+const1.*UU+s(k)*f;
Keyword
- simple simple simple!
- 2-dimension; 3-dimension
- acoustic wave;elastic wave;surface wave
- sponge absorbing boundary condition;Split Perfectly Matched Layer;Convolutional Perfectly Matched Layer
- isotropy,anisotropy media (VTI, HTI)
- Finite Different Modeling
- Jian Ji ;Jian Ji ;Jian Ji . 2Wei ;3Wei . Sheng Bo ;Dan Xing Bo ;Mian Bo ;Xi Shou Bian Jie ;Fen Lie pml;Juan Ji pml. Ge Xiang Tong Xing ;Ge Xiang Yi Xing (VTI,HTI). You Xian Chai Fen Mo Xing .
Notes
3-dimension program cost memory num*nz*nx*ny*8/1024/1024(MB).Care about memory.
while nz=nx=ny=200,num=10; cost memory:610MB.
while nz=nx=ny=600,num=50; cost memory:80GB.
Thinks
RToax https://github.com/Rtoax
Guiting-geo https://github.com/Guiting-geo
Madagascar
seismic unix (SU)
CREWES
