search

daniel-koehn / DENISE-Black-Edition

2D time-domain isotropic (visco)elastic FD modeling and full waveform inversion (FWI) code for P/SV-waves
GNU General Public License v2.0
121 stars 66 forks source link

Error opening field data in DENISE despite reading with Madagascar-based FWI and RTM #45

Closed zhangxiaoshuotttt closed 1 year ago

zhangxiaoshuotttt commented 1 year ago

Hello, @daniel-koehn

First of all, I wanted to express my gratitude to you for recommending the NGP Course “Seismic Full Waveform Inversion” that has helped me get to this point.

I am using DENISE for FWI and RTM based on field seismic data. To prepare the data, I sorted my shot data and velocity field in Madagascar and SegyIO. Following this, I used Madagascar and Python to convert my segy data into su format. I normalized my shot data trace by trace using Madagascar and I ensured that the maximum amplitude was 1, as observed in the examples provided. After these preprocessing steps, I obtained the Source.dat and receiver.dat files and estimated the frequency of the source wavelet.

However, when I try to execute DENISE on my data, I receive an error message stating “Field data could not be opened”. I have checked the data and was able to read and plot it using Madagascar, so I guess that the issue is not with the data. Screenshot from 2023-06-09 23-02-15

When I checked the su file with Madagascar, it appears as shown in the attached picture. It appears that both of these files have incorrect headers. The top file is the forward model by DENISE, and the bottom file is my field data. Screenshot from 2023-06-09 23-11-10

Could you please provide some guidance on what might be causing this issue and how it could be resolved? Any help would be greatly appreciated. If you need any additional information regarding my data, please do not hesitate to let me know. I would be happy to provide more details. Alternatively, you can also email me at ToZhangSS@163.com.

Best regards

Thank you!

ADharaUTEXAS123007 commented 1 year ago

The software uses old su format. Use "suswapbytes" to convert the new su to old su format. DENISE will be able to read it.

zhangxiaoshuotttt commented 1 year ago

The software uses old su format. Use "suswapbytes" to convert the new su to old su format. DENISE will be able to read it.

Thank you very much! Your advice has greatly helped me in solving this problem.

zhangxiaoshuotttt commented 1 year ago

Additionally, I have another solution to this problem. When using Madagascar, we can utilize the sfsegywrite function and set the endian=n parameter.