ááá±á¬á·áá±ážáá¬ážáá°á¡á¬áž áá»á±ážáá°ážáááºáá«áááºá
áááá«ááºáž
áá®áá¬áá¶á¡ááŒá®ážááœáẠá€ááŸááºá
á¯ááᯠáá±ážáá¬ážáá²á·ááŒááºážááŒá
áºáááºá
áá»áœááºá¯ááºááááá·áºááœááºážáá±á«ááºážá ááºááá¯ááááºáá®á á¡áááºážá¡áá¬ážááá±á¬ááºážáá±á¬á¡ááœá²á·á¡á ááºážááá¯áááááŒá¯á á±ááá¯áá«áááºá áá®áá¬áá¶ááááºáááºááŸá¯á á¬áá»ááºááŸá¬ááẠKazan Cathedral ááá¯á¶ááááºááá¯á¡áá¯á¶ážááŒá¯áááºá áá¯áá¬ážááŸáááá¯ážáá»á±á¬ááºážááẠá ááá·áºáá®áá¬á áááºá á¡áááºáááá±áá¬áá»á¬ážáá²á០áá áºáá¯ááŒá áºááŒá®áž áá®áá¬áá¶áá»ááºážáááá·áº áááºážáá០á¡ááœááºááŸááºážáááºážá áœá¬ ááŒááºááœá±á·ááááºá
á
ááá·áºáá®áá¬á
áááºááŒááºááẠá
á®ážááœá¬ážáá±ážááá¹áááá¯ááºáááºáá±á«ááºááœáẠáá»á±á¬ááºážáá¬ážáá»á¬ážááŒá¬ážá០áááºáá±á¬ááºáá»á¬ážá áá«áááºáá°áá»á¬ážááᯠááœá±á·áá¯á¶áá²á·ááŒáááºá áááºážááá¯á·ááᯠáááºážáá»á±á¬ááºááẠááœáá·áºáááŒá¯áá²á·áá±á ááŸááºáá¯á¶áááºáá±á
ááºá¡ááœááºáž á¡ááŸááºáááááºáá±á¬ááºáá
á¹á
ááºážáá±ážáá»á¬áž (á¡áá¯áẠ- ááŸáááºáá¯ááºááŸáááºáá¯ááºáá¯ááºá áá±á¬áááºá Wolfram áááºá¹áá±ááá»á¬ážáá«áá±á¬ á
áá
áºáá¬áá»á¬áž) áá±ážá¡ááºáá²á·áá«áááºá áá±á·áááºá
á¬ááŸáá·áº áá±á¬áºáá®áá±á¬ááºáá»áááºáá»á¬ážááá¯áááºáž áá®áá¬áá¶á¡áá»áááºááá¬ážááœáẠááá·áºááœááºážáá¬ážáááºá á¡ááœá²á·áá²á·áá¶áá¶ááŸá¬ á¡ááá¬ááŸááá²á· áá±á¬áºáá®áá²á· áá¯áá·áºááœá±á¡ááŒá±á¬ááºáž ááŸááºáá¬ážáá¬ážááŒá®ážáá¬áž - á
á¬ážááá¯ááŸá°ážááœá±á á¡áááºážáá±á¬áºáááºá á€áááá«ááºážá¡ááá¯ááºážááŒáá·áºá á¡ááŒá
áºá¡áá»ááºááá¯ááºááá¯ááºá áááºážááá¯á¶á
á¶ááŸáá·áº áááºáá±áá¬ááẠá¡ááŒá¯ááá±á¬áá±á¬ááºáá±á¬ á
áááºáá¶á
á¬ážááŸá¯áá»á¬ážááᯠáá°áá±á¬ááºáá¬ááŒá®ááŒá
áºááŒá±á¬ááºáž á¡áá±ážáá±ážááŒá±á¬ááŒá¬ážááá¯áá«áááºá
áá»áœááºáá±á¬áºáá²á· Kirill Belov á ááŒááºáááºáá¬ážáá²á· á¡á
á®áááºáá¶á
á¬ááᯠâá¡áá¯á¶ážáá»áááá®áááºáá°ááá±áááá¯ááºáᬠááŒá¿áá¬ááœá±ááᯠááŒá±ááŸááºážááá¯á· Wolfram Mathematica ááá¯á¡áá¯á¶ážááŒá¯ááŒááºážááá¯á· áá±á«áºáá«áááºá ááá»ááºáá±áá¬á áá±á¬ááºá
ááºáááºážááœá²ááá¯ááºážááŒá¬ážá
áááºááŒá¬ááŸá¯ ááá¯á·ááá¯áẠâááŸá±ážáá±á¬ááºážááŒá
áºáá»á¬áž á
á®ážáááºážáá¬âá á¡á
á®áááºáá¶á
á¬á á¡ááŒá±á¬ááºážá¡áá¬ááẠá¡ááá¯ááºážááŸá
áºááá¯ááºáž áá«áááºáááº- ááááŠážá
áœá¬á áááŸáááá¯ááºáá±á¬ á¡ááºáááá¯áá®áááºáá»á¬áž á¡áá¯á¶ážááŒá¯ááŸá¯
ááá»ááºááŸá¬ááœá±áá±áž
ááááŠážá áœá¬ áááºááẠáá°áááá°ááá±áááá¯ááºáᬠáá±á·áá¬áá±ážááá®ážááá¯ááᯠááŒá¯áá¯ááºááẠááá¯á¡ááºáááºá áá°áááá°ááá±áááẠáá»á±á¬ááºáá¯á¶ážáá»á¬ážá áá¯ááºááá¯ááºážááá¯ááºáᬠáá¯ááºááá¹áááá»á¬ážááᯠáá±á·áá¬áá±á¬ áááá¹áá¶ááá¬ááŒá áºáááºá áá»á±á¬ááºáá±á¬ááºáá»á¬ážááœáẠáá»áŸááºá á áºá áá¶ááá¯ááºá áá»á±á¬á·áá»á±á¬ááºážáá±á¬ áá¯ááºááá¹áááá»á¬áž ááœá²ááŒá¬ážáá±á¬ááŒá±á¬áá·áºá áá°áááá°ááá±áááá¯ááºáᬠáááºá ááºáááºážáááºážáá»á¬áž ááŸááááº- áá»áŸááºá á áºááŸá¬ááœá±ááŒááºážá áá¶ááá¯ááºáá¬ááºááŸá¬ááœá±ááŒááºážá ááŒá±ááá»ááºááŸá¬ááœá±ááŒááºáž... á€áá±á¬ááºážáá«ážá ááááááºáá¯á¶á¡áá áá»áœááºá¯ááºááá¯á·ááẠááá»ááºá¡áá¹ááá¬ááºááᯠááá¯ááá¯á¡áá±ážá ááẠááœá±ážááœá±ážáá«áááºá ááá»ááºááŸá¬ááœá±áá±ážááẠáá±áá¶ááŸáá·áº ááá¬ááá¬ááºááœá±á·ááŸá¬ááœá±áá±ážá á¡ááááááºážáááºážááŒá áºáááºá á¡ááá¯áá«áááºážáááºážááẠá¡á¶áááºááœááºáá»ááŒá áºá á±áá±á¬ áá¯ááºáá«ááŸá¯áá»á¬ážá ááŸá¯á¶á·áá±á¬áºááŸá¯ááŸáá·áº áá±á·áá¬ááŸá¯á§áááá¬ááᯠáá±ážááœá²á·áá¬ážáá±á¬ áá»á±á¬ááºáá¯á¶ážáá»á¬ážá០áá¯á¶á·ááŒááºááŸá¯ááᯠááŸááºáááºážáááºááŒááºážá¡áá±á«áº á¡ááŒá±áá¶áááºá áá¯ááºáá«ááŸá¯áá»á¬ážááẠáá¯ááºážááŒá±áá±á«áºááœáẠ(ááá¯ááºážáááá¯áẠááá¯á·ááá¯áẠááá±á«ááºááœá²ááá¯ááºáá±á¬ áá¯ááºáá«ááŸá¯ááŒá áºá á±áá±á¬ áá¯ááºáá«ááŸá¯á¡áááºážá¡ááŒá áºáá»á¬áž) ááá¯á·ááá¯áẠáááºááẠ(áá±áá±áááºáá»á¬ážááŒáá·áº) ááœáẠá áááºááŸá¯ááºááŸá¬ážáá±áá«áááºá Elastic vibrations áá»á¬ážááẠááá°áá®áá±á¬ áá¯ááºááá¹ááááŸááá±á¬ á¡ááœáŸá¬áá»á¬ážá áááºááááááºáá»á¬ážááœáẠá¡áááºážááá¯ááºááŒá®áž áá±á¬ááºááŒááºáááºáᬠáá»á±á¬ááºáá¯á¡ááœááºáž ááŒáá·áºááœá¬ážáááºá áá±á¬ááºááŒááºáááºáá±á¬ááŸáá¯ááºážáá»á¬ážááẠáá»ááºááŸá¬ááŒááºááá¯á·ááŒááºááœá¬ážáᬠáá¯ááºážááŒá±áá±á«áºááŸá áá°áááá¯ááºážáá»á¬áž (áá»á¬ážáá±á¬á¡á¬ážááŒáá·áº ááœáá¯ááºáá²ááœáẠááá¯ááºážáá¶á·áá¬ážáá±á¬ áá¶ááá¯ááºáááœá±á·áá»á¬ážááŸá¯á¡áá±á«áºá¡ááŒá±áá¶á áá»áŸááºá á áºááá¯ááºážááá áºáááááá¬áá»á¬áž) ááá¯á·ááá¯áẠáááºáááºá¡ááœááºážááŸá ááá¯ááºáááá¯áá¯ááºážáá»á¬áž (piezoelectric á¡áá»áá¯ážáááºáá±á¬ááºááŸá¯á¡áá±á«áºá¡ááŒá±áá¶á) ááŸááºáááºážáááºáááºá ááŸáá¯ááºážáá»á¬ážáááºáá±á¬ááºáá¬áá±á¬á¡áá«ááœáẠáá°áááá±áá¡ááœáŸá¬áá»á¬ážá áááºáá²ááŸá¯ááᯠáá¯á¶ážááŒááºááá¯ááºáááºá
ááá»ááºáááºáᶠáááºá¹áá±á¬ááœá² áááááá¬
áá±áá±áááºááẠáá»áŸá±á¬á·áá¯ááºáá«ááŸá¯áá»á¬ážááᯠááŸá¯á¶á·áá±á¬áºáá±ážáááºá
ááŸáá¯ááºážáá»á¬ážááẠáá»á±á¬ááºáá¯á¶ážáá¯ááᯠááŒááºááœá¬ážáᬠááá¯ááºáááá¯áá¯ááºážáá»á¬ážááŒáá·áº ááŸááºáááºážáááºáá¬ážáááºá
á
ááá·áºáá®áá¬á
áááºááŸá Blagoveshchensky áá¶áá¬ážá¡áá®ážááŸá áá°áááá¯ááºááá¯ááºážááá¯ááºáá¬á
á
áºáááºáž áá¯áá±áááááºá¹áá±á¬ "Ivan Gubkin"
ááá»ááºá¡áá»ááºááŒáá±á¬áºáááº
áá»á±á¬ááºáá¯á¶ážáá»á¬ážááœáẠááá°áá®áá±á¬ áá¯ááºááá¯ááºážááá¯ááºáᬠáá¯ááºááá¹áááá»á¬ážááŸááááºá ááá»ááºáááºáá¶ááŸá¬ááœá±áá±áž á¡ááœááºá elastic áá¯ááºááá¹áááá»á¬áž ááẠá¡áááá¡á¬ážááŒáá·áº á¡áá±ážááŒá®ážááẠ- elastic vibrations ááŸáá·áº density ááá¯á·á ááŒáá·áºááœá¬ážááŸá¯ á¡ááŸáááºá á¡áááºá á¡ááœáŸá¬ááŸá áºáá¯ááẠáá°áá®áá±á¬ ááá¯á·ááá¯áẠá¡áá¬ážáá°áá¯ááºááá¹áááá»á¬ážááŸááá«á ááŸáá¯ááºážááẠáááºážááá¯á·ááŒá¬ážááŸá áááºááááááºááᯠáááááŒá¯áááááºááá¯ááºáá«á á¡ááœáŸá¬áá»á¬ážááŸá ááŸáá¯ááºážá¡ááŒááºááŸá¯ááºáž ááœá¬ááŒá¬ážáá«á á¡ááœáŸá¬áá»á¬ážá áááºááááááºááœáẠáá±á¬ááºááŒááºáááºááŸá¯ ááŒá áºáá±á«áºáááºááŒá áºáááºá áá¯ááºááá¹áááá»á¬áž ááœá¬ááŒá¬ážáá±áá±á áá±á¬ááºááŒááºáááºááŸá¯ ááŒááºážáááºáá±ááŒá áºáááºá áááºážá ááŒááºážáááºááŸá¯ááᯠáá±á¬ááºááŒááºáááºááááºáž (rc) ááŒáá·áº áá¯á¶ážááŒááºáá«áááº-
Ï ááẠáá»á±á¬ááºáá¬ážááááºáááºážáááŒá
áºááŒá®áž Îœ ááẠááŸáá¯ááºážá¡áá»ááºááŒá
áºááŒá®áž 1 ááŸáá·áº 2 ááẠá¡áá±á«áºááŸáá·áº á¡á±á¬ááºá¡ááœáŸá¬áá»á¬ážááᯠááœáŸááºááŒáááºá
á¡ááá¯ážááŸááºážáá¯á¶ážááŸáá·áº á¡áá¯á¶ážáá»á¬ážáá¯á¶áž ááá»ááºá¡áá»ááºááŒáá±á¬áºáááºáá»á¬ážáá²á០áá áºáá¯ááẠá á°ážá ááºážááŸá¯ááœá±ážáá¯ááºááŸá¯ááºážáá áºáá¯ááŒáá·áº áá±á¬ááºááŒááºáááºááŸá¯ááááºážááá¬ááºážáá áºáá¯á áá±á«ááºážá ááºááŸá¯ááááºá¡ááŒá ẠááŸááºáááºážáááºáá¬ážáá±á¬ ááá»ááºááŒá±áá¬ááᯠááá¯ááºá á¬ážááŒá¯ááá·áºá¡áá« á á¯á ááºážááŸá¯áá¯á¶á á¶ááŒá áºáááºá
áááºááŸá¬áá² s(t) - ááá»ááºááŒá±áá¬áá¶á ááá¯ááá¯áááºááŸá¬á áááºááŸááºáá¬ážáá±á¬ á¡áá¶ááœááºážáá»áááºá¡ááœááºáž áá±á¡á¬ážáá»áŸááºá
á
Ạááá¯á·ááá¯áẠáá°áááá¯ááºážááŒáá·áº ááŸááºáááºážáááºáá¬ážááá»áŸá w(t) - áá±áá±áááºá áá¯ááºáá±ážáá²á· á¡áá»ááºááŒá n(t) - áá»áááºážáá°áá¶áá¶á
á¥ááá¬áá áºáá¯á¡áá±áá²á· áá±á«ááºážá ááºáá¬ážáá²á· ááá»ááºááŒá±áá¬áá±á¬ááºááᯠááœááºááŒáá·áºáá¡á±á¬ááºá ááá»ááºááŸá¬ááœá±áá±ážááœáẠááœááºáá»ááºá áœá¬á¡áá¯á¶ážááŒá¯ááá·áº Ricker pulse ááᯠáááŠážá¡áá»ááºááŒááŸá¯á¡ááŒá Ạá¡áá¯á¶ážááŒá¯áá«áááºá
length=0.050; (*Signal lenght*)
dt=0.001;(*Sample rate of signal*)
t=Range[-length/2,(length)/2,dt];(*Signal time*)
f=35;(*Central frequency*)
wavelet=(1.0-2.0*(Pi^2)*(f^2)*(t^2))*Exp[-(Pi^2)*(f^2)*(t^2)];
ListLinePlot[wavelet, Frame->True,PlotRange->Full,Filling->Axis,PlotStyle->Black,
PlotLabel->Style["Initial wavelet",Black,20],
LabelStyle->Directive[Black,Italic],
FillingStyle->{White,Black},ImageSize->Large,InterpolationOrder->2]
áááŠáž ááá»ááºááŸá¯á¶á·áá±á¬áºááŸá¯
áá»áœááºá¯ááºááá¯á·ááẠ300 ms ááŸáá·áº 600 ms á¡áááºááœáẠáááºááááááºááŸá
áºáá¯ááᯠáááºááŸááºáááºááŒá
áºááŒá®ážá áá±á¬ááºááŒááºáááºááŸá¯ááááºážááááºážáá»á¬ážááẠáá»áááºážáá¶áá«ááºáá»á¬ážááŒá
áºáááá·áºáááº
rcExample=ConstantArray[0,1000];
rcExample[[300]]=RandomReal[{-1,0}];
rcExample[[600]]=RandomReal[{0,1}];
ListPlot[rcExample,Filling->0,Frame->True,Axes->False,PlotStyle->Black,
PlotLabel->Style["Reflection Coefficients",Black,20],
LabelStyle->Directive[Black,Italic]]
áá±á¬ááºááŒááºáááºááááºážáá»á¬áž
ááá»ááºááŒá±áá¬áá±á¬ááºááᯠááœááºáá»ááºááŒá®áž ááŒááŒáá·áºáá¡á±á¬ááºá áá±á¬ááºááŒááºáááºááŸá¯ááááºážááá¬ááºážáá»á¬ážááẠááá°áá®áá±á¬ ááááááºááá¹ááá¬áá»á¬áž ááŸááá±á¬ááŒá±á¬áá·áºá áá»áœááºá¯ááºááá¯á·ááẠááá»ááºáááºážááŒá±á¬ááºážáá±á«áºááœáẠáá±á¬ááºááŒááºáááºááŸá¯ááŸá
áºáá¯ááᯠáááŸááá«áááºá
traceExamle=ListConvolve[wavelet[[1;;;;1]],rcExample];
ListPlot[traceExamle,
PlotStyle->Black,Filling->0,Frame->True,Axes->False,
PlotLabel->Style["Seismic trace",Black,20],
LabelStyle->Directive[Black,Italic]]
ááá¯ááºáጠáá
áºáá¯ááº
á€á¥ááá¬á¡ááœááºá ááŒáá¯áááºá
á¬áááºážááœááºážááẠááá¯á¡ááºááẠ- á¡ááŸááºááááºááœáẠá¡ááœáŸá¬áá»á¬ážá á¡ááááºá¡áááºááᯠáá®áá¬ááŒáá·áº áá¯á¶ážááŒááºááŒá®áž á¡áá»áááºááá¯ááááºážá¡ááœáẠááá»ááºááŒá±áá¬áá¶ááŸá¯ááᯠááœááºáá»ááºáááºá á¡ááááºá¡áááºááᯠáá®áá¬áá»á¬ážáááºááŸááºááŒá®áž á¡ááœáŸá¬áá»á¬ážááŸá á¡ááŸáááºááá¯ááá ááá¯ááºáá±á¬ááºáá»áááºáá»á¬ážááᯠááœááºáá»ááºááŒááºážááẠááá¯ááŸááºáááá·áºáááºá á€ááá
á¹á
ááœááºá áá»áœááºá¯ááºááẠá¡áá»áááºáááºááá¯ážáá±á«áºááœáẠá¡ááœáŸá¬áá»á¬ážááᯠáá»ááºáá»ááºážáááºááŸááºáááºá
ááœááºážáááºáž áá¯áá±ááá¡ááŒá±á¬ááºáž ááŒá±á¬áááºá á¡á²áá®ááᯠáá±á·áá¬ááœá±á·ááŸáááŸá¯ááœá±ááŒá±á¬áá·áº á¡áá¬ážáá° á¡áá»áááºá á®ážáá®áž (ááá»ááºááŒá±áá¬) á¡áá»á¬ážá¡ááŒá¬ážááᯠááŸááºáááºážáááºáá«áááºá á¥ááá¬á¡á¬ážááŒáá·áºá á¡ááŸáẠ25 áá®ááá¯áá®áá¬ááŸáá·áº á¡áá»áẠ15 áá®ááá¯áá®áá¬ááŸá ááá¯ááºáá áºáá¯ááᯠáá±á·áá¬áá±á¬á¡áá«á á¡áá¯ááºááááºáá áºáá¯á á®ááẠ25x25 áá®áá¬ááŸááá±á¬ áá²ááºáá áºáá¯á á®ááᯠá¡ááŸááºá¡áá¬ážááŒá¯ááá·áºáá±áᬠ(ááá¯áá²ááºááᯠáááºáá áºáá¯áá¯áá±á«áºáááº)á áá±á¬ááºáá¯á¶ážáá±áá¬áááºážáá»ááºážááœáẠáá²ááœááºá 600000 áá«áááºáááºááŒá áºáááºá ááá°áá¬áá°áá»ááẠ1 ms ááŸáá·áº ááŸááºáááºážáááºáá»ááẠ5 á áá¹ááá·áºááŒáá·áºá áá±á¬ááºáá¯á¶ážáá±áá¬ááá¯ááºááẠ11 GB ááẠááá¯áá±áááºááŒá áºááŒá®áž áá°áááºáž âáá¯ááºááŒááºážâ áá á¹á ááºážá áá¯áááºááá¬áááẠáá¬ááŸáá·áºáá»á®áá±á¬ áá áºáá«ááá¯ááºáá»á¬áž ááŒá áºááá¯ááºáá«áááºá
áá°ááá¯á·áá²á· áááºááá¯á¡áá¯ááºáá¯ááºáá²á
á¡áá¯áẠGeologyIO
á¡áá¯ááºáááœá¶á·ááŒáá¯ážááá¯ážáááºááŸá¯á
áááºáá²á·áááºá
- ááŒá±áá¯á¶áá±áá¬ááᯠZMAP ááŸáá·áº IRAP áá±á¬áºáááºáá»á¬ážááŒáá·áº áááºááœááºážááŒááºážá
- LAS áá±á¬áºáááºáá±ááœááºážáá»á¬ážááœáẠááá¯ááºážáá¬áá»ááºáá»á¬áž áááºááœááºážááŒááºážá
- ááá»ááºááá¯ááºáá»á¬ážáá±á¬áºáááºá ááá·áºááœááºážááŒááºážááŸáá·áº á¡ááœááº
SEGY
áááºáá±á·áá»áºááᯠááá·áºááœááºážáááºá á
á¯á
ááºážáá¬ážáá±á¬ áááºáá±á·ááºá»á áá±á«ááºážáá¯ááºá
á¬áá»ááºááŸá¬ááŸá ááœáŸááºááŒá¬ážáá»ááºáá»á¬ážááᯠááá¯ááºáá¬ááá«áááºá á¡á±á¬ááºáá«áá¯ááºááᯠáááºááá·áºáá±áá¬ááœááºáááᯠexecute áá¯ááºáá«á
If[PacletInformation["GeologyIO"] === {}, PacletInstall[URLDownload[
"https://wolfr.am/FiQ5oFih",
FileNameJoin[{CreateDirectory[], "GeologyIO-0.2.2.paclet"}]
]]]
ááá¯á·áá±á¬áẠáááºáá±á·áá»áºááᯠáá¯á¶áá±ááá¯ááºááœá²ááœáẠááá·áºááœááºážáááºááŒá áºááŒá®ážá áááŸáááá¯ááºááá·áºáááºážááŒá±á¬ááºážááŸá¬ á¡á±á¬ááºáá«á¡ááá¯ááºážááŒá áºáááºá
FileNameJoin[{$UserBasePacletsDirectory, "Repository"}]
á¥ááá¬á¡áá±ááŒáá·áºá áá»áœááºá¯ááºááá¯á·ááẠpackage á á¡áááá áœááºážáá±á¬ááºáááºáá»á¬ážááᯠááŒááá«áááºá áá±á«áºááá¯ááŸá¯ááẠWolfram áá¬áá¬á áá¬ážááŸá áááºáá±á·áá»áºáá»á¬ážá¡ááœáẠá¡á ááºá¡áá¬á¡á¬ážááŒáá·áº áá¯ááºáá±á¬ááºáááº-
Get["GeologyIO`"]
áááºáá±á·ááºá»ááᯠá¡áá¯á¶ážááŒá¯á áá®ááœááºáá²á·áááºá
á¡áá°ážáááŒáá·áº ááá¯áá²á·ááá¯á·áá±á¬ááá¯ááºááẠâMarmousi.segyâ ááá¯ááºááŒá
áºááẠ- áááºážááẠFrench Petroleum Institute ááŸáá®ááœááºáá¬ážááá·áº áá°áááá±áááá¹áá áá±á«ááºážá
ááºáá¯á¶á
á¶áá
áºáá¯ááŒá
áºáááºá á€áá¯á¶á
á¶ááᯠá¡áá¯á¶ážááŒá¯á developer áá»á¬ážááẠwave field modelingá data processingá seismic trace ááŒá±á¬ááºážááŒááºááŸááºááŒááºáž á¡á
ááŸááááºááá¯á·á¡ááœáẠáááºážááá¯á·á ááá¯ááºááá¯áẠalgorithms ááᯠá
ááºážáááºáá«áááºá Marmousi áá±á¬áºáááºááá¯ááºááá¯ááºá package ááá¯ááºááá¯ááºáá±á«ááºážáá¯ááºáá¯ááºáá¬ážáá²á· repository ááŸá¬ ááááºážáááºážáá¬ážáá«áááºá ááá¯ááºááá¯ááá°áááºá¡ááœáẠá¡á±á¬ááºáá«áá¯ááºááᯠáá¯ááºáá±á¬ááºáá«á
If[Not[FileExistsQ["Marmousi.segy"]],
URLDownload["https://wolfr.am/FiQGh7rk", "Marmousi.segy"];]
marmousi = SEGYImport["Marmousi.segy"]
áááẠ- SEGYData á¡áá¬ááá¹áá¯ááᯠáááºááœááºážáá«á
SEGY áá±á¬áºáááºááœáẠáá±á·áá¬ááœá±á·ááŸááá»ááºáá»á¬ážááá¯ááºáᬠá¡áá»ááºá¡áááºá¡áá»áá¯ážáá»áá¯ážááᯠááááºážáááºážááŒááºáž áá«áááºáááºá ááááŠážá
áœá¬á á€á¡áá¬áá»á¬ážááẠá
á¬áá¬ážááŸááºáá»ááºáá»á¬ážááŒá
áºáááºá áááºážááœáẠá¡áá¯ááºáááºáá±áá¬á ááá¯ááºážáá¬ááŸá¯áá¯ááºáá±á¬ááºáá±á¬ áá¯áá¹ááá®áá»á¬ážá á¡áááºáá»á¬áž á
áááºááá¯á·ááŸáá·áºáááºáááºáá±á¬ á¡áá»ááºá¡áááºáá»á¬áž áá«áááºáááºá áá»áœááºá¯ááºááá¯á·áááá
á¹á
ááœááºá á€áá±á«ááºážá
á®ážá¡á¬áž TextHeader áá±á¬á·ááŒáá·áº áá±á¬ááºážááá¯ááŸá¯áá
áºáá¯ááŒáá·áº áá±á«áºáááºá á€áááºááŸá¬ á¡ááá¯áá»á¯á¶á·áá¬ážáá±á¬ á
á¬áá¬ážáá±á«ááºážá
á®ážááŒá
áºáááº-
Short[marmousi["TextHeader"]]
"Marmousi áá±áá¬á¡á á¯á¶ááᯠá¡ááºá áá®áá»á¯á០áá¯ááºáá¯ááºáá²á·áááº... á¡áááºážáá¯á¶áž á¡ááŒááºááŸá¯ááºáž 1500 m/s ááŸáá·áº á¡áá»á¬ážáá¯á¶áž 5500 m/s)"
"ááŒá±áá¬áá¶áá»á¬áž" áá±á¬á·ááᯠá¡áá¯á¶ážááŒá¯á ááá»ááºááŒá±áá¬áá»á¬ážááᯠáááºáá±á¬ááºááŒááºážááŒáá·áº á¡ááŸááºáááẠáá°áááá±á áá±á¬áºáááºááᯠáááºááŒáááá¯ááºááẠ(áááºáá±á·áá»áºá á¡ááºá¹áá«áááºáá»á¬ážáá²á០áá áºáá¯ááŸá¬ áá±á¬á·áá»á¬ážááẠá¡áá±ážá¡ááœáŸá¬ážáá»á¬ážáá¬ááŒá áºááŒá±á¬ááºáž)á
ArrayPlot[Transpose[marmousi["traces"]], PlotTheme -> "Detailed"]
áá±á¬áºááẠMarmousi
áááºááŸáááœááºá áááºáá±á·ááºá»ááẠááá·áºá¡á¬áž ááŒá®ážáá¬ážáá±á¬ááá¯ááºáá»á¬ážá០á¡á
áááºá¡ááá¯ááºážáá»á¬ážá¡ááœááºáž áá±áá¬áá»á¬ážááᯠáááºááá¯ááºá
á±ááŒá®áž á¡ááœááºá¡á
á¬áž áááºááááºážáá
áºáá«ááá¯ááºá¡ááááŸááá±á¬ ááá¯ááºáá»á¬ážááᯠá
á®áá¶áá±á¬ááºááœááºáá±ážááá¯ááºá
á±áááºááŒá
áºáááºá áááºáá±á·áá»áºááá¯ááºáá±á¬ááºáá»ááºáá»á¬ážááẠ.segy ááá¯á· áá±áá¬áááºááá¯á·ááŒááºážááá¯ááºáᬠáá¯ááºáá±á¬ááºáá»ááºáá»á¬ážáá«ááŸáááŒá®áž ááá¯ááºáá¡áá¯á¶ážááœáẠáá
áºá
áááºáá
áºááá¯ááºážááá·áºááœááºážáá¬ážáááºá
áá®ážááŒá¬ážá¡á¬ážááŒáá·áºá .segy ááá¯ááºáá»á¬ážá ááŸá¯ááºááœá±ážáá±á¬ááœá²á·á ááºážáá¯á¶ááŸáá·áº áá¯ááºáá±á¬ááºááá·áºá¡áá« á¡áá¯ááºááá¯ááºáá±á¬ááºááá¯ááºá áœááºážááᯠáááááŒá¯ááá·áºáááºá áá±á¬á·áá»á¬ážááŸáá·áº á¡ááœáŸááºážáá»á¬ážááᯠá¡áá¯á¶ážááŒá¯á áá áºáŠážáá»ááºážá á® ááŒá±áá¬áá¶áá»á¬ážááŸáá·áº áá±á«ááºážá á®ážáá»á¬ážááᯠáááºáá±á¬ááºááŒáá·áºááŸá¯ááá¯ááºáá¯á¶áá¬áá áááºážááá¯á·ááᯠááŒá±á¬ááºážáá²ááŒá®ážáá±á¬áẠááá¯ááºáá áºáá¯ááá¯á· áá±ážáááºáááºáž ááœáá·áºááŒá¯áá¬ážáá±á¬ááŒá±á¬áá·áº ááŒá áºáááºá GeologyIO á á¡áá±á¬ááºá¡áááºáá±á¬áºááŸá¯ááá¯ááºáᬠáááºážááá¬ááá¯ááºáᬠá¡áá±ážá áááºá¡áá»ááºá¡áááºáá»á¬ážá áœá¬ááẠá€áá±á¬ááºážáá«ážá áá±á¬ááºáá»á±á¬áºááœááºááŒá®áž áá®ážááŒá¬ážáá±á¬áºááŒáá»ááºááŸáá·áº ááá¯ááºáááºáá«áááºá
ááá»ááºááŸá¬ááœá±áá±ážááœáẠáá±á¬ááºá ááºáááºážááœá²ááŒááºážá áááºááŒá¬ááŸá¯ áááºá ááºááŸá¯
Wolfram Mathematica ááœáẠááá»ááºáá±áá¬ááᯠáááºááœááºážááá¯ááºááŸá¯ááẠá ááºážáááºáá±áá¬á¡ááœáẠbuilt-in signal processing function ááᯠá¡áá¯á¶ážááŒá¯ááá¯ááºáááºá ááá»ááºááŒá±áá¬áá áºáá¯á á®ááẠá¡áá»áááºá á®ážáá®ážáá áºáá¯ááᯠááá¯ááºá á¬ážááŒá¯áá±á¬ááŒá±á¬áá·áº áááºážááá¯á·ááá¯áá±á·áá¬ááẠá¡ááááááááá¬áá áºáá¯ááŸá¬ áá±á¬ááºá ááºáááºážááœá²ááŒááºážá áááºááŒá¬ááŒááºážááŒá áºáááºá ááá»ááºáá±áá¬á ááŒáááºááŸá¯ááºážááœá²á·á ááºážááŸá¯ááᯠááá¯ááºážááŒá¬ážá áááºááŒá¬ááẠááá¯á¡ááºáá»ááºáá»á¬ážá¡áááºá á¥ááá¬á¡á¬ážááŒáá·áº á¡á±á¬ááºáá«ááá¯á·ááᯠáá»áœááºá¯ááºááá¯á· á¡áááºáá±ážááá¯ááºáá«áááºá
- ááá°áá®áá±á¬ááŸáá¯ááºážá¡áá»áá¯ážá¡á á¬ážáá»á¬ážááᯠááá°áá®áá±á¬ááŒáááºááŸá¯ááºážááœá²á·á ááºážááŸá¯ááŒáá·áº ááœááºááŒááºááá¹ááá¬ááŸááááºá áááºážááẠááá·áºá¡á¬áž á¡áá¯á¶ážáááºáá±á¬ááŸáá¯ááºážáá»á¬ážááᯠáá®ážáá±á¬ááºážááá¯ážááŒáááºááŸáá·áº á¡ááŸá±á¬áá·áºá¡ááŸááºááŸáá¯ááºážáá»á¬ážááᯠááááŸáááºááá¯ááºá á±áá«áááºá
- porosity ááŸáá·áº saturation áá²á·ááá¯á·áá±á¬ áá»á±á¬ááºáá¯ááºááá¹áááá»á¬ážááẠááŒáááºááŸá¯ááºážááœá²á·á ááºážááŸá¯ááᯠááááá¯ááºá á±ááá¯ááºáááºá áááºážááẠá¡áá±á¬ááºážáá¯á¶ážáá¯ááºááá¹áááá»á¬ážááŸááá±á¬ áá»á±á¬ááºáá»á¬ážááᯠááœá²ááŒá¬ážáááºááŸááºááá¯ááºá á±áááºá
- ááá°áá®áá±á¬ á¡áá°ááŸááá±á¬ á¡ááœáŸá¬áá»á¬ážááẠááá°áá®áá±á¬ ááŒáááºááŸá¯ááºážáá±á¬ááºáá»á¬ážááœáẠááœá²ááœá²áá»ááºáá»á¬ážááᯠááŒá áºá á±áááºá
ááááá¡áá»ááºááẠá€áá±á¬ááºážáá«ážá ááááááºáá¯á¶ááœáẠá¡áááá¡áá»ááºááŒá áºáááºá á¡á±á¬ááºááœáẠááœá²ááŒá¬ážáá±á¬ á¡áá°ááŸááá±á¬ á¡ááœáŸá¬áá áºáá¯ááœáẠááá»ááºáááºáá¶ááá¯ááºááŒá±áá»á¬ážááᯠááœááºáá»ááºáááºá¡ááœáẠáá¯ááºá¡ááá¯ááºážá¡á áá áºáá¯ááŒá áºáááºá á€áá¯á¶á á¶ááᯠá¡ááœáŸá¬áá»á¬ážá áœá¬á០áááºáááºáá±áá±á¬ ááŸáá¯ááºážáá»á¬ážááᯠáá áºáá¯ááŸáá·áºáá áºáᯠááœáŸááºážááŒá¯á¶áá¬ážáá±á¬á¡áá«ááœáẠáááºáá±á¬ááºá áœááºáááºááŸá¯ áááºáá±á¬ááºááŸá¯áá»á¬ážááᯠááœá²ááŒááºážá áááºááŒá¬ááẠá€áá¯á¶á á¶ááᯠááá»ááºá á°ážá ááºážáá±á·áá¬ááŒááºážááœáẠá¡á ááºá¡áá¬á¡á¬ážááŒáá·áº áá±á·áá¬áá²á·áááºá
nx=200;(* Number of grid points in X direction*)
ny=200;(* Number of grid points in Y direction*)
T=2;(*Total propagation time*)
(*Velocity and density*)
modellv=Table[4000,{i,1,ny},{j,1,nx}];(* P-wave velocity in m/s*)
rho=Table[2200,{i,1,ny},{j,1,nx}];(* Density in g/cm^3, used constant density*)
Table[modellv[[150-Round[i*0.5];;,i]]=4500;,{i,1,200}];
Table[modellv[[;;70,i]]=4500;,{i,1,200}];
(*Plotting model*)
MatrixPlot[modellv,PlotLabel->Style["Model of layer",Black,20],
LabelStyle->Directive[Black,Italic]]
pinch-out ááœá²á·á
ááºážááŸá¯áá¯á¶á
á¶
áááºá¡ááœááºážááŸá ááŸáá¯ááºážá¡ááŒááºááŸá¯ááºážááẠ4500 m/s á áááºá¡ááŒáẠ4000 m/s ááŸáááŒá®áž ááááºáááºážáááŸá¬ 2200 g/cm³ ááŒá
áºááẠá ááá¯ááá¯á·áá±á¬áá¯á¶á
á¶á¡ááœááºá áá»áœááºá¯ááºááá¯á·ááẠáá±á¬ááºááŒááºáááºááŸá¯ááááºážááááºážáá»á¬ážááŸáá·áº ááá»ááºááŒá±áá¬áá»á¬ážááᯠááœááºáá»ááºáá«áááºá
rc=Table[N[(modellv[[All,i]]-PadLeft[modellv[[All,i]],201,4000][[1;;200]])/(modellv[[All,i]]+PadLeft[modellv[[All,i]],201,4500][[1;;200]])],{i,1,200}];
traces=Table[ListConvolve[wavelet[[1;;;;1]],rc[[i]]],{i,1,200}];
starttrace=10;
endtrace=200;
steptrace=10;
trasenum=Range[starttrace,endtrace,steptrace];
traserenum=Range[Length@trasenum];
tracedist=0.5;
Rotate[Show[
Reverse[Table[
ListLinePlot[traces[[trasenum[[i]]]]*50+trasenum[[i]]*tracedist,Filling->{1->{trasenum[[i]]*tracedist,{RGBColor[0.97,0.93,0.68],Black}}},PlotStyle->Directive[Gray,Thin],PlotRange->Full,InterpolationOrder->2,Axes->False,Background->RGBColor[0.97,0.93,0.68]],
{i,1,Length@trasenum}]],ListLinePlot[Transpose[{ConstantArray[45,80],Range[80]}],PlotStyle->Red],PlotRange->All,Frame->True],270Degree]
áááºáá±á¬áºáááºá¡ááœáẠááá»ááºááŒá±áá¬áá»á¬áž
á€áá¯á¶ááœááºááŒáá¬ážááá·áº ááá»ááºááŒá±áá¬áá»á¬ážá á¡á
á®á¡áá®ááᯠááá»ááºá¡ááá¯ááºážáá¯áá±á«áºáááºá áááºááœá±á·ááŒááºáááá·áºá¡ááá¯ááºážá áá±á¬ááºááŒááºáááºáá±á¬ááŸáá¯ááºážáá»á¬ážá áá»á®ááŒáá±ááŒá®ááẠá¡á
á±á¬ááá¯ááºážá áááºááŸááºáá¬ážááá·áº áá±á¬áºáááºááŸáá·áº ááŸááºážááŸááºážáááºážáááºáž ááá¯ááºáá®áá±á¬ááŒá±á¬áá·áº áááºážáá¡áááá¹áá¬ááºááœáá·áºááá¯áá»ááºááᯠá¡ááá¯ááá¯ááááá¯ááºáá±á¬á¡ááá·áºááœáẠáá¯ááºáá±á¬ááºááá¯ááºáááºá áá²ááœááºá
áá»á¬ážááᯠá¡áá±ážá
áááºááœá²ááŒááºážá
áááºááŒá¬áá«á á á០áá ááá·áºá¡áá áá²ááœááºá
áá»á¬áž ááœá²ááŒá¬ážááŒááºážáááŸááááºááᯠáááááŒá¯áááááá·áºááẠ- ááœá²á·á
ááºážááŒááºážááá±á«ááºááá¯ážááŸáá·áºá¡á±á¬ááºááŒá±á០áá±á¬ááºááŒááºáááºááŸá¯ááẠáá
áºáá¯ááŸáá·áºáá
áºáá¯ááááºáá«á 1st trace ááŸá
áááºá áá±á¬ááºááŒááºáááºááŸá¯áá»á¬áž á
áááºáááºáá±á¬ááºá
áœááºáááºáá¬áá«áááºá áá±á¬áºáááºááœááºááŸááá±á¬áºáááºážá áá±á¬ááºááŒááºáááºááŸá¯ááááºážááááºážáá»á¬ážááẠá¡áá»á¬ážááá¯ááºáááŒá±á¬ááºážáá²áá« - ááá»ááºááŒá±áá¬áá»á¬ážááẠááœá²á·á
ááºážááŸá¯áá¡áá°ááᯠááŒá±á¬ááºážáá²ááœá¬ážáááºááŸáá·áºá¡áá»áŸ áááºážááá¯á·áááŒááºážáááºááŸá¯ááᯠááŒá±á¬ááºážáá²á
á±áááºá
ááœá²á·á ááºážááŸá¯áá¡áá±á«áºááá¯ááºážáááºááááááºááŸáá±á¬ááºááŒááºáááºááŸá¯ááá¬áááá¯áá¯á¶ážáááºááŒáá·áºááŒáá«á áá¯á·á áááºážááŒá±á¬ááºáž 60 ááŸá áááºá áá±á¬ááºááŒááºáááºááŸá¯áááŒááºážáááºááŸá¯á áááºáá¬ááŒá®áž 70 áááºážááŒá±á¬ááºážááœááºáááºážáááºá¡áá»á¬ážáá¯á¶ážááŒá áºáá¬áááºá á€áááºááŸá¬ á¡ááá¯ážááŸáá·áº á¡ááœáŸá¬áá»á¬ážá á¡á±á¬ááºááŒá±á០ááŸáá¯ááºážáá»á¬ážá ááŸá±á¬áá·áºááŸááºááŸá¯ááᯠáááºááŸá¬ážá á±ááŒá®áž á¡áá»áá¯á·áá±á¬ á¡ááŒá±á¡áá±áá»á¬ážááœáẠááá»ááºááŸááºáááºážááŸá áááá¬áááºááŸá¬ážáá±á¬ ááœá²ááœá²áá»ááºáá»á¬ážááᯠááŒá áºáá±á«áºá á±áá«áááºá
ListLinePlot[GaussianFilter[Abs[traces[[All,46]]],3][[;;;;2]],
InterpolationOrder->2,Frame->True,PlotStyle->Black,
PlotLabel->Style["Amplitude of reflection",Black,20],
LabelStyle->Directive[Black,Italic],
PlotRange->All]
áááºáá¡áá±á«áºáááºá¡á
áœááºážá០áá±á¬ááºááŒááºáááºááá·áºááŸáá¯ááºážáá»á¬ážá ááá¬ááááááº
á¡áá»ááºááŒááŸá¯ááẠááŒáááºááŸá¯ááºážáááá·áºáá±á¬á¡áá«á ááŒá®ážáá¬ážáá±á¬ááœá²á·á
ááºážáá¯á¶á¡áá°áá»á¬ážááœáẠá¡ááŸá±á¬áá·áºá¡ááŸááºáá»á¬ážá
áááºáá±á«áºáá¬ááŒá®áž ááŒáááºááŸá¯ááºážááŒáá·áºááá·áºá¡áá»ááºááŒááŸá¯ááœááºá áá±ážáááºáá±á¬á¡áá°ááœáẠáááºáá±á¬ááºá
áœááºáááºááŸá¯ááŒá
áºáááºááŒá±á¬ááºáž áá¯áá¹áááááºáá«áááºá á¡á±á¬ááºáá« áá¯ááºá¡ááá¯á¡ááœá¬áá»á¬ážááẠááŒáááºááŸá¯ááºáž 35 Hzá 55 Hz ááŸáá·áº 85 Hz ááŸááá±á¬ á¡áá»ááºááŒááŸá¯ááᯠáááºáá®ážáá±ážáááºá
waveletSet=Table[(1.0-2.0*(Pi^2)*(f^2)*(t^2))*Exp[-(Pi^2)*(f^2)*(t^2)],
{f,{35,55,85}}];
ListLinePlot[waveletSet,PlotRange->Full,PlotStyle->Black,Frame->True,
PlotLabel->Style["Set of wavelets",Black,20],
LabelStyle->Directive[Black,Italic],
ImageSize->Large,InterpolationOrder->2]
ááŒáááºááŸá¯ááºáž 35 Hzá 55Hzá 85Hz ááŸááá±á¬ áááºážááŒá
áºá¡áá»ááºááŒá¡á
á¯áá
áºáá¯
ááá»ááºááŒá±áá¬áá»á¬ážááᯠááœááºáá»ááºááŒááºážááŸáá·áº áá±á¬ááºááŒááºáááºáá±áá±á¬ ááŸáá¯ááºážááœááºáá»á¬ážá ááááºáá
áºáá»á¬ážááᯠááœááºáá»ááºááŒááºážááŒáá·áºá ááá°áá®áá±á¬ ááŒáááºááŸá¯ááºážáá»á¬ážá¡ááœáẠááœá²ááœá²ááŸá¯áá
áºáá¯á¡á¬áž ááá°áá®áá±á¬ ááœá²á·á
ááºážá¡áá°áá»á¬ážááœáẠááœá±á·ááŸáááá¯ááºáááºá
tracesSet=Table[ListConvolve[waveletSet[[j]][[1;;;;1]],rc[[i]]],{j,1,3},{i,1,200}];
lowFreq=ListLinePlot[GaussianFilter[Abs[tracesSet[[1]][[All,46]]],3][[;;;;2]],InterpolationOrder->2,PlotStyle->Black,PlotRange->All];
medFreq=ListLinePlot[GaussianFilter[Abs[tracesSet[[2]][[All,46]]],3][[;;;;2]],InterpolationOrder->2,PlotStyle->Black,PlotRange->All];
highFreq=ListLinePlot[GaussianFilter[Abs[tracesSet[[3]][[All,46]]],3][[;;;;2]],InterpolationOrder->2,PlotStyle->Black,PlotRange->All];
Show[lowFreq,medFreq,highFreq,PlotRange->{{0,100},All},
PlotLabel->Style["Amplitudes of reflection",Black,20],
LabelStyle->Directive[Black,Italic],
Frame->True]
ááá°áá®áá±á¬ááŒáááºááŸá¯ááºážáá»á¬ážá¡ááœáẠáááºáá¡áá±á«áºáááºá¡á
áœááºážá០áá±á¬ááºááŒááºáááºááá·áºááŸáá¯ááºážáá»á¬ážá áá»ááºááŒá±á¬áá±á¬ááááºáá»á¬áž
áá±áá¶ááŸá¬ááœá±áá±ážááœáẠá¡áááá¡áá¯ááºáá»á¬ážáá²á០áá
áºáá¯ááŸá¬ áá±ááœááºážáááºááŒááºážá¡ááœáẠá¡áá¬ážá¡áá¬á¡ááŸááá¯á¶ážáá±áá¬áá»á¬ážááᯠá¡áá²ááŒááºááẠ(ááá¯ááá¯áááºááŸá¬ áá±ááœááºážáá»á¬ážáááºáá±á¬ááºááá·áºáá±áá¬áá»á¬ážááœáẠááŒá±ááá»ááºáá±á·áá¬ááœá±á·ááŸááá»ááºáá»á¬ážá០ááŒá
áºáá±á«áºáá¬ááŒááºážá á¡áá°á¡áá«ážááŸáá·áº áááºáááºá áá±á¬ááºáá»ááºááœá²ááá¯ááºááŸá¯ááẠá¡ááœááºá¡áá¯á¶ážáááºáá«áááºá ááá¯áá°)á ááá¯á·á¡ááŒááºá áá°áááá±áááá¹áááœáẠááŒá
áºáá±á«áºáá¬ááá·áº á¡áá°ááᯠáááááá¬áᬠááŒá±á¬ááºážáá²ááœá¬ážá
á±ááá·áº á¥áá«ááºááŒá
áºá
á±áá±á¬ á¡áá¬ááá¹áá¯áá»á¬ážáááºáž ááŸáááá¯ááºáááºá áááºážá áááºážááá¯á·ááᯠáá±á·áá¬áááºá¡ááœáẠáá±á¬ááºá
ááºáááºážááœá²ááŒááºážá
áááºááŒá¬ááŸá¯á¡á¬áž áááá±á¬ááºáá±á¬áááááá¬áá
áºáᯠááŒá
áºá
á±áááºá áá±á¬ááºážáá«ážá áá±á¬ááºá¡ááá¯ááºážááœáẠá€áá²á·ááá¯á·áá±á¬ áá°áááá±áááá¯ááºáᬠá¡áá¬ááá¹áá¯áá»á¬ážááᯠááá¯ááá¯á¡áá±ážá
áááºáá¯á¶ážáááºáá«áááºá
ááá¹ááá±á áááºážááá¯á·ááᯠáááºáááºáááá²á·ááááºážá áááºážááá¯á·ááœáẠáá¬ááá¯ááŸá¬ááááºáááºážá
áá±á¬ááºážáá«ážááœáẠááœá²ááŒááºážá áááºááŒá¬áá¬ážáá±á¬ áá á¹á ááºážáá»á¬ážááᯠá¡áá±á¬ááºááá¯ááºáá±ážáá®ážáá¬ážááœáẠáááŸááá²á·áááºá ááŒáœááºážáá»ááºáááŸá áá°ááá¯ááºážááááŒááá·áºá¡ááá¯ááºáž ááá¯áá±áááẠáá»áœááºá¯ááºááá¯á·ááá¯ááºáá¶á á¡ááá áá±áá¶áá¯ááºáá¯ááºáá¬áá±áááŒá áºáááºá ááœááºáá²á·ááá·áºáá¬á á¯ááŸá Ạ60 ááœáẠáá±áááœáẠááá¯ááºáá»á¬ážá áááºááŒáœá áœá¬ ááœá¶á·ááŒáá¯ážááá¯ážáááºááŸá¯ááᯠá áááºáá²á·áááºá áá±áá¶ááá¯ááºáá»á¬ážááŸá¬ááœá±ááŒááºážá á¡ááááááºážáááºážááŸá¬ ááá»ááºáááºáá¶ááŸá¬ááœá±áá±ážááŒá áºáááºá á€áááºááŒá±á ááŒáá¯ááºáá¯áá¬ááºáá¯á¶áá»á¬ážááᯠááŒáá·áºááááºááŸá¬ á áááºáááºá á¬ážá áá¬áá±á¬ááºážáááºá áá±ážáááºáá±á¬á¡ááá¯ááºážá¡áá¬ááŒáá·áº á ááá·áºáááºáá»á¬ážááŸáá·áº áá±áááºá¡áá»á¬ážá¡ááŒá¬ážááᯠááŸááºáá¬ážááá¯ááºáááºá ááŒá±áá¯á¶ááᯠáá»á²á·ááŒááºážááŒáá·áºá á¡á á¯ááá¯ááºá¡ááŒá¯á¶ááá¯áẠáá°ážáá±á¬áºááá·áºáá±áá¬áá»á¬ážááᯠááŒááºááœá±á·ááá¯ááºááŒá®ážá ááŒá±áá¯á¶ááᯠááá·áºáááºáá»ááºá¡áá áá»á²á·ááŒááºážááŒáá·áº ááá»ááºáááºáá¶ááá¯ááºááá·áº áááá¯ááá¯ááºáá»á¬ážá ááŸááºážáááºážáá»ááºáá»á¬ážááá¯áááºáž ááá¯ááºážááŒá¬ážááá¯ááºáá«áááºá áá±á·áá¬á á±á¬áá·áºááŒáá·áºááŸá¯áá»á¬áž ááŒá¯áá¯ááºáá²á·áááºá
Yandex ááŒá±áá¯á¶áá»á¬ážá ááŒáá¯ááºáá¯áá¯á¶ - Noyabrsk ááŒáá¯á·á§áááá¬
á¡ááœááºáá»á¬ážáá²á០áá±ááœááºáž pads ááœááºáááº
á¡áá±á¬ááºááá¯ááºáá±ážáá®ážáá¬ážá áá»ááºááŒáá·áºáá±á¬ á¡áááºááŸáá¯ááºážááœáẠ- 1 áá®ááá¯áá®áá¬á០5 áá®ááá¯áá®áá¬á០áá±áá¶áá±á¬ááºáá±á¬ áá»á±á¬ááºáá¯á¶ážáá»á¬áž ááŒá
áºáá±á«áºáááºá Jurassic ááŸáá·áº Cretaceous áá±ááºáá»á¬ážááœáẠáá±áá¶áá»á¬ážáá«áááºáá±á¬ á¡ááááá¯áááºááᯠááœá²á·á
ááºážáá²á·áááºá Jurassic áá±ááºááᯠáá¬áááºáá°áá¯ááºááŸááºááœá±ááá± áá±á¬áºáá±á¬áºáá»á¬ážáá»á¬áž ááááŒáá«áááºá
áááºáá±á¬ááº
Jurassic áá¬á
Jurassic áá±ááºááœááºá á¡áá±á¬áẠSiberia ááááºááŒá±ááẠáááºáááºáááºážááá¯ážáááºáž (ááŒá
áºáá»á¬ážááŸáá·áº áááºáááºááááºááá¯ááºážááá¯ááŒááºá áá¯ááºážááŒá±) ááŒá
áºáááºááᯠáááááŒá¯áá«á áá¬áá®á¥áá¯á áááºáá±á¬áá·áºáááºáá¬ááŸááá±á¬ááŒá±á¬áá·áºá ááá¯á¡áá»áááºá áá¯á¶ááŸááºááŸá¯áááºážááẠá€áá²á·ááá¯á·áááºááá¯ááºáááº-
Jurassic Siberia
á€áá¯á¶ááœááºá áá»áœááºá¯ááºááá¯á·á¡ááœáẠá¡áá±ážááŒá®ážáá±á¬á¡áá¬ááŸá¬ ááááá
á¹áá¬ááºáá»á¬ážááŸáá·áº ááŸááºáá»á¬ážááá¯ááºáá±á¬áºáááºáž áá±á¬ááºáá¶ááœáẠááŒá
áºááá¯á¶ááŒá
áºáááºá ááŒá
áºááẠá
á±á¬á
á±á¬á áááºááá·áºáá²á·áá±á¬ áá°áááá±áááá¯ááºáᬠá¡áá¬ááá¹áá¯áá
áºáá¯ááŒá
áºáááºá á¡ááŸááºááŸá¬ ááŒá
áºáá»á¬ážá áá¯ááºáá±á¬ááºáá»ááºááẠáá±á¬ááºážááœááºá
áœá¬ á
á®á
á®áá¬ážáá±á¬ áá²áá»á±á¬ááºáá»á¬áž á
á¯áá¯á¶áá¬á
á±ááŒá®áž áááºážáá±á¬áẠáá±áá¶á¡ááœáẠááŸá±á¬ááºáááºáá
áºáᯠááŒá
áºáá¬áááºááŒá
áºáááºá á€áá±ááŸá±á¬ááºáááºáá»á¬ážááẠáá°ážáá°ážáááºážáááºáž ááŸá¯ááºááœá±ážáá±á¬áá¯á¶ááá¹áá¬áẠ(ááŒá
áºááŒááºážááŒááºáá²á·ááá¯á·) ááœá²ááŒá¬ážááá¯ááºááŒá®áž áááºážááá¯á·ááœáẠá¡áá°ááŸáááá¯ááºááẠ- áááºážáá¬ážáá»á¬ážá¡áá®ážááœáẠá¡áá°ááẠáá±ážáááºáá±á¬áºáááºáž áá±áááºážááŒá±á¬ááºážá á¡áááºáááá¯ááŸáá·áº ááá¯áá®ážááẠááá¯á·ááá¯áẠááœá±á·áá±á¬ááºáá±á¬áá±áá¬áá»á¬ážááœáẠááá¯ážáá¬áá«áááºá ááá¯á·ááŒá±á¬áá·áº Jurassic ááœáẠááŒá
áºáá±á«áºáá¬áá±á¬ ááŒá
áºáá»á¬ážááẠááá¯á¡áá« áá¯á¶ážáá®ááá¯áá®áá¬ááá·áº á¡áááºááœáẠááŸááá±ááŒá®áž áá±áá¶ááŸá±á¬ááºáááºáá»á¬áž ááŸá¬ááœá±ááẠá¡áá¬áá»á¬áž ááŒá
áºáá¬áááºá
ááá¹ááá±á áá¯ááºáá±á¬ááºááŒááºážááŸáá·áº á áááºáá°ážáá¯á¶áá±á¬áºááŒááºážá
áá±á¬ááºážáá«ážááœááºááŒááá¬ážááá·áº ááá»ááºááá¯ááºáá¬áá á¹á ááºážáá»á¬ážááŸáá·áºáááºáááºá áá»ááºáá»ááºážááŒáá¯áááºá á¬áááºážááœááºážááŒáá«á áá¯á· - ááœá²ááŒááºážá áááºááŒá¬ááŸá¯á¡ááœááºá¡áá¯á¶ážááŒá¯ááá·áºáá±áá¬ááá¬áááẠáááá¬áááºááŸá¬ážáá±á¬ááŒá±á¬áá·áº - ááá»ááºááŒá±áá¬áá¶áá»á¬ážááá°áááºážá¡á á¯á¡áá±ážáá¡ááá¯ááºážáá áºááá¯ááºážáᬠáá±á¬ááºážáá«ážáá á¬áá¬ážááœáẠááá·áºááœááºážáá¬ážáááºá á€á¡áá¬ááẠáááºáá°áááᯠá¡áááºáá±á¬áºááŒáá« ááœááºáá»ááºááŸá¯áá»á¬ážááᯠááŒááºáááºáááºáá®ážááá¯ááºáááºááŒá áºáááºá
ááá»ááºáá±áá¬ááŸáá·áº á¡áá¯ááºáá¯ááºáá±á¬á¡áá«á áá°áááá°ááá±áááá¬ááŸááºáá áºáŠážááẠá¡áá°ážááŒá¯áá±á¬á·ááºáá²ááºááᯠá¡áá¯á¶ážááŒá¯áá±á·ááŸáááẠ(á¥ááᬠPetrel ááá¯á·ááá¯áẠParadigm) ááẠááá·áºá¡á¬áž áá±áá¬á¡áá»áá¯ážá¡á á¬ážá¡áá»áá¯ážáá»áá¯ážááœá²ááŒááºážá áááºááŒá¬ááá¯ááºááŒá®áž á¡áááºááŒá±áá±á¬ááááºáá áºá¡ááºáá¬áá±á·á áºáá áºáá¯áá«ááŸááá±á¬ ááá¯ážáááºááŸá¯áá»á¬ážá áœá¬ááŸááá±á¬ áá¯ááºáááºážáá±á«ááºážáá±á¬ááºáá»á¬ážá áœá¬ááᯠá¡áá¯á¶ážááŒá¯áá«áááºá á¡á¬ážáá¯á¶ážá¡áááºááŒá±áá±á¬áºáááºážá á€áá±á¬á·ááºáá²ááºá¡áá»áá¯ážá¡á á¬ážáá»á¬ážááŸá¬áááºáž á¡á¬ážáááºážáá»ááºáá»á¬ážááŸáááẠ- á¥ááá¬á¡á¬ážááŒáá·áºá áááºááŒáááºáá±á¬áá¬ážááŸááºážáá»á¬ážááœáẠáá±ááºáá® algorithms áá»á¬ážááᯠá¡áá±á¬ááºá¡áááºáá±á¬áºááŒááºážááẠá¡áá»áááºáá»á¬ážá áœá¬áá°áááŒá®áž á¡ááá¯á¡áá»á±á¬ááºááœááºáá»ááºááŸá¯áá»á¬ážááŒá¯áá¯ááºááẠááŒá áºááá¯ááºááŒá±áá»á¬ážááᯠáá»á¬ážáá±á¬á¡á¬ážááŒáá·áº ááá·áºáááºáá¬ážáá«áááºá ááá¯ááá¯á·áá±á¬á¡ááŒá±á¡áá±áá»áá¯ážááœááºá ááœááºáá»á°áá¬áááºá¹áá»á¬á áá áºáá»á¬ážááŸáá·áº á¡ááá·áºááŒáá·áºáááá¯ááááºážáááºážáá¬áá¬á áá¬ážáá»á¬ážááᯠá¡áá¯á¶ážááŒá¯ááẠá¡ááœááºá¡áááºááŒá±áá¬áᬠáá»ááºááŒáá·áºáá±á¬ algorithmic á¡ááŒá±áá¶ááᯠá¡áá¯á¶ážááŒá¯ááá¯ááºááŒá®áž áá áºáá»áááºáááºážááœáẠáá¯ááºááá¯ážáá¯ááºá ááºáá»á¬ážá áœá¬ááᯠáá¯ááºáá±á¬ááºááá¯ááºáááºááŒá áºáááºá á€áááºááŸá¬ Wolfram Mathematica ááŸá ááá»ááºááá¯ááºáᬠá¡áá»ááºá¡áááºáá»á¬ážááŸáá·áº áá¯ááºáá±á¬ááºááẠá¡áá¯á¶ážááŒá¯ááá·áº áááá¬áááŒá áºáááºá áá±áá¬ááŸáá·áºá¡ááŒááºá¡ááŸááºá¡ááŒááºá¡ááŸááºáá¯ááºáá±á¬ááºááŸá¯á¡ááœáẠááŒáœááºááá±á¬áá¯ááºáá±á¬ááºááá¯ááºá áœááºážááá¯áá±ážáá¬ážááŒááºážááẠáááá·áºáá»á±á¬áºáá« - áá±áá°áá»áááºáá¶áá¬ážáá±á¬áá±á¬áºáááºááŸáááºááŒááºážáá±áá»á¬á á±áááºá ááá¯áá»ááºáá±á¬ algorithms áá»á¬ážááá¯á¡áá¯á¶ážááŒá¯áᬠááŒááºááá±á¬áºáááºááá¯á·ááŒááºáááºááŒááºážá¡ááœáẠááá¯á¡áá±ážááŒá®ážáá«áááºá
á¡ááá¯ááŒá¯áá¬ážáá±á¬ á¡á
á®á¡á
á¥áºá¡ááá¯ááºážá áá»áœááºá¯ááºááá¯á·ááẠáá°áááºážááá»ááºáááºáá¶áá±áá¬ááᯠáááºááŒá®áž áááºážááá¯á·ááᯠááŒááá«áááºá
Get["GeologyIO`"]
seismic3DZipPath = "seismic3D.zip";
seismic3DSEGYPath = "seismic3D.sgy";
If[FileExistsQ[seismic3DZipPath], DeleteFile[seismic3DZipPath]];
If[FileExistsQ[seismic3DSEGYPath], DeleteFile[seismic3DSEGYPath]];
URLDownload["https://wolfr.am/FiQIuZuH", seismic3DZipPath];
ExtractArchive[seismic3DZipPath];
seismic3DSEGY = SEGYImport[seismic3DSEGYPath]
á€áááºážááŒáá·áº áá±á«ááºážáá¯ááºáá¯ááºááŒá®áž áááºááœááºážááá·áº áá±áá¬ááẠ10 ááŸáá·áº 5 áá®ááá¯áá®áᬠá¡ááá¯ááºážá¡áá¬ááŸááá±á¬ á§áááá¬áá±á«áºááœáẠááŸááºáááºážáááºáá¬ážáá±á¬ áááºážááŒá±á¬ááºážáá»á¬ážááŒá áºáááºá áá±áá¬ááᯠáá¯á¶ážáááºááŒáẠááá»ááºá á áºáááºážáááºážáááºážááᯠá¡áá¯á¶ážááŒá¯á ááá°áá«á (ááŸáá¯ááºážáá»á¬ážááᯠáá°áááá°ááááá¯ááá¯ááºáá áºáá¯á á®ááœááºááá¯ááºáá² á§áááá¬áá áºáá¯áá¯á¶ážááᯠáá áºááŒáá¯ááºáááºáááºáž ááŸááºáááºážáááºáá¬ážáááº)á ááá»ááºááá¯ááºáᬠá¡áá»ááºá¡ááẠcubes áá»á¬ážááᯠáááŸáááá¯ááºáááºá á€á¡áá¬áá»á¬ážááẠáá°áááá±ááááºáááºážáá»ááºááᯠá¡áá±ážá áááºáá±á·áá¬ááœáá·áºááŒá¯ááá·áº áá±á«ááºááá¯ááºááŸáá·áº á¡áá»á¬ážááá¯áẠá¡ááá¯ááºážáá»á¬ážááŒá áºáááºá áá¯á¶ážáááºáá¬ážáá±á¬ á¥ááá¬ááœááºá áá»áœááºá¯ááºááá¯á·ááẠáá¯á¶ážáááºááŒááºáá±áá¬ááᯠááá¯ááºááœááºáá±áá«áááºá á€áá²á·ááá¯á·áá±á¬ á á¬áá¬ážáá±á«ááºážá á®ážá០á¡áá»ááºá¡áááºá¡áá»áá¯á·ááᯠáá»áœááºá¯ááºááá¯á· áááŸáááá¯ááºáá«áááºá
StringPartition[seismic3DSEGY["textheader"], 80] // TableForm
C 1 á€ááẠGEOLOGYIO áááºáá±á·áá»áºá ááºážáááºááŸá¯á¡ááœáẠááá¯ááºááŒááá¯ááºááŒá áºáááºá
ááᯠC 2
ááᯠC 3
ááᯠC 4
C 5 DATE á¡áá¯á¶ážááŒá¯áá°á¡áááº- WOLFRAM á¡áá¯á¶ážááŒá¯áá°
C 6 á á áºáááºážá¡áááº- SIBERIA ááá áºáá±áá¬
C 7 ááá¯ááºá¡áá»áá¯ážá¡á á¬áž 3D SeismIC ááá¬á
ááᯠC 8
ááᯠC 9
C10 Z á¡áá»áá¯ážá¡á á¬áž- ááá 2200M áá±á¬ááºáá¯á¶áž 2400M
á€áá±áá¬á¡ááœá²ááẠáá±áá¬ááœá²ááŒááºážá áááºááŒá¬ááŸá¯á á¡áááá¡ááá·áºáá»á¬ážááᯠááá¯ááºááŒááẠáá»áœááºá¯ááºááá¯á·á¡ááœáẠáá¯á¶áá±á¬ááºáá«áááá·áºáááºá ááá¯ááºááŸáááŒá±áá¬áá»á¬ážááᯠáááºááá¯ááºááŸááºáááºážáááºáá¬ážááŒá®áž áááºážááá¯á·áá áºáá¯á á®ááẠá¡á±á¬ááºáá«áá¯á¶ááá¹áá¬ááºááŸáá·áºáá°ááẠ- áááºážááẠáá±á«ááºááá¯ááºáááºááá¯ážáá áºáá»áŸá±á¬áẠá¡áááºážáááºááá¯ážáá áºáá»áŸá±á¬áẠáá±á¬ááºááŒááºáááºáá±á¬ááŸáá¯ááºážáá»á¬ážá áá»ááºáá»ááºááŒáá·áºááŒáá·áºáá»á¶á·ááŸá¶á·ááŸá¯ááŒá áºáááºá
ListLinePlot[seismic3DSEGY["traces"][[100]], InterpolationOrder -> 2,
PlotStyle -> Black, PlotLabel -> Style["Seismic trace", Black, 20],
LabelStyle -> Directive[Black, Italic], PlotRange -> All,
Frame -> True, ImageSize -> 1200, AspectRatio -> 1/5]
ááá»ááºáááºáᶠá¡ááá¯ááºážááœá±áá²á ááŒá±áá¬ááœá±á
áá±á·áá¬áá¬ážáá±á¬ á§áááá¬á áŠážáááºáá¬áá
áºáá¯á
á®ááœáẠááŒá±áá¬áá»á¬ážáááºáá»áŸááŸááááºááᯠááááŸáááŒááºážááŒáá·áºá áááºááẠáá¯á¶ážáááºááŒááºáá±áᬠáááºážáá»ááºážáá
áºáá¯ááᯠáááºáá®ážááŒá®áž Image3D[] áá¯ááºáá±á¬ááºáá»ááºááᯠá¡áá¯á¶ážááŒá¯á ááŒáááá¯ááºáááº
traces=seismic3DSEGY["traces"];
startIL=1050;EndIL=2000;stepIL=2; (*кППÑЎОМаÑа Ð¥ МаÑала О кПМÑа ÑÑÑЌкО О Ñаг ÑÑаÑÑ*)
startXL=1165;EndXL=1615;stepXL=2; (*кППÑЎОМаÑа Y МаÑала О кПМÑа ÑÑÑЌкО О Ñаг ÑÑаÑÑ*)
numIL=(EndIL-startIL)/stepIL+1; (*кПлОÑеÑÑвП ÑÑаÑÑ Ð¿ÐŸ ÐŸÐžÑ Ð¥*)
numXL=(EndXL-startXL)/stepIL+1; (*кПлОÑеÑÑвП ÑÑаÑÑ Ð¿ÐŸ ÐŸÐžÑ Y*)
Image3D[ArrayReshape[Abs[traces/Max[Abs[traces[[All,1;;;;4]]]]],{numIL,numXL,101}],ViewPoint->{-1, 0, 0},Background->RGBColor[0,0,0]]
ááá»ááºáá±áᬠcube á XNUMXD áá¯á¶á (áá±á«ááºááá¯ááºáááºááá¯áž - á¡ááááºá¡áááº)
á
áááºáá«áááºá
á¬ážáá±á¬ áá°áááá±áááá¯ááºáᬠá¡ááºá¹áá«áááºáá»á¬ážááẠááŒááºážáááºáá±á¬ ááá»ááºááá¯ááºáᬠááœá²ááœá²áá»ááºáá»á¬ážááᯠáááºáá®ážáá«áá ááœáá·áºáááºážááŒááºáá¬ááŸá¯ááŸááá±á¬ áá¯á¶ááááºáá±á¬ááºáááááá¬áá»á¬ážááᯠá¡áá¯á¶ážááŒá¯ááá¯ááºáááºá ááá¯ááºáá°ážááŸá¯á "á¡áá±ážááá«áá±á¬" á§áááá¬áá»á¬ážááᯠáááŒááºááá¯ááºáá±á¬ááŒá±á¬áá·áº ááœá²ááœá²áá»ááºáá»á¬ážááá¯áᬠááŒááºááá¯ááºáááºááŒá
áºáááºá Wolfram Mathematica ááœáẠáááºážááᯠá¡áá¯á¶ážááŒá¯á áá¯ááºáá±á¬ááºááá¯ááºáááºá
data = ArrayReshape[Abs[traces/Max[Abs[traces[[All,1;;;;4]]]]],{numIL,numXL,101}];
Graphics3D[{Opacity[0.1], Raster3D[data, ColorFunction->"RainbowOpacity"]},
Boxed->False, SphericalRegion->True, ImageSize->840, Background->None]
Opacity[] ááŸáá·áº Raster3D[] áá¯ááºáá±á¬ááºáá»ááºáá»á¬ážááᯠá¡áá¯á¶ážááŒá¯á ááá»ááºáá±áᬠcube áá¯á¶
áá±á«ááºážá
ááºá¥ááá¬ááœááºáá²á·ááá¯á·á áá°áááºáž cube áá¡ááá¯ááºážáá»á¬ážááœáẠáá°áááá±ááááºááááááºá¡áá»áá¯á· (á¡ááœáŸá¬áá»á¬áž) ááᯠááœá²ááŒá¬ážáá±á¬áááºáá¬ááŸá¯ááŒáá·áº ááœá²ááŒá¬ážáááºááŸááºááá¯ááºáááºá
áá±á¬ááºá ááºáááºážááœá²ááŒááºážá áááºááŒá¬ááŸá¯á¡ááœáẠá¡ááááááááá¬ááŸá¬ Fourier á¡ááœááºááŒá±á¬ááºážááŒááºáž ááŒá áºáááºá áááºážáá¡áá°á¡áá®ááŒáá·áºá ááŒá±áá¬áá¶áá áºáá¯á á® ááá¯á·ááá¯áẠááŒá±áá¬áá¶á¡á¯ááºá á¯áá áºáá¯á á®á áá»ááºáá»ááºááŒáá·áºááŒáá·áº-ááŒáááºááŸá¯ááºážááᯠá¡áá²ááŒááºááá¯ááºáááºá ááá¯á·áá±á¬áºá áá±áá¬ááᯠááŒáááºááŸá¯ááºážááá¯ááááºážááá¯á· ááœáŸá²ááŒá±á¬ááºážááŒá®ážáá±á¬ááºá ááŒáááºááŸá¯ááºážááŒá±á¬ááºážáá²ááŸá¯ááẠáááºááá·áºá¡áá»áááºáá»á¬áž (áááºáá²áá±á¬áá±áá¬ááœááºáááºáá«) ááŸáá·áºáááºáááºááá·áº á¡áá»ááºá¡ááẠáá»á±á¬ááºáá¯á¶ážááœá¬ážáá«áááºá á¡áá»ááẠ(á¡ááááºá¡áááº) áááºááá¯ážáá±á«áºááœáẠá¡áá»ááºááŒááŒá±á¬ááºážáá²ááŸá¯áá»á¬ážááᯠáá±áá á¶áááºááŸááºááá¯ááºá á±áááºá ááŒáááºážáá±á«ááºááŸááá±á¬ Fourier á¡ááœááºááŒá±á¬ááºážááŒááºážááŸáá·áº ááŸáá¯ááºážáá¯á¶ážááŒáá¯ááœá²ááŒááºážááá¯á·ááᯠá¡áá¯á¶ážááŒá¯áá«áááºá á€áá±á¬ááºážáá«ážááẠwavelet decomposition ááá¯á¡áá¯á¶ážááŒá¯áááºá Wavelet ááœá²ááŒááºážá áááºááŒá¬ááŒááºážáááºážááá¬ááᯠ90 áá¯ááŸá áºáá»á¬ážááœáẠááá»ááºááŸá¬ááœá±áá±ážááœáẠáááºááŒáœá áœá¬á¡áá¯á¶ážááŒá¯áá¬áá²á·áááºá Windowed Fourier á¡ááœááºááŒá±á¬ááºážááŒááºážááẠá¡á¬ážáá¬áá»ááºááŸá¬ á¡áá»áááºááá¯ááºáž ááŒááºáááºááŒááºáá¬ážááŸá¯ ááá¯áá±á¬ááºážáááºáᯠáá°ááá«áááºá
á¡á±á¬ááºáá±á¬áºááŒáá« áá¯ááºá¡ááá¯ááºážá¡á ááᯠá¡áá¯á¶ážááŒá¯á ááá»ááºááŒá±áá¬áá¶áá»á¬ážáá²á០áá áºáá¯ááᯠá¡á áááºá¡ááá¯ááºážáá áºáá¯áá»ááºážá á®á¡ááŒá áºááá¯á· ááŒáá¯ááœá²ááœá¬ážá á±ááá¯ááºáááº-
cwd=ContinuousWaveletTransform[seismicSection["traces"][[100]]]
Show[
ListLinePlot[Re[cwd[[1]]],PlotRange->All],
ListLinePlot[seismicSection["traces"][[100]],
PlotStyle->Black,PlotRange->All],ImageSize->{1500,500},AspectRatio->Full,
PlotLabel->Style["Wavelet decomposition",Black,32],
LabelStyle->Directive[Black,Italic],
PlotRange->All,
Frame->True]
áá²ááœááºá
ááᯠá¡á
áááºá¡ááá¯ááºážáá»á¬ážá¡ááŒá
áºááá¯á· ááŒáá¯ááœá²á
á±áááºá
ááá°áá®áá±á¬ ááŸáá¯ááºážáá±á¬ááºááŸááá»áááºáá»á¬ážááœáẠáá±á¬ááºááŒááºáááºááŸá¯á
áœááºážá¡ááºááᯠáááºáá²á·ááá¯á· ááŒáá·áºáá±áááºááᯠá¡áá²ááŒááºáááºá scalograms (analogous to a spectrogram) ááᯠá¡áá¯á¶ážááŒá¯áááºá á
ááºážáááºážá¡áá áááºááœá±á·ááœáẠá¡á
áááºá¡ááá¯ááºážá¡á¬ážáá¯á¶ážááᯠááœá²ááŒááºážá
áááºááŒá¬ááẠáááá¯á¡ááºáá«á áá¯á¶ááŸááºá¡á¬ážááŒáá·áºá á¡áááá·áºá á¡áááºááŸáá·áº ááŒáááºááŸá¯ááºážááŒáá·áº á¡á
áááºá¡ááá¯ááºážáá»á¬ážááᯠááœá±ážáá»ááºáá¬ážáááºá
freq=(500/(#*contWD["Wavelet"]["FourierFactor"]))&/@(Thread[{Range[contWD["Octaves"]],1}]/.contWD["Scales"])//Round;
ticks=Transpose[{Range[Length[freq]],freq}];
WaveletScalogram[contWD,Frame->True,FrameTicks->{{ticks,Automatic},Automatic},FrameTicksStyle->Directive[Orange,12],
FrameLabel->{"Time","Frequency(Hz)"},LabelStyle->Directive[Black,Bold,14],
ColorFunction->"RustTones",ImageSize->Large]
áá¯á¶ááŒááºážá áá¯ááºáá±á¬ááºáá»ááºááááº
Wolfram Language ááẠwavelet transformation á¡ááœáẠfunction ááá¯á¡áá¯á¶ážááŒá¯áááºá
tracesCWD=Table[Map[Hilbert[#,0]&,Re[ContinuousWaveletTransform[traces[[i]]][[1]]][[{13,15,18}]]],{i,1,Length@traces}];
function ááá¯á¡áá¯á¶ážááŒá¯ááŒá®ážáá±á¬ááº
ááááºáá»á¬ážááᯠáá¯ááºáá±á¬áºááá¯ááºáááºáá»á¬ážááŸáá·áº áá»áŸáá±ááẠááá¯á¡ááºáá«á ááá¯á·ááá¯áẠáá±á¬ááºáááºá¡á¬áž áá±ážáá±á¬ááºááẠááá¯á¡ááºáá«áá GeologyIO áááºáá±á·áá»áºá SEGYExport[] áá¯ááºáá±á¬ááºáá»ááºááᯠááẠá¡áá¯á¶ážááŒá¯ááá¯ááºáá«áááºá
outputdata=seismic3DSEGY;
outputdata["traces",1;;-1]=tracesCWD[[All,3]];
outputdata["textheader"]="Wavelet Decomposition Result";
outputdata["binaryheader","NumberDataTraces"]=Length[tracesCWD[[All,3]]];
SEGYExport["D:result.segy",outputdata];
ဠcubes áá¯á¶ážáᯠ(ááŒáááºááŸá¯ááºážáááá·áºá á¡áááºá¡áááºááŸáá·áº ááŒáááºááŸá¯ááºážááŒáá·áº á¡á
áááºá¡ááá¯ááºážáá»á¬áž) ááŒáá·áº RGB áá±á¬á
ááºááŒááºážááᯠáá¯á¶ááŸááºá¡á¬ážááŒáá·áº áá±áá¬ááᯠá¡áá°áááœááŒááºáá±á¬ááºááẠá¡áá¯á¶ážááŒá¯áá«áááºá á¡á
áááºá¡ááá¯ááºážáá
áºáá¯á
á®ááá¯áááºážáááá¯ááºááá¯ááºá¡áá±á¬áẠ- á¡áá®á á¡á
áááºážá á¡ááŒá¬á Wolfram Mathematica ááœáẠá€áá¯ááºáá±á¬ááºáá»ááºááᯠá¡áá¯á¶ážááŒá¯á áá¯ááºáá±á¬ááºááá¯ááºáááºá
ááááºááŸá¬ áá°áááá±áááá¯ááºáᬠá¡áááá¹áá¬ááºááœáá·áºááá¯ááá¯ááºááá·áº áá¯á¶ááááºáá»á¬ážááŒá áºáááºá á¡ááá¯ááºážááœáẠááŸááºáááºážáááºáá¬ážáá±á¬ áá»ááºážáá»á¬ážááẠáá±ááŸá±á¬ááºáááºáá»á¬áž ááŒá áºááá¯ááºááŒá®áž áá±áá¶á¡áááºáá»á¬ážáá«ááŸááá±á¬ paleochannel áá»á¬ážááᯠá¡áá±ážá áááºáá±á¬áºááŒááá¯ááºá á±áá«áááºá ááá¯áá²á·ááá¯á·áá±á¬ááŒá áºá áá áºá áá±ááºáá®áá±á¬ áááºá ááºáá á¹á ááºážáá»á¬ážááᯠááŸá¬ááœá±ááŒááºážááŸáá·áº ááœá²ááŒááºážá áááºááŒá¬ááŒááºážááẠáá»áœááºá¯ááºááá¯á·á¡á¬áž áááºážááœááºáá»á¬ážá á¡áá¬ážá¡áá¬á¡áá±á¬ááºážáá¯á¶áž á¡á áááºá¡ááá¯ááºážáá»á¬ážááᯠáá¯á¶ážááŒááºááá¯ááºá á±áá«áááºá áááºážááŒá±á¬ááºážáá»á¬ážááẠáá±á¬ááºážááœááºá áœá¬ áá²áá»á±á¬áẠá¡áá°á¡áá«ážáá»á¬ážááŒáá·áº áááºááŸá¬ážááŒá®áž áá®ááá¯ááŸá±á¬ááºááẠáá±á¬ááºážááœááºáá±á¬ áá±ááŸá±á¬ááºáááºáá áºáá¯ááŒá áºáááºá "áá¬" ááœá²ááœá²áá»ááºáá»á¬ážá¡ááŒááºáááºááŸá á§áááá¬áá»á¬ážááẠáá±ááºáá®áá±ááœáŸááºážááœááºááŒááºá¡áááºáá»á¬ážááŸáá·áº áááºáá°áááºá áá±ááœáŸááºážááœááºááŒááºá¡áááºáá»á¬ážááᯠá¡áááá¡á¬ážááŒáá·áº ááœáŸá¶á·á á±ážáá»á±á¬ááºáá»á¬ážááŒáá·áº ááá¯ááºá á¬ážááŒá¯ááŒá®áž áááºážáá¯ááºáá»á¬ážááá¯á· áá°ážáá±á¬áºááŒááºážááŸá¬ áááá±á¬ááºááŸá¯ áááŸááá±á
áá±áᬠcube á RGB á¡áá»ááºá á¡áááºáááá¯ááœáẠ(áááá¯ááááºáááºááá¯á·á¡áááºážáááº) ááœá±á·áá±á¬ááºáá±áá±á¬ááŒá
áºááá¯ááŒá±áá¬áá¶ááá¯ááºáááºá
áá±áᬠcube á RGB á¡áá»ááºá áááºáááºááŒááºážááŸá¬ ááœá±á·áá±á¬ááºáá±áá²á· ááŒá
áºááᯠááŒá±áá¬áá¶ááá¯ááºáá«áááºá
á¡áá»áá¯á·ááá
á¹á
áá»á¬ážááœáẠááá»ááºáá±áá¬áá»á¬ážá á¡áááºá¡ááœá±ážááẠáááááá¬áᬠááŒááºáááºááŒááºáá¬ážáá±á¬ áá¯á¶áá»á¬ážááᯠáááŸáá
á±áá«áááºá áááºážááẠáá°áá¶áá¶áá»áŸá±á¬á·áá»áá±áž á¡ááºáááá¯áá®áááºá០á¡áá¯á¶ážááŒá¯ááá·áº á
ááºáááááá¬áá»á¬áž áááºáááºá¡áá¯ááºáááºážá
áá
áºá¡áá±á«áº áá°áááºáááºá ááá¯ááá¯á·áá±á¬á¡ááŒá±á¡áá±áá»áá¯ážááœáẠááŒá
áºááŒá±á¬ááºážá¡ááá¯ááºážá¡á
áá»á¬ážáá¬áá ááá¯ážáá»á²á·ááŒá
áºáá»á¬ážá¡á¬ážáá¯á¶ážááá¯áááºáž ááŒááºááá¯ááºáááºá
ááá»ááºáááºáá¶áá±áᬠcube (á¡áá»á¬ážááá¯ááºá¡áá»ááº) á á¡á
áááºá¡ááá¯ááºážáá¯á¶ážáá¯ááᯠRGB áá±á¬á
ááºááŒááºážá á¡ááẠá áá®ááá¯áá®áá¬ááá·áºá
Saratov á¡áá®ážááŸá Volga ááŒá
áºáááŒáá¯ááºáá¯áá¯á¶ááááº
áá±á¬ááºáá»ááº
Wolfram Mathematica ááẠááá»ááºáá±áá¬ááᯠááá¯ááºážááŒá¬ážá
áááºááŒá¬ááá¯ááºááŒá®áž áá¬ááºááá¹áá¯ááŸá¬ááœá±ááŒááºážááŸáá·áº áááºááá¯ááºááá·áº ááŒá¿áá¬áá»á¬ážááᯠááŒá±ááŸááºážááá¯ááºááŒá®áž GeologyIO áááºáá±á·ááºá»ááẠá€áá¯ááºáááºážá
ááºááᯠááá¯ááá¯á¡áááºááŒá±á
á±áááºá ááá»ááºáá±áá¬á áááºáá±á¬ááºáá¯á¶ááŸá¬ ááœááºáá»ááºááŸá¯áá»á¬ážááᯠá¡ááŸáááºááŒáŸáá·áºááẠbuilt-in áááºážáááºážáá»á¬ážááᯠá¡áá¯á¶ážááŒá¯ááŒááºáž (
source: www.habr.com