vkazei / fastfwi Goto Github PK
View Code? Open in Web Editor NEWSimple frequency domain full-waveform inversion (FWI) regularized by Sobolev space norm
License: MIT License
Simple frequency domain full-waveform inversion (FWI) regularized by Sobolev space norm
License: MIT License
Dear Dr Vladimir Kazei
I use your published program fastFWI to run some simulations. And compare received data with 2D analytical results, modeling frequency is 90 kHz. The first model I used is a 1 meter by 1 meter homogeneous model, with transducers deployed along a circle and each transducers acted as both transmitters and receivers
The real part of received data at transducer locations obtained by finite difference modeling (FDM) and 2D analytical green functions are shown as follows
It can be seen that FDM data and analytical results matched very well.
However, when I increase the size of model to 2 meter (x direction) by 3 meter (y direction). And use parallel array setup, as shown in the following figure (grid spacing 2mm)
The real part of FDM wavefield and analytical wavefield (obtained by green functions) are shown as follows
It shows that with the increase of model size, the analytical results (homogeneous model) and FDM results begin to differ significantly. Such discrepancy influence the imaging results, leading to the algorithm failure at some circumstances. But I cannot figure out the possible reasons.
What is the possible reason that causes discrepancy between the FDM modeling wavefield and analytical results? Is such discrepancy caused by modeling grid size or some other reasons?
Hope to hear from your answers!
when I run the matlab code, it inform me cannot find the function F in line 134, is there need some other package?
Ps: the file is marmousi_ii_FWI.m
未定义函数或变量 'lbfgsb_wrapper'。
出错 lbfgsb (line 175)
lbfgsb_wrapper( m, x, l, u, nbd, f, g,
factr, pgtol, wa, iwa, task,iprint,...
出错 fwiFunc (line 61)
[mk, farsh, info] = lbfgsb(fh, l, u, opts );
出错 BP_tooth_FWI (line 159)
fwiResult = fwiFunc(iFWI, m.Init, D,
model, opts);
Hi,
I'm trying to run the "marmousi_II_FWI.m" example, but it shows that:
Undefined function or variable 'F'.
Error in marmousi_II_FWI (line 140)
DClean = F(m.Base,model);
Could you tell me how to solve this problem? Thanks!
Possible solution
add opts.scale_model -- flag
modify fwiFunc so that it rescales the model if needed
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.