Waad ku mahadsan tahay qoraaga blogga warbixintiisa
Horudhac
Qoraalkan waxa la qoray ka dib shirku waxaana ku jira warbixin kooban oo aan bixiyay. Dhacdadan ayaa bishii June ka dhacday St. Petersburg. Anigoo tixgalinaya in aan ka shaqeeyo meel ka mid ah goobta shirka, ma caawin kari waayay ka qaybgalka xafladan. 2016 iyo 2017, waxaan dhageystay warbixinada shirka, sanadkan waxaan ka jeediyay bandhig. Marka hore, mawduuc xiiso leh (waxay iila muuqataa aniga) ayaa soo baxay, kaas oo aan ku horumarinayno , iyo marka labaad, ka dib markii daraasad dheer oo ku saabsan sharciga Federaalka Ruushka ee ku saabsan siyaasadda cunaqabataynta, ee ganacsiga aan ka shaqeeyo, ilaa laba liisanka ayaa soo muuqday. .
Inta aanan u gudbin mowduuca khudbadeyda, waxaan jecelahay in aan xuso qabanqaabada wanaagsan ee xaflada. Bogga booqashada shirku wuxuu isticmaalaa sawirka Kazan Cathedral. Cathedral waa mid ka mid ah meelaha ugu muhiimsan ee St. Petersburg, waxaana si cad looga dhex arkayaa hoolka uu shirku ka dhacay.
Albaabka laga galo Jaamacadda Dhaqaalaha ee St. Inta lagu guda jiro diiwaangelinta, xusuuso yaryar ayaa la bixiyay ( toy - bir biligleynaya, qalin, dhejisyo leh calaamadaha Wolfram). Qado iyo nasasho qaxwo ayaa sidoo kale lagu daray jadwalka shirka. Waxaan horeyba u xusay qaxwaha macaan iyo pies darbiga kooxda - cunto karisyada ayaa aad u fiican. Qaybtan hordhaca ah, waxaan jeclaan lahaa inaan ku nuuxnuuxsado in dhacdada lafteeda, qaabkeeda iyo meesha ay ku taalba ay durba keenayaan shucuur wanaagsan.
Warbixinta ay diyaariyeen aniga iyo Kirill Belov waxaa loogu magac daray "Isticmaalka Wolfram Mathematica si loo xalliyo mashaakilaadka geophysics-ka la dabaqay. Falanqaynta Spectral ee xogta seismic ama "meeshii ay ku ordi jireen wabiyadii hore." Nuxurka warbixintu waxa ay dabooshaa laba qaybood: marka hore, isticmaalka algorithms ee laga heli karo si loo falanqeeyo xogta juquraafiga, iyo marka labaad, tani waa sida loo geliyo xogta juquraafiga Wolfram Mathematica.
Sahaminta seismic
Marka hore waxaad u baahan tahay inaad safar gaaban ku sameyso geophysics-ka. Geophysics-ka waa sayniska baraa sifooyinka jireed ee dhagaxyada. Waa hagaag, tan iyo dhagaxyadu waxay leeyihiin sifooyin kala duwan: koronto, magnetic, elastic, waxaa jira habab u dhigma geophysics: rajada korantada, magnetic prospecting, seismic prospecting ... macnaha guud ee maqaalkan, waxaan kaliya ka wada hadli doonaa sahaminta seismic si faahfaahsan. Sahaminta seismic waa habka ugu muhiimsan ee lagu baadho saliidda iyo gaasta. Habka wuxuu ku salaysan yahay kicinta gariirka laastikada iyo duubista ku xigta ee jawaabta dhagxaanta ka kooban goobta daraasadda. Gariirku waa ku farxsan yahay dhulka (oo leh ilo gariir firfircoon ama aan qarxin) ama badda (qoryaha hawada leh). Gariirada laastikada ah waxay ku faafaan cufnaanta dhagaxa, iyagoo dib u celinaya oo ka tarjumaya xuduudaha lakabyada leh sifooyin kala duwan. Mowjadaha muuqda ayaa ku soo noqda oogada waxaana lagu duubaa geophones dhulka (sida caadiga ah aaladaha korantada ee ku saleysan dhaqdhaqaaqa magnetka ee ku dhex jira gariiradda) ama taleefoonnada biyaha ee badda (oo ku saleysan saameynta piezoelectric). Waqtiga imaatinka hirarka, mid ayaa xukumi kara qoto dheer ee lakabyada juqraafiyeed.
Qalabka jiidista ee seismic
Qoriga hawadu waxa uu kiciyaa gariirka laastikada ah
Mowjadaha ayaa dhex mara cufka dhagaxa waxaana lagu duubaa taleefoonnada biyaha qaada
Markabka cilmi-baadhista Juquraafiga "Ivan Gubkin" oo ku yaal tiirarka u dhow buundada Blagoveshchensky ee St. Petersburg
Qaabka calaamada seismic
Dhagaxyadu waxay leeyihiin astaamo jidheed oo kala duwan. Sahaminta seismic, guryaha laastikada ayaa ugu horayn muhiim ah - xawaaraha faafinta gariirka laastikada iyo cufnaanta. Haddii laba lakab ay leeyihiin hanti isku mid ah ama isku mid ah, markaa mawjada "ma ogaan doonto" xadka u dhexeeya. Haddii xawaaraha hirarka ee lakabyadu kala duwan yihiin, markaa milicsiga ayaa ka dhici doona xadka lakabyada. Inta badan farqiga u dhexeeya guryaha, ayaa ka sii xoog badan milicsiga. Xoogeeda waxaa lagu go'aamin doonaa isugeynta milicsiga (rc):
halka Ο ay tahay cufnaanta dhagaxa, Ξ½ waa xawaaraha hirarka, 1 iyo 2 waxay muujinayaan lakabyada sare iyo hoose.
Mid ka mid ah kuwa ugu fudud uguna badan ee la isticmaalo moodooyinka calaamadaha seismic-ka ayaa ah qaabka convolution, marka raadadka seismic-ga la duubay lagu matalo natiijada isku dhafka isku xigxiga isbahaysiga milicsiga leh garaaca garaaca wadnaha:
xagee (t) - raad seismic, i.e. wax kasta oo lagu duubay hydrophone ama geophone waqti go'an duubista, w(t) - signalka ka dhashay qoriga hawada, n(t) - qaylada random.
Aan xisaabinno raad-raac seismic-ka-samaysan tusaale ahaan. Waxaan u isticmaali doonaa garaaca garaaca Ricker, oo si weyn loogu isticmaalo sahaminta seismic, sida calaamadda bilowga ah.
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]Dareenka seismic ee hore
Waxaan dejin doonaa laba xuduud oo qoto dheer oo ah 300 ms iyo 600 ms, iyo isugeynta milicsiga waxay noqon doontaa tirooyinka aan tooska ahayn.
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]]Isku-xigta isku-xirnaanta milicsiga
Aynu xisaabino oo soo bandhigno raadadka seismic-ka. Mar haddii isugeynta milicsiga ay leeyihiin calaamado kala duwan, waxaan helnaa laba fikradood oo beddelaadda raadraaca seismic-ka.
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]]Jid la mid ah
Tusaale ahaan, waxaa lagama maarmaan ah in la sameeyo boos celin - xaqiiqda, qoto dheer ee lakabyada ayaa la go'aamiyaa, dabcan, mitirka, iyo xisaabinta raadadka seismic waxay ku dhacdaa wakhtiga domainka. Waxay noqon doontaa mid aad u sax ah in moolka la dhigo mitirka lana xisaabiyo wakhtiyada imaatinka iyadoo la ogaanayo xawaaraha lakabyada. Xaaladdan oo kale, waxaan isla markiiba dhigay lakabyada dhidibka wakhtiga.
Haddii aan ka hadalno cilmi-baarista goobta, markaa natiijada indha-indheyntan oo kale ayaa la diiwaangeliyaa tiro aad u badan oo taxane waqti isku mid ah (raadkii seismic). Tusaale ahaan, marka la baranayo goob dhererkeedu yahay 25 km iyo 15 km, halkaas oo natiijada shaqada, raad kastaa wuxuu tilmaamayaa unug cabbirkiisu yahay 25x25 mitir (unug noocaas ah waxaa loo yaqaannaa bin), diyaarinta xogta ugu dambeysa waxay ka koobnaan doontaa 600000 raad. Waqtiga muunad ee 1 ms iyo wakhtiga duubista ee 5 ilbiriqsi, faylka xogta kama dambaysta ah wuxuu ka badnaan doonaa 11 GB, mugga alaabta asalka ah ee "ceeriin" waxay noqon kartaa boqolaal gigabytes.
Sida loola shaqeeyo ?
Xirmo
Horumarinta xirmada ayaa bilaabatay derbiga VK ee kooxda taageerada ku hadasha Ruushka. Waxaa mahad leh jawaabaha bulshada, xal degdeg ah ayaa loo helay. Natiijo ahaan, waxay u kortay horumar halis ah. U dhiganta Xitaa waxaa lagu calaamadeeyay dhexdhexaadiyeyaasha. Hadda, xirmadu waxay taageertaa la shaqaynta noocyada xogta soo socota ee si firfircoon loogu isticmaalo warshadaha juqraafiga:
- soo dejinta xogta khariidada ee qaababka ZMAP iyo IRAP
- soo dejinta cabbirada ee ceelasha qaabka LAS
- galinta iyo soo saarida qaabka faylalka seismic
Si loo rakibo xirmada, waa inaad raacdaa tilmaamaha ku yaal bogga soo dejinta ee xirmada la soo ururiyey, i.e. samee koodka soo socda midkasta :
If[PacletInformation["GeologyIO"] === {}, PacletInstall[URLDownload[
"https://wolfr.am/FiQ5oFih",
FileNameJoin[{CreateDirectory[], "GeologyIO-0.2.2.paclet"}]
]]]Taas ka dib xirmada waxaa lagu dhejin doonaa galka caadiga ah, dariiqa loo heli karo sida soo socota:
FileNameJoin[{$UserBasePacletsDirectory, "Repository"}]Tusaale ahaan, waxaan muujin doonaa awoodaha ugu muhiimsan ee xirmada. Wicitaanku si caadi ah ayaa loogu sameeyaa baakooyinka ku qoran Luqadda Wolfram:
Get["GeologyIO`"]Xirmada waxaa lagu horumariyaa iyadoo la isticmaalayo . Tani waxay kuu ogolaaneysaa inaad raacdo shaqada ugu weyn ee xirmada oo leh dukumentiyo, taas oo marka la eego qaabka soo jeedinta aan ka duwanayn dukumeentiyada Wolfram Mathematica lafteeda, iyo inaad bixiso xirmada faylasha imtixaanka ee ugu horreeya ee isbarashada.
Faylka noocan oo kale ah, gaar ahaan, waa faylka "Marmousi.segy" - kani waa qaabka synthetic ee qaybta juqraafiga, kaas oo uu sameeyay Machadka Batroolka Faransiiska. Isticmaalka moodelkan, horumariyayaashu waxay tijaabiyaan algorithms-kooda qaabaynta goobta hirarka, habaynta xogta, rogida raadadka seismic, iwm. Qaabka Marmousi laftiisa waxa lagu kaydiyaa kaydka meesha xidhmada lafteeda laga soo dejiyey. Si aad u heshid faylka, socodsii koodka soo socda:
If[Not[FileExistsQ["Marmousi.segy"]],
URLDownload["https://wolfr.am/FiQGh7rk", "Marmousi.segy"];]
marmousi = SEGYImport["Marmousi.segy"]Natiijada soo dejinta - shay SEGYData
Qaabka SEGY wuxuu ku lug leeyahay kaydinta macluumaadka kala duwan ee ku saabsan indho-indheynta. Marka hore, kuwani waa faallooyin qoraal ah. Tan waxaa ku jira macluumaadka ku saabsan goobta shaqada, magacyada shirkadaha sameeyay cabbirada, iwm. Xaaladeena, madaxan waxa lagu magacaabaa codsi wata furaha TextHeader. Waa kan madaxa qoraalka oo la soo gaabiyo:
Short[marmousi["TextHeader"]]"Xogta Marmousi waxaa laga soo saaray Machadka ... xawaaraha ugu yar ee 1500 m/s iyo ugu badnaan 5500 m/s)"
Waxaad soo bandhigi kartaa qaabka dhabta ah ee juqraafiga adiga oo gelaya raadadka seismic-ka adiga oo isticmaalaya furaha "raadka" (mid ka mid ah sifooyinka xirmada ayaa ah in furayaashu ay yihiin kuwo aan kiis dareen lahayn):
ArrayPlot[Transpose[marmousi["traces"]], PlotTheme -> "Detailed"]Model Marmousi
Waqtigan xaadirka ah, xirmada waxay sidoo kale kuu ogolaaneysaa inaad ku shubto xogta qaybo ka mid ah faylal waaweyn, taasoo suurtogal ka dhigaysa in la farsameeyo faylasha kuwaas oo cabbirkoodu gaari karo tobanaan gigabytes. Xirmada shaqadeeda waxa kale oo ka mid ah hawlaha u dhoofinta xogta .segy iyo qayb ku lifaaqida dhamaadka faylka.
Si gooni ah, waxaa habboon in la xuso shaqeynta xirmada marka la shaqeynayo qaabdhismeedka adag ee faylasha .segy. Maadaama ay kuu ogolaanayso inaadan helin oo kaliya raadadka shakhsi ahaaneed iyo madaxyada adigoo isticmaalaya furayaasha iyo tusmooyinka, laakiin sidoo kale inaad bedesho iyaga ka dibna u qor faylka. Qaar badan oo ka mid ah tafaasiisha farsamada ee hirgelinta GeologyIO ayaa ka baxsan baaxadda maqaalkan waxayna u badan tahay inay mudan yihiin sharraxaad gaar ah.
Muhiimadda falanqaynta spectral ee sahaminta seismic
Awoodda lagu soo dejiyo xogta seismic-ka Wolfram Mathematica waxay kuu ogolaanaysaa inaad adeegsato shaqaynta habaynta calaamadaha ku dhex dhisan ee xogta tijaabada ah. Maadaama raad kasta oo seismic ahi uu ka dhigan yahay taxane wakhti ah, mid ka mid ah qalabka ugu muhiimsan ee lagu baranayo iyaga waa falanqaynta spectral. Waxaa ka mid ah shuruudaha lagu falanqeynayo inta jeer ee xogta seismic-ka, waxaan magacaabi karnaa, tusaale ahaan, kuwan:
- Noocyada kala duwan ee hirarka waxaa lagu gartaa halabuurka soo noqnoqda ee kala duwan. Tani waxay kuu ogolaaneysaa inaad muujiso hirarka faa'iidada leh oo aad xakameyso hirarka faragelinta.
- Guryaha dhagaxa sida porosity iyo saturation waxay saameyn karaan inta jeer ee ka kooban. Tani waxay suurtogal ka dhigaysaa in la aqoonsado dhagaxyada leh hantida ugu fiican.
- Lakabyo leh dhumucyo kala duwan ayaa keena cillado kala duwan oo soo noqnoqda.
Qodobka saddexaad waa kan ugu weyn ee dulucda maqaalkan. Hoos waxaa ku yaal jajab kood ah oo loogu talagalay xisaabinta raadadka seismic ee kiiska lakabka leh dhumucdiisuna kala duwan tahay - moodeel jeex ah. Qaabkan waxaa dhaqan ahaan lagu daraaseeyaa sahaminta seismic si loo falanqeeyo saamaynta faragelinta marka mawjadaha ka muuqda lakabyo badan ay is dul saaran yihiin.
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]]Qaabka samaynta qanjaruufo
Xawaaraha mawjadda gudaha jeexdintu waa 4500 m/s, ka baxsan galka 4000 m/s, cufnaantana waxa loo malaynayaa inay joogto 2200 g/cmΒ³. Qaabka noocan oo kale ah, waxaanu xisaabinnaa isugeynta milicsiga iyo raadadka seismic.
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]Raadinta seismic ee moodka wedge
Isku xigxiga raadadka seismic ee shaxanka lagu muujiyay waxa loo yaqaan qaybta seismic. Sida aad arki karto, tafsiirkeeda ayaa sidoo kale lagu fulin karaa heer dareen leh, maadaama joomatari ee hirarka muuqda ay si cad u waafaqsan tahay qaabka hore loo cayimay. Haddii aad si faahfaahsan u falanqeyso raadadka, waxaad ogaan doontaa in raadadka laga bilaabo 1 ilaa 30 aysan kala duwanayn - milicsiga saqafka sare ee samaynta iyo ka hooseeya midba midka kale ma isku dhejiyo. Laga bilaabo raadraaca 31aad, milicsigu wuxuu bilaabayaa inuu farageliyo. Iyo, in kasta oo qaabka, isku-xirnaanta milicsiga aysan isbeddelin si toos ah - raadadka seismic waxay beddelaan xooggooda marka dhumucda lakabka isbeddelo.
Aynu tixgelinno baaxadda milicsiga ee ka imanaya soohdinta sare ee samaynta. Laga bilaabo dariiqa 60-aad, xoogga milicsiga wuxuu bilaabmaa inuu kordho oo dariiqa 70aad wuxuu noqdaa ugu badnaan. Tani waa sida faragelinta mowjadaha ka imanaya saqafka iyo hoose ee lakabyada ay isu muujiyaan, taas oo u horseedaysa xaaladaha qaarkood cillado muhiim ah oo ku saabsan rikoodhka seismic-ka.
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]Sawirka baaxadda hirarka muuqata ee ka imanaya cidhifka sare ee jeexjeexa
Waa wax macquul ah in marka calaamaduhu uu hooseeyo-soo noqnoqoshada, faragelintu waxay bilaabataa inay ka soo baxdo dhumucyada dhismaha ee waaweyn, iyo xaaladda calaamadda soo noqnoqda sare, faragelintu waxay ku dhacdaa dhumuc yar. Koodhka soo socdaa wuxuu abuuraa calaamad leh soo noqnoqoshada 35 Hz, 55 Hz iyo 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]Calaamadaha isha oo leh soo noqnoqoshada 35 Hz, 55Hz, 85Hz
Anagoo xisaabinayna raadadka seismic-ka iyo jaantusyada garaafyada baaxadda hirarka milicsanaya, waxaynu arki karnaa in soo noqnoqda kala duwanaansho aan caadi ahayn lagu arkay dhumucyo kala duwan.
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]Sawirrada baaxadda hirarka muuqata ee ka imanaya cidhifka sare ee jeexjeexyada kala duwan
Awoodda lagu soo gebogebeynayo dhumucda dhismaha ee natiijooyinka indha-indheynta seismic waa mid aad u faa'iido leh, sababtoo ah mid ka mid ah hawlaha ugu muhiimsan ee sahaminta saliidda ayaa ah in la qiimeeyo qodobbada ugu rajo-gelinta badan ee dhigista ceelka (ie, meelaha la sameeyay ka dhumuc weyn). Intaa waxaa dheer, qaybta juqraafiga waxaa jiri kara walxo asalkoodu uu keeno isbeddel xooggan oo dhumucdiisuna tahay sameynta. Tani waxay ka dhigaysaa falanqaynta spectral qalab wax ku ool ah oo lagu barto iyaga. Qaybta soo socota ee maqaalka waxaan si faahfaahsan uga fiirsan doonaa walxaha juqraafiga ah.
Xogta tijaabada ah. Halkee ka keentay iyo waxa laga raadiyo iyaga?
Qalabka lagu falanqeeyay maqaalka ayaa laga helay Galbeedka Siberian. Gobolka, sida uu qof kastaa og yahay oo aan laga reebin malaha, waa gobolka ugu weyn ee laga soo saaro shidaalka dalkeena. Horumarinta firfircoon ee kaydka ayaa ka bilaabmay gobolka 60-meeyadii qarnigii la soo dhaafay. Habka ugu muhiimsan ee lagu raadinayo kaydka saliida waa sahaminta seismic. Waxaa xiiso leh in la eego sawirada dayax-gacmeedka ee dhulkan. Miisaan yar, waxaad ogaan kartaa tiro aad u badan oo dhiiqo iyo harooyin ah, adoo ballaarinaya khariidada, waxaad arki kartaa goobaha ceelasha laga qodayo kooxda, iyo adigoo khariirada ku ballaarinaya xadka, waxaad sidoo kale kala saari kartaa meelaha bannaan ee profiles-ka ay weheliyaan seismic indho-indhayn ayaa la sameeyay.
Sawirka dayax-gacmeedka ee khariidadaha Yandex - agagaarka magaalada Noyabrsk
Shabakad suuf ah oo ceelasha ah oo ku yaal mid ka mid ah beeraha
Dhagaxyada saliidda leh ee Galbeedka Siberian waxay ku dhacaan qoto dheer oo kala duwan - min 1 km ilaa 5 km. Mugga ugu weyn ee dhagaxyada ay ku jiraan saliidda ayaa la sameeyay xilliyadii Jurassic iyo Cretaceous. Xilliga Jurassic waxaa laga yaabaa in dad badani u yaqaanaan filimka magaca isku midka ah. aad ayuu uga duwanaa kan casriga ah. Encyclopedia Britannica waxa uu leeyahay taxane paleomaps ah oo lagu garto xilli kasta.
Hada
xilliga Jurassic
Fadlan ogow in wakhtiyadii Jurassic, dhulka Galbeedka Siberian uu ahaa xeeb badeed (dhul ay ka gudbaan webiyada iyo badda gacmeed). Maadaama cimiladu ahayd mid raaxo leh, waxaan u qaadan karnaa in muuqaalka caadiga ah ee wakhtigaas uu u ekaa sidan:
Jurassic Siberia
Sawirkan, waxa muhiimka noo ah maaha wax aad u badan xoolaha iyo shimbiraha, laakiin sawirka webiga ee gadaasha. Wabigu waa isla shayga juqraafiga ah ee aan mar hore istaagnay. Xaqiiqdu waxay tahay in dhaqdhaqaaqa wabiyadu uu ogolaado in si fiican loo kala sooco dhagaxyada ciiddu, taas oo markaa noqon doonta kaydka saliidda. Kaydadkani waxay yeelan karaan qaab cakiran oo adag (sida sariirta webiga) waxayna leeyihiin dhumuc kala duwan - bangiyada agtooda dhumucdiisu waa yar tahay, laakiin waxay u dhowdahay bartamaha kanaalka ama meelaha dhexdhexaadka ah way kordhisaa. Haddaba, wabiyada laga sameeyay Jurassic waxay hadda ku jiraan qoto dheer ilaa saddex kilomitir waana shayga raadinta kaydka saliidda.
Xogta tijaabada ah. Habaynta iyo aragga
Aynu isla markiiba samayno boos celin ku saabsan walxaha seismic-ka ah ee lagu muujiyay maqaalka - sababtoo ah xaqiiqda ah in qadarka xogta loo isticmaalo falanqaynta ay muhiim tahay - kaliya qayb ka mid ah jaangooyooyinka asalka ah ee raadraaca seismic ayaa lagu daray qoraalka maqaalka. Tani waxay u oggolaan doontaa qof kasta inuu soo saaro xisaabaadka kor ku xusan.
Markaad la shaqeyneyso xogta seismic, geophysicist wuxuu badanaa isticmaalaa software khaas ah (waxaa jira dhowr hoggaamiyeyaal warshadeed kuwaas oo horumarkooda si firfircoon loo isticmaalo, tusaale ahaan Petrel ama Paradigm), taas oo kuu ogolaaneysa inaad falanqeyso noocyada kala duwan ee xogta waxayna leedahay interface garaafyo ku habboon. In kasta oo ay ku habboon yihiin, noocyada softiweerku sidoo kale waxay leeyihiin cilladahooda - tusaale ahaan, hirgelinta algorithms casriga ah ee noocyada xasilloon waxay qaadataa waqti badan, iyo suurtogalnimada xisaabinta otomatiga ah ayaa badanaa xaddidan. Xaaladdan oo kale, waxay noqotaa mid aad u habboon in la isticmaalo hababka xisaabta kumbuyuutarka iyo luqadaha barnaamijyada heerka sare ah, taas oo u oggolaanaysa isticmaalka saldhigga ballaaran ee algorithm iyo, isla mar ahaantaana, qaadashada wax badan oo joogto ah. Tani waa mabda'a loo isticmaalo in lagu shaqeeyo xogta seismic ee Wolfram Mathematica. Waa wax aan habooneyn in la qoro shaqeyn qani ah oo loogu talagalay shaqada isdhexgalka ee xogta - waxaa aad muhiim u ah in la hubiyo soo dejinta qaabka guud ee la aqbali karo, iyada oo la adeegsanayo algorithms-ka la rabo iyaga oo dib ugu soo celinaya qaab dibadda ah.
Iyadoo la raacayo qorshaha la soo jeediyay, waxaanu ku shubi doonaa xogta asalka ah ee seismic waxaanan soo bandhigi doonaa :
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]Xogta la soo dejiyey ee lagu soo dejiyey habkan ayaa ah waddooyinka laga diiwaan geliyey dhul cabbirkiisu yahay 10 ilaa 5 kiiloomitir. Haddii xogta lagu helo habka sahanka seismic-ga saddex-geesoodka ah (mawjadaha lama diiwaan gelinin muuqaalada juquraafi ahaaneed ee shakhsi ahaaneed, laakiin dhammaan aagga isku mar), waxa suurtogal ah in la helo xogta seismic. Kuwani waa walxo saddex-cabbir ah, qaybo toosan iyo kuwo toosan oo u oggolaanaya daraasad faahfaahsan oo ku saabsan deegaanka juqraafi. Tusaalaha la tixgaliyay, waxaanu la macaamilaynaa xogta saddex-geesoodka ah. Waxaan ka heli karnaa xoogaa macluumaad ah madaxa qoraalka, sida tan
StringPartition[seismic3DSEGY["textheader"], 80] // TableFormC 1 Tani waa faylalka DEMO ee Tijaabada Xidhmada GEOLOGYIO
C 2
C 3
C 4
C 5 TAARIIKH MAGACA ISMAALIYEED: WOLFRAM USER
C 6 MAGACYADA SUUGAANTA: MEEL SIBERIA ah
C 7 FILE NOOCA 3D mugga seismic
C 8
C 9
C10 Z RANGE: 2200M KOWAAD 2400M
Qalabkan xogta ayaa nagu filnaan doona si aan u muujino marxaladaha ugu muhiimsan ee falanqaynta xogta. Raad-raacyada faylka ayaa si isdaba-joog ah loo diiwaangeliyaa mid kasta oo iyaga ka mid ahina waxay u eg yihiin shaxanka soo socda - tani waa qaybinta baaxadda hirarka milicsiga ee dhinaca dhidibka toosan ( dhidibka qoto dheer).
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]Mid ka mid ah raadraaca qaybta seismic-ka
Ogaanshaha inta raadad ee ku yaal jiho kasta oo aagga la bartay, waxaad abuuri kartaa xog saddex-geesood ah oo aad soo bandhigi kartaa adigoo isticmaalaya 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]]Sawirka XNUMXD ee cube xogta seismic. ( dhidibka taagan - qoto dheer)
Haddii sifooyinka juqraafiga ee xiisaha ay abuuraan cillado seismic ah oo aad u daran, markaa qalab muuqaal ah oo leh daahfurnaan ayaa la isticmaali karaa. Meelaha "aan muhiimka ahayn" ee duubista ayaa laga dhigi karaa kuwo aan la arki karin, oo ka tagaya kaliya cillado muuqda. Wolfram Mathematica tan waxaa lagu samayn karaa iyadoo la isticmaalayo ΠΈ .
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]Sawirka seismic data cube iyadoo la isticmaalayo Opacity[] iyo Raster3D[] functional
Sida tusaalaha synthetic, qaybaha cube-ga asalka ah mid ayaa lagu aqoonsan karaa qaar ka mid ah xudduudaha juqraafiga (lakabyada) oo leh gargaar doorsooma.
Qalabka ugu muhiimsan ee falanqaynta spectral waa isbedelka Fourier. Caawinteeda, waxaad qiimeyn kartaa xajmiga-soo noqnoqda ee raad kasta ama koox kasta oo raad raac ah. Si kastaba ha ahaatee, ka dib markii xogta loo wareejiyo domainka soo noqnoqda, macluumaadka ayaa lumaya wakhtigee (akhri qoto dheer) inta jeer ee isbeddelka. Si loo awoodo in loo habeeyo isbeddelada ishaarada ee dhidibka wakhtiga (qoto dheer), is beddelka Fourier ee daaqada leh iyo furfurnaanta mawjada ayaa la isticmaalaa. Maqaalkani waxa uu isticmaalaa mawjada burburka. Tiknoolajiyada falanqaynta mawjada ayaa bilaabatay in si firfircoon loogu isticmaalo sahaminta seismic 90-meeyadii. Faa'iidada ka badan isbeddelka Fourier ee daaqada la furay ayaa loo arkaa inay tahay xallinta wakhtiga wanaagsan.
Adigoo isticmaalaya jajabka koodka soo socda, waxaad u kala jajabin kartaa mid ka mid ah raadadka seismic-ka qaybo gaar ah:
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]Burburinta raadka qaybo ka mid ah
Si loo qiimeeyo sida tamarta milicsiga loo qaybiyo wakhtiyada imaatinka hirarka kala duwan, scalograms (oo la mid ah spectrogram) ayaa la isticmaalaa. Sida caadiga ah, ficil ahaan looma baahna in la falanqeeyo dhammaan qaybaha. Caadi ahaan, qaybaha hoose, dhexe iyo kuwa soo noqnoqda ayaa la doortaa.
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]Scalogram Natiijada shaqada
Luqadda Wolfram waxay u isticmaashaa shaqada beddelka mawjada . Codsiga shaqadan ee dhammaan qaybaha raadraaca ayaa la fulin doonaa iyadoo la adeegsanayo shaqada . Halkan waxaa mudan in la xuso mid ka mid ah awoodaha Wolfram Mathematica - awoodda isticmaalka isbarbardhigga . Tusaalaha kore, ma jirto baahi loo qabo isbarbardhigga - mugga xogtu maaha mid weyn, laakiin marka la shaqeynayo xogta tijaabada ah ee ka kooban boqollaal kun oo raad raac ah, tani waa lagama maarmaan.
tracesCWD=Table[Map[Hilbert[#,0]&,Re[ContinuousWaveletTransform[traces[[i]]][[1]]][[{13,15,18}]]],{i,1,Length@traces}]; Kadib markaad codsato shaqada Xogta cusubi waxay u muuqataa mid u dhiganta wareegyada la doortay. Tusaalaha kore, soo noqnoqoshadani waa: 38Hz, 33Hz, 27Hz. Doorashada soo noqnoqoshada inta badan waxaa lagu fuliyaa iyadoo lagu saleynayo baaritaanka - waxay helayaan khariidado wax ku ool ah oo loogu talagalay isku-dhafka soo noqnoqda ee kala duwan waxayna doortaan midka ugu macluumaadka badan marka laga eego aragtida juqraafiyeed.
Haddii aad u baahan tahay inaad natiijooyinka la wadaagto asxaabta ama aad siiso macaamiisha, waxaad isticmaali kartaa shaqada SEGYExport[] ee xirmada GeologyIO
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];Saddex ka mid ah xabadahan (soo noqnoqoshada hoose, dhexe-dhexe, iyo qaybaha soo noqnoqda sare), isku darka RGB waxaa caadi ahaan loo isticmaalaa in lagu sawiro xogta. Qayb kasta waxaa loo qoondeeyey midabkeeda - casaan, cagaar, buluug. Wolfram Mathematica tan waxaa lagu samayn karaa iyadoo la isticmaalayo shaqada .
Natiijadu waa sawiro laga samayn karo tafsiir juquraafi. Dhexdhexaadiyeyaasha ku qoran qaybta ayaa suurtogal ka dhigaya in la kala saaro paleochannels, kuwaas oo u badan tahay inay noqdaan kaydad oo ay ku jiraan kayd shidaal. Raadinta iyo falanqaynta analoguesyada casriga ah ee nidaamka webiga noocan oo kale ah ayaa noo ogolaanaya inaan go'aamino qaybaha ugu rajo-gelinta badan ee dhexdhexaadiyeyaasha. Kanaalada laftooda waxaa lagu gartaa lakabyo qaro weyn oo ciid si fiican loo kala soocay waana kayd wanaagsan oo saliid ah. Meelaha ka baxsan cilladaha "shareerka" waxay la mid yihiin kaydadka casriga ah ee bannaanka daadku. Kaydka bannaanka daadku waxa inta badan metelaya dhagxaan dhoobo ah oo qodista aagaggani waxay noqon doontaa mid aan waxtar lahayn.
Qaybta RGB ee cube xogta. Bartamaha (wax yar dhanka bidix ee xarunta) waxaad ka raadsan kartaa webiga dhexdhexaadka ah.
Qaybta RGB ee cube xogta. Dhinaca bidix waxaad ka raadsan kartaa webiga dhexdhexaadka ah.
Xaaladaha qaarkood, tayada xogta seismic waxay u oggolaanaysaa sawiro cad cad. Tani waxay kuxirantahay habka shaqada goobta, qalabka loo isticmaalo algorithm dhimista qaylada. Xaaladahan oo kale, ma aha oo kaliya qaybo ka mid ah hababka wabiga, laakiin sidoo kale dhammaan webiyada paleo-ballaaran.
Isku darka RGB ee saddex qaybood oo ka mid ah cube xogta seismic ah (jeex-jeex). Qoto dheer qiyaastii 2 km.
Sawirka dayax-gacmeedka ee wabiga Volga ee u dhow Saratov
gunaanad
Wolfram Mathematica waxay falanqeyn kartaa xogta seismic waxayna xallin kartaa dhibaatooyinka la dabaqay ee la xiriira raadinta macdanta, iyo xirmada GeologyIO waxay ka dhigaysaa habkan mid ku habboon. Qaab dhismeedka xogta seismic waa sida iyadoo la adeegsanayo habab la dhisay si loo dardargeliyo xisaabinta (, ,β¦) waa mid aad u hufan oo kuu ogolaanaya inaad ka baaraandegto tiro badan oo xog ah. Ilaa xad ballaaran, tan waxaa fududeeyay astaamaha kaydinta xogta ee xirmada GeologyIO. Jid ahaan, xirmada waxaa loo isticmaali karaa ma aha oo kaliya goobta sahaminta seismic ee codsatay. Ku dhawaad ββnoocyada xogta la mid ah ayaa loo isticmaalaa in radar dhulka la geliyo iyo seismology.Haddii aad hayso soo jeedin ku saabsan sida loo hagaajiyo natiijada, kuwaas oo algorithms falanqaynta signal ka Wolfram Xisaabta arsenal waa lagu dabaqi karaa xogtan, ama haddii aad leedahay naqdin, fadlan ka tago a faalo.
Source: www.habr.com