Wolfram Mathematica ក្នុងភូគព្ភវិទ្យា

សូមអរគុណដល់អ្នកនិពន្ធប្លក់ លោក Anton Ekimenko សម្រាប់របាយការណ៍របស់គាត់។

សេចក្តីណែនាំ

កំណត់សម្គាល់នេះត្រូវបានសរសេរនៅក្រោយសន្និសីទ Wolfram សន្និសីទបច្ចេកវិទ្យារុស្ស៊ី និងមានសេចក្តីសង្ខេបនៃរបាយការណ៍ដែលខ្ញុំបានផ្តល់ឱ្យ។ ព្រឹត្តិការណ៍​នេះ​បាន​កើត​ឡើង​ក្នុង​ខែ​មិថុនា​ក្នុង​ទីក្រុង St. ដោយពិចារណាថាខ្ញុំធ្វើការប្លុកពីគេហទំព័រសន្និសីទ ខ្ញុំមិនអាចជួយចូលរួមព្រឹត្តិការណ៍នេះបានទេ។ នៅឆ្នាំ 2016 និង 2017 ខ្ញុំបានស្តាប់របាយការណ៍សន្និសីទ ហើយឆ្នាំនេះខ្ញុំបានធ្វើបទបង្ហាញមួយ។ ទីមួយ ប្រធានបទគួរឱ្យចាប់អារម្មណ៍មួយ (វាហាក់ដូចជាខ្ញុំ) បានបង្ហាញខ្លួន ដែលយើងកំពុងអភិវឌ្ឍជាមួយ Kirill Belovហើយទីពីរ បន្ទាប់ពីការសិក្សាយ៉ាងយូរនៃច្បាប់របស់សហព័ន្ធរុស្ស៊ី ទាក់ទងនឹងគោលនយោបាយដាក់ទណ្ឌកម្ម នៅឯសហគ្រាសដែលខ្ញុំធ្វើការ អាជ្ញាប័ណ្ណជាច្រើនបានបង្ហាញខ្លួន។ Wolfram គណិតវិទ្យា.

មុននឹងបន្តទៅប្រធានបទនៃសុន្ទរកថារបស់ខ្ញុំ ខ្ញុំចង់កត់សម្គាល់ការរៀបចំដ៏ល្អនៃព្រឹត្តិការណ៍នេះ។ ទំព័រទស្សនាសន្និសីទប្រើរូបភាពនៃវិហារ Kazan ។ វិហារនេះគឺជាកន្លែងទាក់ទាញដ៏សំខាន់មួយរបស់ St. Petersburg ហើយអាចមើលឃើញយ៉ាងច្បាស់ពីសាលដែលសន្និសីទបានធ្វើឡើង។

នៅច្រកចូលសាកលវិទ្យាល័យសេដ្ឋកិច្ច St. Petersburg អ្នកចូលរួមត្រូវបានជួបដោយជំនួយការក្នុងចំណោមនិស្សិត - ពួកគេមិនអនុញ្ញាតឱ្យពួកគេវង្វេង។ ក្នុងអំឡុងពេលចុះឈ្មោះ វត្ថុអនុស្សាវរីយ៍តូចៗត្រូវបានផ្តល់ឱ្យ (ប្រដាប់ក្មេងលេងមួយ - លោតភ្លឹបភ្លែតៗ ប៊ិច ស្ទីកឃ័រដែលមាននិមិត្តសញ្ញា Wolfram)។ អាហារថ្ងៃត្រង់ និងការសម្រាកកាហ្វេក៏ត្រូវបានបញ្ចូលក្នុងកាលវិភាគសន្និសីទផងដែរ។ ខ្ញុំបានកត់សម្គាល់រួចហើយអំពីកាហ្វេឆ្ងាញ់ និងនំនៅលើជញ្ជាំងរបស់ក្រុម - មេចុងភៅគឺអស្ចារ្យណាស់។ ជាមួយនឹងផ្នែកណែនាំនេះ ខ្ញុំចង់បញ្ជាក់ថា ព្រឹត្តិការណ៍ផ្ទាល់ ទម្រង់ និងទីតាំងរបស់វាកំពុងនាំមកនូវអារម្មណ៍វិជ្ជមានរួចហើយ។

របាយការណ៍ដែលត្រូវបានរៀបចំឡើងដោយខ្ញុំ និង Kirill Belov ត្រូវបានគេហៅថា "ការប្រើប្រាស់ Wolfram Mathematica ដើម្បីដោះស្រាយបញ្ហានៅក្នុងការអនុវត្តភូគព្ភសាស្ត្រ។ ការវិភាគវិសាលគមនៃទិន្នន័យរញ្ជួយដី ឬ "កន្លែងដែលទន្លេបុរាណបានរត់" ។ ខ្លឹមសារនៃរបាយការណ៍មានពីរផ្នែក៖ ទីមួយ ការប្រើប្រាស់ក្បួនដោះស្រាយដែលមាននៅក្នុង Wolfram គណិតវិទ្យា សម្រាប់ការវិភាគទិន្នន័យភូមិសាស្ត្រ និងទីពីរ នេះជារបៀបដាក់ទិន្នន័យភូមិសាស្ត្រទៅក្នុង Wolfram Mathematica ។

ការរុករករញ្ជួយដី

ដំបូង​អ្នក​ត្រូវ​ធ្វើ​ដំណើរ​កម្សាន្ត​ខ្លីៗ​ទៅ​ក្នុង​ភូមិសាស្ត្រ​។ ភូគព្ភសាស្ត្រ គឺជាវិទ្យាសាស្ត្រដែលសិក្សាពីលក្ខណៈរូបវន្តនៃថ្ម។ ជាការប្រសើរណាស់ ដោយសារថ្មមានលក្ខណៈសម្បត្តិផ្សេងៗគ្នា៖ អគ្គិសនី ម៉ាញេទិក ការបត់បែន មានវិធីសាស្រ្តដែលត្រូវគ្នានៃភូគព្ភសាស្ត្រ៖ ការទស្សន៍ទាយអគ្គិសនី ការទស្សន៍ទាយម៉ាញេទិក ការរំពឹងទុកនៃការរញ្ជួយដី... នៅក្នុងបរិបទនៃអត្ថបទនេះ យើងនឹងពិភាក្សាតែអំពីការរំពឹងទុកនៃការរញ្ជួយដីនៅក្នុងលម្អិតបន្ថែមទៀតប៉ុណ្ណោះ។ ការរុករករញ្ជួយដី គឺជាវិធីសាស្ត្រសំខាន់ក្នុងការស្វែងរកប្រេង និងឧស្ម័ន។ វិធីសាស្រ្តគឺផ្អែកលើការរំភើបនៃរំញ័រយឺត និងការកត់ត្រាជាបន្តបន្ទាប់នៃការឆ្លើយតបពីថ្មដែលបង្កើតតំបន់សិក្សា។ រំញ័ររំភើបនៅលើដី (ជាមួយឌីណាមិច ឬប្រភពរំញ័រដែលមិនផ្ទុះនៃរំញ័រយឺត) ឬនៅសមុទ្រ (ជាមួយកាំភ្លើងខ្យល់)។ រំញ័រ Elastic បន្តពូជតាមរយៈម៉ាស់ថ្ម ដែលត្រូវបានឆ្លុះបញ្ចាំង និងឆ្លុះបញ្ចាំងនៅព្រំដែននៃស្រទាប់ដែលមានលក្ខណៈសម្បត្តិខុសៗគ្នា។ រលកដែលឆ្លុះត្រឡប់មកផ្ទៃវិញ ហើយត្រូវបានកត់ត្រាដោយភូគព្ភសាស្ត្រនៅលើដី (ជាធម្មតាឧបករណ៍អេឡិចត្រូឌីណាមិកផ្អែកលើចលនារបស់មេដែកដែលផ្អាកនៅក្នុងឧបករណ៏) ឬអ៊ីដ្រូហ្វូននៅក្នុងសមុទ្រ (ផ្អែកលើឥទ្ធិពល piezoelectric) ។ នៅពេលនៃការមកដល់នៃរលក មនុស្សម្នាក់អាចវិនិច្ឆ័យជម្រៅនៃស្រទាប់ភូមិសាស្ត្រ។

ឧបករណ៍អូសនាវារញ្ជួយដី


កាំភ្លើងខ្យល់រំភើបរំញ័រយឺត


រលក​ឆ្លង​កាត់​ដុំ​ថ្ម ហើយ​ត្រូវ​បាន​កត់ត្រា​ដោយ​អ៊ីដ្រូហ្វូន


នាវាស្រាវជ្រាវភូមិសាស្ត្រ "Ivan Gubkin" នៅផែក្បែរស្ពាន Blagoveshchensky ក្នុងទីក្រុង St.

គំរូសញ្ញារញ្ជួយ

ថ្មមានលក្ខណៈសម្បត្តិរូបវ័ន្តខុសៗគ្នា។ សម្រាប់ការរុករករញ្ជួយដី លក្ខណៈសម្បត្តិបត់បែនមានសារៈសំខាន់ជាចម្បង - ល្បឿននៃការសាយភាយនៃរំញ័រយឺត និងដង់ស៊ីតេ។ ប្រសិនបើស្រទាប់ពីរមានលក្ខណៈសម្បត្តិដូចគ្នា ឬស្រដៀងគ្នា នោះរលក "នឹងមិនកត់សម្គាល់" ព្រំដែនរវាងពួកវាទេ។ ប្រសិនបើល្បឿនរលកក្នុងស្រទាប់ខុសគ្នា នោះការឆ្លុះបញ្ចាំងនឹងកើតឡើងនៅព្រំដែននៃស្រទាប់។ ភាពខុសគ្នានៃលក្ខណៈសម្បត្តិកាន់តែច្រើន ការឆ្លុះបញ្ចាំងកាន់តែខ្លាំង។ អាំងតង់ស៊ីតេរបស់វានឹងត្រូវបានកំណត់ដោយមេគុណឆ្លុះបញ្ចាំង (rc)៖

ដែល ρ គឺជាដង់ស៊ីតេថ្ម ν គឺជាល្បឿនរលក 1 និង 2 បង្ហាញពីស្រទាប់ខាងលើ និងខាងក្រោម។

គំរូសញ្ញារញ្ជួយដីដ៏សាមញ្ញបំផុត និងប្រើញឹកញាប់បំផុតមួយគឺគំរូ convolution នៅពេលដែលដានរញ្ជួយដីដែលបានកត់ត្រាត្រូវបានតំណាងថាជាលទ្ធផលនៃការ convolution នៃលំដាប់នៃមេគុណឆ្លុះបញ្ចាំងជាមួយនឹងជីពចរស៊ើបអង្កេតមួយ:

កន្លែងណា s(t) - ដានរញ្ជួយដី, ឧ។ អ្វី​គ្រប់​យ៉ាង​ដែល​ត្រូវ​បាន​ថត​ដោយ​អ៊ីដ្រូ​ហ្វូន ឬ​ភូមិសាស្ត្រ​ក្នុង​អំឡុង​ពេល​ថត​ថេរ w(t) - សញ្ញាដែលបង្កើតដោយកាំភ្លើងខ្យល់ n(t) - សំលេងរំខានចៃដន្យ។

ចូរយើងគណនាដានរញ្ជួយដីសំយោគជាឧទាហរណ៍មួយ។ យើងនឹងប្រើជីពចរ Ricker ដែលត្រូវបានគេប្រើយ៉ាងទូលំទូលាយក្នុងការរុករករញ្ជួយដី ជាសញ្ញាដំបូង។

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 ហើយបរិមាណនៃវត្ថុធាតុដើម "ដើម" អាចមានរាប់រយជីហ្គាបៃ។

របៀបធ្វើការជាមួយពួកគេ។ Wolfram គណិតវិទ្យា?

កញ្ចប់ ភូគព្ភសាស្ត្រ

ការអភិវឌ្ឍន៍កញ្ចប់បានចាប់ផ្តើម សំណួរ នៅលើជញ្ជាំង VK នៃក្រុមគាំទ្រនិយាយភាសារុស្ស៊ី។ អរគុណចំពោះការឆ្លើយតបរបស់សហគមន៍ ដំណោះស្រាយមួយត្រូវបានរកឃើញយ៉ាងឆាប់រហ័ស។ ហើយ​ជា​លទ្ធផល វា​បាន​រីក​ចម្រើន​ទៅ​ជា​ការ​អភិវឌ្ឍ​ដ៏​ធ្ងន់ធ្ងរ។ ដែលត្រូវគ្នា។ ប្រកាសជញ្ជាំងសហគមន៍ Wolfram វាត្រូវបានសម្គាល់ដោយអ្នកសម្របសម្រួល។ បច្ចុប្បន្ន កញ្ចប់គាំទ្រការធ្វើការជាមួយប្រភេទទិន្នន័យខាងក្រោម ដែលត្រូវបានប្រើយ៉ាងសកម្មនៅក្នុងឧស្សាហកម្មភូគព្ភសាស្ត្រ៖

1. ការនាំចូលទិន្នន័យផែនទីក្នុងទម្រង់ ZMAP និង IRAP
2. ការនាំចូលការវាស់វែងក្នុងទម្រង់ LAS អណ្តូង
3. ការបញ្ចូល និងលទ្ធផលនៃទ្រង់ទ្រាយឯកសាររញ្ជួយ ស៊ីជី

ដើម្បីដំឡើងកញ្ចប់ អ្នកត្រូវតែធ្វើតាមការណែនាំនៅលើទំព័រទាញយកនៃកញ្ចប់ដែលបានជួបប្រជុំគ្នា i.e. ប្រតិបត្តិកូដខាងក្រោមនៅក្នុងណាមួយ។ សៀវភៅកត់ត្រាគណិតវិទ្យា:

If[PacletInformation["GeologyIO"] === {}, PacletInstall[URLDownload[
    "https://wolfr.am/FiQ5oFih", 
    FileNameJoin[{CreateDirectory[], "GeologyIO-0.2.2.paclet"}]
]]]

បន្ទាប់ពីនោះកញ្ចប់នឹងត្រូវបានដំឡើងនៅក្នុងថតលំនាំដើម ផ្លូវទៅកាន់ដែលអាចទទួលបានដូចខាងក្រោម៖

FileNameJoin[{$UserBasePacletsDirectory, "Repository"}]

ជាឧទាហរណ៍ យើងនឹងបង្ហាញពីសមត្ថភាពចម្បងនៃកញ្ចប់។ ការហៅទូរសព្ទត្រូវបានធ្វើឡើងជាប្រពៃណីសម្រាប់កញ្ចប់ជាភាសា Wolfram៖

Get["GeologyIO`"]

កញ្ចប់ត្រូវបានបង្កើតឡើងដោយប្រើ Wolfram Workbench. នេះអនុញ្ញាតឱ្យអ្នកអមដំណើរមុខងារសំខាន់នៃកញ្ចប់ជាមួយនឹងឯកសារ ដែលនៅក្នុងលក្ខខណ្ឌនៃទម្រង់បទបង្ហាញមិនខុសគ្នាពីឯកសាររបស់ Wolfram Mathematica ខ្លួនវាផ្ទាល់ និងដើម្បីផ្តល់កញ្ចប់ជាមួយនឹងឯកសារសាកល្បងសម្រាប់អ្នកស្គាល់គ្នាដំបូង។

ជាពិសេសឯកសារបែបនេះគឺជាឯកសារ "Marmousi.segy" - នេះគឺជាគំរូសំយោគនៃផ្នែកភូមិសាស្ត្រដែលត្រូវបានបង្កើតឡើងដោយវិទ្យាស្ថានប្រេងបារាំង។ ដោយប្រើគំរូនេះ អ្នកអភិវឌ្ឍន៍សាកល្បងក្បួនដោះស្រាយផ្ទាល់ខ្លួនរបស់ពួកគេសម្រាប់គំរូវាលរលក ដំណើរការទិន្នន័យ ការបញ្ច្រាសដានរញ្ជួយ ជាដើម។ ម៉ូដែល Marmousi ខ្លួនវាត្រូវបានរក្សាទុកនៅក្នុងឃ្លាំងពីកន្លែងដែលកញ្ចប់ខ្លួនវាត្រូវបានទាញយក។ ដើម្បីទទួលបានឯកសារ សូមដំណើរការកូដខាងក្រោម៖

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 អនុញ្ញាតឱ្យអ្នកប្រើមុខងារដំណើរការសញ្ញាដែលមានស្រាប់សម្រាប់ទិន្នន័យពិសោធន៍។ ដោយសារដានរញ្ជួយដីនីមួយៗតំណាងឱ្យស៊េរីពេលវេលា ឧបករណ៍សំខាន់មួយសម្រាប់សិក្សាពួកវាគឺការវិភាគវិសាលគម។ ក្នុងចំណោមតម្រូវការជាមុនសម្រាប់ការវិភាគសមាសភាពប្រេកង់នៃទិន្នន័យរញ្ជួយដី យើងអាចដាក់ឈ្មោះឧទាហរណ៍ដូចខាងក្រោមៈ

1. ប្រភេទផ្សេងគ្នានៃរលកត្រូវបានកំណត់លក្ខណៈដោយសមាសភាពប្រេកង់ផ្សេងគ្នា។ នេះអនុញ្ញាតឱ្យអ្នករំលេចរលកមានប្រយោជន៍ និងទប់ស្កាត់រលកជ្រៀតជ្រែក។
2. លក្ខណៈសម្បត្តិថ្មដូចជា porosity និង saturation អាចប៉ះពាល់ដល់សមាសភាពប្រេកង់។ នេះធ្វើឱ្យវាអាចកំណត់អត្តសញ្ញាណថ្មដែលមានលក្ខណៈសម្បត្តិល្អបំផុត។
3. ស្រទាប់ដែលមានកម្រាស់ខុសៗគ្នាបណ្តាលឱ្យមានភាពមិនប្រក្រតីនៅក្នុងជួរប្រេកង់ផ្សេងៗគ្នា។

ចំណុចទីបីគឺជាចំណុចសំខាន់នៅក្នុងបរិបទនៃអត្ថបទនេះ។ ខាងក្រោមនេះគឺជាបំណែកកូដសម្រាប់គណនាដានរញ្ជួយដីនៅក្នុងករណីនៃស្រទាប់ដែលមានកម្រាស់ខុសគ្នា - គំរូក្រូចឆ្មារ។ គំរូនេះត្រូវបានសិក្សាជាប្រពៃណីនៅក្នុងការរុករករញ្ជួយដី ដើម្បីវិភាគឥទ្ធិពលជ្រៀតជ្រែក នៅពេលដែលរលកដែលឆ្លុះបញ្ចាំងពីស្រទាប់ជាច្រើនត្រូវបានដាក់លើគ្នាទៅវិញទៅមក។

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]]

គំរូនៃការបង្កើតស្នាមប្រេះ


ល្បឿនរលកនៅខាងក្នុងក្រូចឆ្មារគឺ 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]

ដានរញ្ជួយដីសម្រាប់គំរូក្រូចឆ្មារ


លំដាប់​នៃ​ដាន​រញ្ជួយ​ដី​ដែល​បង្ហាញ​ក្នុង​រូប​នេះ​ត្រូវ​បាន​គេ​ហៅ​ថា​ផ្នែក​រញ្ជួយ​ដី។ ដូចដែលអ្នកអាចឃើញការបកស្រាយរបស់វាក៏អាចត្រូវបានអនុវត្តនៅលើកម្រិតវិចារណញាណមួយចាប់តាំងពីធរណីមាត្រនៃរលកដែលឆ្លុះបញ្ចាំងយ៉ាងច្បាស់ត្រូវគ្នាទៅនឹងគំរូដែលបានបញ្ជាក់ពីមុន។ ប្រសិនបើអ្នកវិភាគដានឱ្យបានលម្អិតអ្នកនឹងសម្គាល់ឃើញថាដានពីលេខ 1 ដល់ប្រហែល 30 មិនខុសគ្នាទេ - ការឆ្លុះបញ្ចាំងពីដំបូលនៃការបង្កើតនិងពីខាងក្រោមមិនត្រួតលើគ្នា។ ចាប់ផ្តើមពីដានទី 31 ការឆ្លុះបញ្ចាំងចាប់ផ្តើមជ្រៀតជ្រែក។ ហើយទោះបីជានៅក្នុងគំរូក៏ដោយមេគុណនៃការឆ្លុះបញ្ចាំងមិនផ្លាស់ប្តូរផ្ដេកទេ - ដានរញ្ជួយដីផ្លាស់ប្តូរអាំងតង់ស៊ីតេរបស់វានៅពេលដែលកម្រាស់នៃការបង្កើតផ្លាស់ប្តូរ។

ចូរយើងពិចារណាពីទំហំនៃការឆ្លុះបញ្ចាំងពីព្រំដែនខាងលើនៃការបង្កើត។ ចាប់ផ្តើមពីផ្លូវទី 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


បណ្តាញនៃអណ្តូងទឹកនៅវាលមួយ។


ថ្មដែលផ្ទុកប្រេងនៃស៊ីបេរីខាងលិចកើតឡើងក្នុងជម្រៅធំទូលាយ - ពី 1 គីឡូម៉ែត្រទៅ 5 គីឡូម៉ែត្រ។ បរិមាណសំខាន់នៃថ្មដែលមានប្រេងត្រូវបានបង្កើតឡើងនៅសម័យ Jurassic និង Cretaceous ។ សម័យ Jurassic ប្រហែលជាត្រូវបានគេស្គាល់ដោយមនុស្សជាច្រើនពីខ្សែភាពយន្តដែលមានឈ្មោះដូចគ្នា។ អាកាសធាតុ Jurassic មានភាពខុសគ្នាយ៉ាងខ្លាំងពីសម័យទំនើប។ សព្វវចនាធិប្បាយ Britannica មានស៊េរីនៃ paleomaps ដែលកំណត់លក្ខណៈនៃយុគសម័យ helogical នីមួយៗ។

តឹងតែងបច្ចុប្បន្ន

រយៈពេល Jurassic


សូមចំណាំថានៅក្នុងសម័យ Jurassic ទឹកដីនៃភាគខាងលិចស៊ីបេរីគឺជាឆ្នេរសមុទ្រ (ដីឆ្លងកាត់ដោយទន្លេនិងសមុទ្ររាក់) ។ ដោយសារអាកាសធាតុមានផាសុកភាព យើងអាចសន្មត់ថាទេសភាពធម្មតានៅសម័យនោះមើលទៅដូចនេះ៖

Jurassic ស៊ីបេរី


នៅក្នុងរូបភាពនេះ អ្វីដែលសំខាន់សម្រាប់យើងគឺមិនមែនសត្វ និងសត្វស្លាបច្រើននោះទេ ប៉ុន្តែជារូបភាពនៃទន្លេនៅផ្ទៃខាងក្រោយ។ ទន្លេគឺជាវត្ថុភូគព្ភសាស្ត្រដូចគ្នាដែលយើងបានឈប់នៅមុននេះ។ ការពិតគឺថាសកម្មភាពនៃទន្លេអនុញ្ញាតឱ្យថ្មភក់ដែលបានតម្រៀបយ៉ាងល្អដើម្បីកកកុញដែលបន្ទាប់មកនឹងក្លាយជាអាងស្តុកទឹកសម្រាប់ប្រេង។ អាងស្តុកទឹកទាំងនេះអាចមានរាងស្មុគ្រស្មាញ និងចម្លែក (ដូចជាគ្រែទន្លេ) ហើយពួកវាមានកំរាស់អថេរ - នៅជិតច្រាំងទន្លេ កំរាស់គឺតូច ប៉ុន្តែខិតទៅជិតកណ្តាលនៃឆានែល ឬនៅតំបន់ផ្លូវវាងវាកើនឡើង។ ដូច្នេះ ទន្លេដែលបង្កើតនៅក្នុង Jurassic ឥឡូវនេះមានជម្រៅប្រហែល XNUMX គីឡូម៉ែត្រ ហើយជាវត្ថុនៃការស្វែងរកអាងស្តុកប្រេង។

ទិន្នន័យពិសោធន៍។ ដំណើរការនិងការមើលឃើញ

អនុញ្ញាតឱ្យយើងធ្វើការកក់ទុកភ្លាមៗទាក់ទងនឹងសម្ភារៈរញ្ជួយដីដែលបានបង្ហាញនៅក្នុងអត្ថបទ - ដោយសារតែការពិតដែលថាបរិមាណនៃទិន្នន័យដែលបានប្រើសម្រាប់ការវិភាគគឺមានសារៈសំខាន់ - មានតែបំណែកនៃដាននៃការរញ្ជួយដីដើមប៉ុណ្ណោះដែលត្រូវបានរួមបញ្ចូលនៅក្នុងអត្ថបទនៃអត្ថបទ។ នេះនឹងអនុញ្ញាតឱ្យនរណាម្នាក់បង្កើតឡើងវិញនូវការគណនាខាងលើ។

នៅពេលធ្វើការជាមួយទិន្នន័យរញ្ជួយដី ភូគព្ភវិទូជាធម្មតាប្រើកម្មវិធីឯកទេស (មានអ្នកដឹកនាំឧស្សាហកម្មជាច្រើនដែលការអភិវឌ្ឍន៍ត្រូវបានប្រើប្រាស់យ៉ាងសកម្ម ឧទាហរណ៍ Petrel ឬ Paradigm) ដែលអនុញ្ញាតឱ្យអ្នកវិភាគប្រភេទទិន្នន័យផ្សេងៗគ្នា និងមានចំណុចប្រទាក់ក្រាហ្វិកងាយស្រួល។ ទោះបីជាមានភាពងាយស្រួលទាំងអស់ក៏ដោយក៏ប្រភេទនៃកម្មវិធីទាំងនេះក៏មានគុណវិបត្តិរបស់ពួកគេផងដែរ - ឧទាហរណ៍ការអនុវត្តក្បួនដោះស្រាយទំនើបនៅក្នុងកំណែដែលមានស្ថេរភាពត្រូវការពេលវេលាច្រើនហើយលទ្ធភាពសម្រាប់ការគណនាដោយស្វ័យប្រវត្តិជាធម្មតាមានកម្រិត។ ក្នុងស្ថានភាពបែបនេះ វាងាយស្រួលប្រើក្នុងការប្រើប្រាស់ប្រព័ន្ធគណិតវិទ្យាកុំព្យូទ័រ និងភាសាសរសេរកម្មវិធីកម្រិតខ្ពស់ ដែលអនុញ្ញាតឱ្យប្រើមូលដ្ឋានក្បួនដោះស្រាយដ៏ធំទូលាយ ហើយក្នុងពេលជាមួយគ្នានេះ អនុវត្តទម្លាប់ច្រើន។ នេះគឺជាគោលការណ៍ដែលប្រើដើម្បីធ្វើការជាមួយទិន្នន័យរញ្ជួយដីនៅក្នុង Wolfram Mathematica ។ វាមិនសមស្របទេក្នុងការសរសេរមុខងារសម្បូរបែបសម្រាប់ការងារអន្តរកម្មជាមួយទិន្នន័យ - វាមានសារៈសំខាន់ជាងក្នុងការធានាការផ្ទុកពីទម្រង់ដែលទទួលយកជាទូទៅ អនុវត្តក្បួនដោះស្រាយដែលចង់បានទៅពួកវា និងផ្ទុកពួកវាត្រឡប់ទៅទម្រង់ខាងក្រៅវិញ។

អនុវត្តតាមគ្រោងការណ៍ដែលបានស្នើឡើង យើងនឹងផ្ទុកទិន្នន័យរញ្ជួយដីដើម ហើយបង្ហាញពួកវានៅក្នុង Wolfram គណិតវិទ្យា:

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]

ទិន្នន័យ​ដែល​បាន​ទាញយក​និង​នាំចូល​តាម​វិធី​នេះ​គឺ​ជា​ផ្លូវ​ដែល​បាន​កត់ត្រា​នៅលើ​ផ្ទៃដី​ដែល​មាន​ទំហំ​ ១០ គុណនឹង ៥ គីឡូម៉ែត្រ។ ប្រសិនបើទិន្នន័យត្រូវបានទទួលដោយប្រើវិធីសាស្ត្រស្ទង់មតិរញ្ជួយដីបីវិមាត្រ (រលកត្រូវបានកត់ត្រាមិននៅតាមបណ្តោយទម្រង់ភូមិសាស្ត្រនីមួយៗ ប៉ុន្តែនៅលើផ្ទៃដីទាំងមូលក្នុងពេលដំណាលគ្នា) វាអាចទទួលបានគូបទិន្នន័យរញ្ជួយដី។ ទាំងនេះគឺជាវត្ថុបីវិមាត្រ ផ្នែកបញ្ឈរ និងផ្ដេក ដែលអនុញ្ញាតឱ្យមានការសិក្សាលម្អិតអំពីបរិស្ថានភូមិសាស្ត្រ។ នៅក្នុងឧទាហរណ៍ដែលបានពិចារណា យើងកំពុងដោះស្រាយជាមួយទិន្នន័យបីវិមាត្រ។ យើងអាចទទួលបានព័ត៌មានមួយចំនួនពីបឋមកថាអត្ថបទដូចនេះ

StringPartition[seismic3DSEGY["textheader"], 80] // TableForm

C 1 នេះគឺជាឯកសារសាកល្បងសម្រាប់ការធ្វើតេស្តកញ្ចប់ភូមិសាស្ត្រ
C 2
C 3
C 4
C 5 DATE ឈ្មោះអ្នកប្រើប្រាស់៖ WOLFRAM USER
C 6 ឈ្មោះស្ទង់មតិ៖ កន្លែងណាមួយនៅស៊ីបេរី
C 7 ប្រភេទឯកសារ 3D SEISMIC volume
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]]

រូបភាព XNUMXD នៃគូបទិន្នន័យរញ្ជួយ។ (អ័ក្សបញ្ឈរ - ជម្រៅ)


ប្រសិនបើលក្ខណៈភូមិសាស្ត្រនៃចំណាប់អារម្មណ៍បង្កើតភាពខុសប្រក្រតីនៃការរញ្ជួយដីខ្លាំងនោះ ឧបករណ៍មើលឃើញដែលមានតម្លាភាពអាចប្រើប្រាស់បាន។ តំបន់ "មិនសំខាន់" នៃការថតអាចត្រូវបានធ្វើឱ្យមើលមិនឃើញដោយបន្សល់ទុកតែភាពមិនប្រក្រតីដែលអាចមើលឃើញ។ នៅក្នុង Wolfram Mathematica នេះអាចត្រូវបានធ្វើដោយប្រើ ភាពស្រអាប់[] и Raster3D[].

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[]

ដូចនៅក្នុងឧទាហរណ៍សំយោគ នៅផ្នែកនៃគូបដើម គេអាចកំណត់ព្រំដែនភូមិសាស្ត្រមួយចំនួន (ស្រទាប់) ជាមួយនឹងភាពធូរស្រាលអថេរ។

ឧបករណ៍សំខាន់សម្រាប់ការវិភាគវិសាលគមគឺការបំប្លែង Fourier ។ ដោយមានជំនួយរបស់វា អ្នកអាចវាយតម្លៃវិសាលគមប្រេកង់អំព្លីទីត នៃដាននីមួយៗ ឬក្រុមនៃដាន។ ទោះយ៉ាងណាក៏ដោយបន្ទាប់ពីការផ្ទេរទិន្នន័យទៅដែនប្រេកង់ព័ត៌មានត្រូវបានបាត់បង់អំពីពេលវេលា (អាននៅជម្រៅអ្វី) ប្រេកង់ផ្លាស់ប្តូរ។ ដើម្បីអាចធ្វើមូលដ្ឋានីយកម្មការផ្លាស់ប្តូរសញ្ញានៅលើអ័ក្សពេលវេលា (ជម្រៅ) ការបំប្លែង Fourier ដែលមានបង្អួច និងការបំបែករលកធាតុអាកាសត្រូវបានប្រើ។ អត្ថបទនេះប្រើការបំបែករលក។ បច្ចេកវិទ្យាការវិភាគរលកធាតុអាកាសបានចាប់ផ្តើមប្រើប្រាស់យ៉ាងសកម្មក្នុងការរុករករញ្ជួយដីក្នុងទសវត្សរ៍ទី 90 ។ អត្ថប្រយោជន៍លើការបំប្លែង 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]

ការបំបែកដានចូលទៅក្នុងសមាសធាតុ


ដើម្បីវាយតម្លៃពីរបៀបដែលថាមពលឆ្លុះបញ្ចាំងត្រូវបានចែកចាយនៅពេលវេលាមកដល់រលកផ្សេងៗគ្នា មាត្រដ្ឋាន (អាណាឡូកទៅនឹងវិសាលគម) ត្រូវបានប្រើ។ តាមក្បួនមួយនៅក្នុងការអនុវត្តវាមិនចាំបាច់វិភាគសមាសធាតុទាំងអស់ទេ។ ជាធម្មតា សមាសធាតុប្រេកង់ទាប កណ្តាល និងខ្ពស់ត្រូវបានជ្រើសរើស។

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]

មាត្រដ្ឋាន។ លទ្ធផលមុខងារ WaveletScalogram[]


ភាសា Wolfram ប្រើមុខងារសម្រាប់ការបំប្លែងរលក ContinuousWaveletTransform[]. ហើយការអនុវត្តមុខងារនេះទៅសំណុំដានទាំងមូលនឹងត្រូវបានអនុវត្តដោយប្រើមុខងារ តារាង[]. នៅទីនេះវាមានតម្លៃកត់សម្គាល់ចំណុចខ្លាំងមួយនៃ Wolfram Mathematica - សមត្ថភាពក្នុងការប្រើភាពស្របគ្នា។ តារាងប៉ារ៉ាឡែល[]. ក្នុងឧទាហរណ៍ខាងលើ មិនចាំបាច់មានការប៉ារ៉ាឡែលទេ - បរិមាណទិន្នន័យមិនធំទេ ប៉ុន្តែនៅពេលធ្វើការជាមួយសំណុំទិន្នន័យពិសោធន៍ដែលមានដានរាប់រយរាប់ពាន់ នេះគឺជាការចាំបាច់។

tracesCWD=Table[Map[Hilbert[#,0]&,Re[ContinuousWaveletTransform[traces[[i]]][[1]]][[{13,15,18}]]],{i,1,Length@traces}]; 

បន្ទាប់ពីអនុវត្តមុខងារ ContinuousWaveletTransform[] សំណុំទិន្នន័យថ្មីលេចឡើងដែលត្រូវគ្នាទៅនឹងប្រេកង់ដែលបានជ្រើសរើស។ នៅក្នុងឧទាហរណ៍ខាងលើ ប្រេកង់ទាំងនេះគឺ: 38Hz, 33Hz, 27Hz ។ ជម្រើសនៃប្រេកង់ត្រូវបានអនុវត្តជាញឹកញាប់បំផុតនៅលើមូលដ្ឋាននៃការធ្វើតេស្ត - ពួកគេទទួលបានផែនទីដែលមានប្រសិទ្ធភាពសម្រាប់ការផ្សំប្រេកង់ផ្សេងៗគ្នាហើយជ្រើសរើសព័ត៌មានដែលមានព័ត៌មានច្រើនបំផុតពីទស្សនៈរបស់អ្នកភូគព្ភវិទូ។

ប្រសិនបើអ្នកត្រូវការចែករំលែកលទ្ធផលជាមួយសហសេវិក ឬផ្តល់វាដល់អតិថិជន អ្នកអាចប្រើមុខងារ SEGYExport[] នៃកញ្ចប់ 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];

ជាមួយនឹងគូបទាំងបីនេះ (ប្រេកង់ទាប ពាក់កណ្តាលប្រេកង់ និងសមាសធាតុប្រេកង់ខ្ពស់) ការលាយបញ្ចូលគ្នា RGB ត្រូវបានប្រើជាធម្មតា ដើម្បីមើលឃើញទិន្នន័យរួមគ្នា។ សមាសធាតុនីមួយៗត្រូវបានកំណត់ពណ៌ផ្ទាល់ខ្លួនរបស់វា - ក្រហមបៃតងខៀវ។ នៅក្នុង Wolfram Mathematica នេះអាចត្រូវបានធ្វើដោយប្រើមុខងារ ColorCombine[].

លទ្ធផលគឺជារូបភាពដែលការបកស្រាយភូគព្ភសាស្ត្រអាចត្រូវបានធ្វើឡើង។ មគ្គុទ្ទេសក៍ដែលត្រូវបានកត់ត្រានៅលើផ្នែកនេះធ្វើឱ្យវាអាចធ្វើទៅបានដើម្បីកំណត់ paleochannels ដែលទំនងជាអាងស្តុកទឹក និងមានស្តុកប្រេង។ ការស្វែងរកនិងការវិភាគនៃ analogues ទំនើបនៃប្រព័ន្ធទន្លេបែបនេះអនុញ្ញាតឱ្យយើងកំណត់ផ្នែកដែលរំពឹងទុកបំផុតនៃ meanders ។ ឆានែលខ្លួនឯងត្រូវបានកំណត់លក្ខណៈដោយស្រទាប់ក្រាស់នៃថ្មភក់ដែលបានតម្រៀបយ៉ាងល្អហើយជាអាងស្តុកទឹកដ៏ល្អសម្រាប់ប្រេង។ តំបន់នៅខាងក្រៅភាពមិនធម្មតានៃ "ចរ" គឺស្រដៀងទៅនឹងប្រាក់បញ្ញើទឹកជំនន់ទំនើប។ ស្រទាប់ទំនាបទឹកជំនន់ត្រូវបានតំណាងជាចម្បងដោយថ្ម clayey ហើយការខួងចូលទៅក្នុងតំបន់ទាំងនេះនឹងមិនមានប្រសិទ្ធភាពទេ។

បំណែក RGB នៃគូបទិន្នន័យ។ នៅកណ្តាល (បន្តិចទៅខាងឆ្វេងនៃកណ្តាល) អ្នកអាចតាមដានទន្លេដែលវង្វេង។

បំណែក RGB នៃគូបទិន្នន័យ។ នៅផ្នែកខាងឆ្វេងអ្នកអាចតាមដានទន្លេដែលវង្វេង។


ក្នុងករណីខ្លះ គុណភាពនៃទិន្នន័យរញ្ជួយអនុញ្ញាតឱ្យមានរូបភាពកាន់តែច្បាស់។ នេះអាស្រ័យលើវិធីសាស្រ្តការងារវាល ឧបករណ៍ដែលប្រើដោយក្បួនដោះស្រាយកាត់បន្ថយសំលេងរំខាន។ ក្នុងករណីបែបនេះ មិនត្រឹមតែបំណែកនៃប្រព័ន្ធទន្លេប៉ុណ្ណោះដែលអាចមើលឃើញប៉ុណ្ណោះទេ ប៉ុន្តែក៏មានទន្លេប៉ាលេអូដែលលាតសន្ធឹងផងដែរ។

ការលាយ RGB នៃសមាសធាតុបីនៃដុំទិន្នន័យរញ្ជួយដី (ចំណិតផ្តេក)។ ជម្រៅប្រហែល 2 គីឡូម៉ែត្រ។

រូបភាពផ្កាយរណបនៃទន្លេ Volga នៅជិត Saratov

សេចក្តីសន្និដ្ឋាន

Wolfram Mathematica អនុញ្ញាតឱ្យអ្នកវិភាគទិន្នន័យរញ្ជួយដី និងដោះស្រាយបញ្ហាដែលបានអនុវត្តទាក់ទងនឹងការរុករករ៉ែ ហើយកញ្ចប់ GeologyIO ធ្វើឱ្យដំណើរការនេះកាន់តែងាយស្រួល។ រចនាសម្ព័ននៃទិន្នន័យរញ្ជួយដីគឺបែបនោះដោយប្រើវិធីសាស្រ្តដែលភ្ជាប់មកជាមួយដើម្បីបង្កើនល្បឿននៃការគណនា (តារាងប៉ារ៉ាឡែល[], ParallelDo[],…) មានប្រសិទ្ធភាពខ្ពស់ និងអនុញ្ញាតឱ្យអ្នកដំណើរការទិន្នន័យយ៉ាងច្រើន។ ក្នុងកម្រិតធំ នេះត្រូវបានសម្របសម្រួលដោយលក្ខណៈពិសេសនៃការផ្ទុកទិន្នន័យនៃកញ្ចប់ GeologyIO ។ ដោយវិធីនេះកញ្ចប់អាចត្រូវបានប្រើមិនត្រឹមតែនៅក្នុងវិស័យនៃការរុករករញ្ជួយដីដែលបានអនុវត្តនោះទេ។ ទិន្នន័យស្ទើរតែដូចគ្នាត្រូវបានប្រើប្រាស់ក្នុងរ៉ាដាដែលជ្រាបចូលដី និងរញ្ជួយដី។ ប្រសិនបើអ្នកមានការណែនាំអំពីរបៀបកែលម្អលទ្ធផល ដែលក្បួនដោះស្រាយការវិភាគសញ្ញាពីឃ្លាំងអាវុធ Wolfram Mathematica អាចអនុវត្តបានចំពោះទិន្នន័យបែបនេះ ឬប្រសិនបើអ្នកមានមតិរិះគន់ណាមួយ សូម ទុកមតិយោបល់។

ប្រភព: www.habr.com