If you mean the SEGY seismic data format I can recomend to use this project https://github.com/KirillBelovTest/CustomImportExport This package allows you to import and export seismic data into the Mathematica, to work with the trace headers and text headers. It supports both the data format floating-point and integer data. You can dowload repository and find some examples. The comments are on Russian but a think you can understand the code easily. If you would like to improve this project, can you tell us about your wishes and to participate in development.

Thank you for sharing this interesting body of work. I installed the package as per the instructions and am exploring the examples. As it happens the first one I tried, from the file "ExampleOfUse.md" has a minor issue. The code:

Get["CustomImportExport`"];
SetDirectory[$CustomImportExportDirectory];
file = FileNameJoin[{"CustomImportExport", "Resources", 
    "MarmousiModel.segy"}];
data = CustomImport[file, "SEGY"];
traces = data["Traces"];Head[traces]  (* check the form of traces *) 

SEGYElement

So the following ArrayPlot fails unless we take the second part of "traces", and now it works very nicely.

ArrayPlot[Transpose[traces[[2]]], AspectRatio -> 0.5, 
 ImageSize -> Large, PlotLegends -> Automatic, 
 FrameTicks -> Automatic, 
 PlotLabel -> "?????? Marmousi \n ?????? ????????? ?????????? ????"]

Thanks for finding the bug! I fixed this. In fact, there is a small inaccuracy in the documentation. Now you can use the function ArrayPlot like this:

data = CustomImport["path/to/file.segy", "SGY"]

ArrayPlot[data]
ArrayPlot[data["Traces"]]

More examples in the file CustomImportExport.nb

Thank you Anton for sharing these codes with me. This is exactly what I need, to be able to read SEGY data into the Mathematica. However, I will appreciate if you could help with English version of step-by-step guide to use it. I am a new user of Mathematica. Thank you once again

Thank you Anton, I have been able to do that successfully.

It was a small mistake in the documentation. Try this:

data = CustomImport["path/to/file.segy", "SGY"]

ArrayPlot[data]
ArrayPlot[data["Traces"][]]

(* data["Traces"][] - returns array of numbers *)

Or you can download latest version of the package (I updated this several days ago) and code from the ExampleOfUse will be work.

Dear Anton, I am working on velocity analysis (stacking velocity) using mathematica. However, I applied my code on the traces but it did not work. Of course I know why it did not work because my data has to be in CMP. Can the customExport expository do sorting? If not could you please advice on how to sort my data into CMP. Thank you.