chenyk1990 / passive_imaging

This is an open-source package for passive seismic imaging. The methodology is detailed in Chen, Y., O.M. Saad, M. Bai, X. Liu, and S. Fomel, 2021, A compact program for 3D passive seismic source-location imaging, Seismological Research Letters. A PROOF version of the paper can be found at https://drive.google.com/file/d/1jAkbbAsGbpBqjzzjVogDcN7vnmmrgfqk/view?usp=sharing
4 stars 0 forks source link

how to implement it to 2d case #1

Closed datseismo closed 2 years ago

datseismo commented 2 years ago

Can this software be useful for 2d case? If so can you please suggest the correct way to do it.Thanks in advance.

chenyk1990 commented 2 years ago

It should be, but I need to prepare an example. There is an existing 2D version, which is much simpler. I need to find it out and share it with you.

datseismo @.***> 于2021年12月12日周日 11:54写道:

Can this software be useful for 2d case? If so can you please suggest the correct way to do it.Thanks in advance.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWEXEOCZ32TOADBCXFOTUQTOVNANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

chenyk1990 commented 2 years ago

Here is the 2D example.

datseismo @.***> 于2021年12月12日周日 11:54写道:

Can this software be useful for 2d case? If so can you please suggest the correct way to do it.Thanks in advance.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or unsubscribe. Triage notifications on the go with GitHub Mobile for iOS or Android.

chenyk1990 commented 2 years ago

https://github.com/chenyk1990/reproducible_research/tree/master/nonMada/microtri2d 2D example

chenyk1990 commented 2 years ago

Sure.

Setting a velocity model is easy, but I at least need to know some basics of this test dataset. For example, how it is formulated? what's the background of it? what's the target of this test? etc.. Some references of yours or related works are helpful.

datseismo @.***> 于2021年12月13日周一 11:56写道:

I got you point @chenyk1990 https://github.com/chenyk1990. I want to implement this programme to the 2d data attached https://i.fluffy.cc/6464gPqZzH8JhZb76JPdB6GnpKV0D7mb.html Can you please help me doing this, i want to set the velocity model as it is given in https://github.com/chenyk1990/reproducible_research/tree/master/nonMada/microtri2d. Thanks in advance.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-992727058, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWETS4TP23SMGZJVQRFLUQYXWJANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

chenyk1990 commented 2 years ago

That would be easy. Omar (cc'ed) will prepare a 2D imaging example for you.

datseismo @.***> 于2021年12月13日周一 12:07写道:

Target is to imaging the microseismic source directly from the 2d seismic data and here data is generated using scalo_CAE.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub, or unsubscribe. Triage notifications on the go with GitHub Mobile for iOS or Android.

chenyk1990 commented 2 years ago

When it is ready, I will let you know. Omar just went back to Egypt, needs time to set up his workstation.

On Tue, Dec 14, 2021 at 12:49 AM datseismo @.***> wrote:

Dear @chenyk1990 https://github.com/chenyk1990, can you please let me know when Omar(cc'ed) will prepare a 2D example.

— You are receiving this because you were mentioned.

Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-993211839, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWESJJAW4MRMAXYTCAILUQ3SIFANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

datseismo commented 2 years ago

Thanks...

chenyk1990 commented 2 years ago

Omar will hopefully be able to provide you with an example for his paper "Automatic waveform-based source-location imaging using deep learning extracted microseismic signals" since next Sunday.

"Additionally, can you suggest how to create a 1d velocity model and generate seismic data using it in Madagascar." I don't know why you need a 1D model. Could you please explain?

datseismo @.***> 于2021年12月24日周五 07:10写道:

Hii Dear Chen, can you please let me know when the package will be updated ( i requested before).

Additionally, can you suggest how to create a 1d velocity model and generate seismic data using it in Madagascar.

I think this part of the programme you provided creates 2d velocity model not 1d velocity model,

Flow('top',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=2100') Flow('mid',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=2500') Flow('bot',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=1800') Flow('vel','top mid bot','cat axis=1 ${SOURCES[1:2]} | smooth rect1=3 |window f1=30 n1=125|put o1=0 | put d2=%g'%d2).

Thanks in advance.

— Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-1000833718, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWEVYFFP2QCCJONE7AV3USRWLTANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you were mentioned.Message ID: @.***>

chenyk1990 commented 2 years ago

Just choose one trace in the 2d vel then you will have a 1d vel. Using sfwindowvel1.rsf

On Fri, Dec 24, 2021 at 11:01 AM datseismo @.***> wrote:

Dear Chen, My seismic data are generated using 1d velocity model and arranged the data in the 2d fashion and i want to implement rtm for the 2d data sets using my velocity model which is 1d.

Or Does the rtm run only with 2d/3d velocity model?..Thanks.

— Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-1000899172, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWEUKCMHRO4NM3CUQ4CTUSSRNTANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you were mentioned.Message ID: @.***>

chenyk1990 commented 2 years ago

Yes, correct. It is simplest way.

On Fri, Dec 24, 2021 at 11:10 AM datseismo @.***> wrote:

You mean first i have to generate 2d velocity model using Flow('top',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=2100') Flow('mid',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=2500') Flow('bot',None,'spike n1=80 n2=115 d1=7.5 d2=7.5 mag=1800') Flow('vel','top mid bot','cat axis=1 ${SOURCES[1:2]} | smooth rect1=3 |window f1=30 n1=125|put o1=0 | put d2=%g'%d2).

After that i have to do sfwindow vel1.rsf

— Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-1000902082, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWEVEFENLUNRGA6OXLSTUSSSR3ANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you were mentioned.Message ID: @.***>

chenyk1990 commented 2 years ago

Sure, no problem.

But I get more confused about your goal. Let me understand what you want. Please explain clearly what you really want to model.

  1. Explain why you want a 1D vel. I can't get why (1D vel is needed). You can suggest a reference or paper or a research proposal for me so that I know the background of your problem.

  2. What you want is not related to the 3D imaging package, it is related to basic programming of seismic modeling, or basic use of Madagascar, so, there is no need to ask questions this way. It is not efficient. Instead, you can directly contact me by email @.***

I believe what you want is easy, but I just need what you really want.

Best regards, Yangkang Chen

datseismo @.***> 于2021年12月25日周六 08:03写道:

Dear Chen sir, I tried this programme with 1d model(sfwindow) as you suggested and successfully run up to the step plotting the velocity model. Afterwords I am getting error on generating data and so....

Errors: mod2d.exe: Using pseudo-spectral... mod2d.exe: Forward modeling mod2d.exe: build/api/c/alloc.c: cannot allocate 18446744073709551144 bytes: Cannot allocate memory scons: *** [data0.rsf] Error 1 scons: building terminated because of errors.

I think somewhere I am unable to recognize the error as i am new to the Madagascar programming, so it would be very nice if you can suggest a solution for running the example code for 1d model. Thank you.

— Reply to this email directly, view it on GitHub https://github.com/chenyk1990/passive_imaging/issues/1#issuecomment-1001021819, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADHCWET2FXAP74XP7M4GAJTUSXFMPANCNFSM5J4QBPEQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you were mentioned.Message ID: @.***>

datseismo commented 2 years ago

Sorry for asking the questions here...Thanks.