Ushbu tadqiqot metan gidrati o‘sishi tajribalaridan hisoblangan ko‘rinma yoki samarali kinetik parametrni sodda va izohlanadigan ma’lumot modeli bilan bashorat qilish mumkin yoki yo‘qligini o‘rganadi. Tadqiqotchi avval chop etilgan eksperimental tadqiqotdan olingan 30 ta kuzatuvdan foydalangan; bosim, gidrat muvozanat harorati va ortiqcha sovishni ko‘p o‘zgaruvchili chiziqli regressiya modelining kirishlari sifatida tanlagan. Yetti xil o‘zgaruvchi kombinatsiyasi bir xil 10 qatlamli kross-validatsiya usuli bilan taqqoslangan.

Eng muvaffaqiyatli model bo‘lgan Model 7 bosim, muvozanat harorati va ortiqcha sovishni birgalikda ishlatgan. Ushbu model uchun kross-validatsiyalangan R² qiymati 0,9682, ildiz o‘rtacha kvadratik xato 3092,5 W/m²K va o‘rtacha mutlaq xato 2429,1 W/m²K deb bildirilgan. Xuddi shu ma’lumot va validatsiya usuli bilan baholangan Arrhenius modeli esa R² = 0,9540, RMSE = 3712,0 W/m²K va MAE = 2967,0 W/m²K natijalarini bergan. Regressiya modeli qo‘shimcha operatsion o‘zgaruvchilarni o‘z ichiga olgani uchun o‘rganilgan ma’lumotlar to‘plamida pastroq bashorat xatosini ta’minlagan.

Biroq bashorat qilinayotgan kattalik molekulyar miqyosda bevosita o‘lchangan sof reaksiya konstantasi emas. Parametr eksperimental interfeys o‘sish tezligi, gidrat hosil bo‘lish issiqligi, gidrat zichligi va issiqlik uzatish koeffitsientidan hosil qilingan. Shu sabab u termodinamik harakatlantiruvchi kuch, issiqlik va massa uzatish, interfeys hodisalari hamda tajriba qurilmasining umumiy ta’sirlarini o‘z ichiga olgan ko‘rinma parametr sifatida baholanishi kerak. Model ham fizik gidrat o‘sish qonuni emas, faqat o‘xshash sharoitdagi tajribalar uchun ma’lumotlar to‘plamiga xos empirik yondashuvni taqdim etadi.

Turkiya nuqtai nazaridan baholash: Tadqiqot Turkiyada tabiiy gaz uzatish quvurlari, suv osti quvur tizimlari, oqim ishonchliligi, yuqori bosimli jarayonlar va energetika muhandisligi sohalarida kichik eksperimental ma’lumotlar to‘plamlarini shaffof modellar bilan baholash uchun foydali usul namunasini beradi. Yondashuvni Turkiyada qo‘llash uchun mahalliy gaz tarkiblari, quvur va reaktor geometriyalari, oqim sharoitlari, sho‘rlik, ingibitorlar, aralashtirish darajasi va issiqlik-massa uzatish sharoitlarini qamrab oladigan mustaqil eksperimental ma’lumotlar to‘planishi kerak. Model turli tizimlarga ko‘chirilishidan oldin tashqi ma’lumotlar bilan tekshirilishi, takroriy yoki ichma-ich kross-validatsiya qo‘llanishi va fizik o‘zgaruvchilar bilan kengaytirilishi lozim. Ushbu tadqiqotdan Turkiyadagi istalgan quvur liniyasining gidrat tiqilish xavfi, xavfsiz ish chegarasi yoki ingibitor ehtiyoji bevosita hisoblab chiqilmaydi.

Batafsil tushuntirish

Tadqiqot ko‘rib chiqayotgan muammo nima?

Yuqori bosim va past harorat sharoitida suv molekulalari vodorod bog‘lari orqali qafasga o‘xshash kristall tuzilmalar hosil qilib, metan kabi kichik gaz molekulalarini shu qafaslar ichida tutib qolishi mumkin. Hosil bo‘lgan qattiq kristall tuzilmalar gaz gidrati deb ataladi.

Tabiiy gaz uzoq quvurlar orqali tashiladigan tizimlarda suyuqlikning atrof-muhitga issiqlik yo‘qotishi yoki bosim pasayishi sabab sovishi tizimni gidratlar barqaror bo‘ladigan harorat va bosim hududiga olib kelishi mumkin. Gidrat kristallarining o‘sishi va aglomeratsiyasi quvur kesimini toraytirishi, bosim yo‘qotishini oshirishi va keyingi bosqichda quvur tiqilib qolishiga olib kelishi mumkin.

Shu sabab faqat gidratlarning qaysi harorat va bosimda termodinamik barqaror ekanini bilish yetarli emas. Gidratning qanday tezlikda yadrolanishi, o‘sishi va qancha vaqtda operatsion xavf tug‘dirishi mumkinligini tushunish ham zarur.

Metan gidrati hosil bo‘lishi qaysi bosqichlardan o‘tadi?

Tadqiqotning 4-sahifadagi 1-rasmi tipik gaz sarfi egri chizig‘ini to‘rtta sifat bosqichiga ajratadi:

1. I-bosqich — Gazning erishi va induksiya: Gaz suv fazasida eriydi; barqaror gidrat yadrolari makroskopik ravishda sezilguncha induksiya davri kechadi.
2. II-bosqich — Tez o‘sish: Barqaror yadrolar o‘sadi va gaz tezlik bilan gidrat fazasiga olinadi. Bu gaz sarfi egri chizig‘ining eng tik hududidir.
3. III-bosqich — Sekinlashuvchi o‘sish: Kristall tuzilmalar o‘sishda davom etadi, biroq gaz sarfi tezligi pasaya boshlaydi.
4. IV-bosqich — Barqaror yoki cheklangan holat: Gaz sarfi juda sekinlashadi. Bu bosqich ko‘pincha massa uzatish cheklovlari bilan bog‘lanadi.

Tadqiqot ayniqsa gidrat o‘sish tajribalaridan hosil qilingan kinetik parametrga qaratilgan. Biroq haqiqiy tajribada interfeys reaksiyasi, gazning suyuqlik ichida tashilishi, kristall sirtiga diffuziya va hosil bo‘lish issiqligining muhitdan chiqarilishi bir vaqtning o‘zida yuz berishi mumkin.

“Kinetik parametr” nega ko‘rinma yoki samarali parametr?

Ideal intrinsik kinetik konstanta faqat gidrat-suv-gaz interfeysida molekulyar reaksiya tezligini ifodalashi kerak. O‘rganilgan tadqiqotdagi k parametri esa bevosita molekulyar o‘lchovdan olinmagan. U avval chop etilgan tajribalarda kuzatilgan interfeys harakat tezligi va issiqlik muvozanati ifodalaridan hisoblangan.

Shu sabab k qiymatida quyidagi ta’sirlar birgalikda bo‘lishi mumkin:

• Gidrat yuzasidagi haqiqiy reaksiya kinetikasi,
• Gazning suyuq fazada erishi va tashilishi,
• Interfeysdagi massa uzatish qarshiligi,
• Gidrat hosil bo‘lishidan ajralgan issiqlikning chiqarilishi,
• Tajriba idishi geometriyasi va aralashtirish sharoitlari,
• Bosim va fugasitelik hosil qiladigan termodinamik harakatlantiruvchi kuch,
• Gidrat plyonkasi yoki g‘ovak tuzilma ichida rivojlanadigan transport cheklovlari.

Tadqiqotchi shu sabab model “intrinsik reaksiya konstantasi”ni emas, balki eksperimental tizimning umumiy ta’sirlarini aks ettiruvchi ko‘rinma kinetik parametrni bashorat qilishini alohida ta’kidlaydi.

Arrhenius yondashuvining roli

An’anaviy kinetik tadqiqotlarda haroratga bog‘liqlik ko‘pincha Arrhenius tenglamasi bilan ifodalanadi:

\[ k=A\exp\left(-\frac{E_a}{RT_{\mathrm{eq}}}\right) \]

Bu yerda:

• k: Kinetik parametr yoki tezlik konstantasi,
• A: Arrhenius old eksponensial ko‘paytuvchisi,
• E_a: Aktivatsiya energiyasi,
• R: Universal gaz konstantasi,
• T_eq: Kelvin bo‘yicha muvozanat harorati.

Arrhenius modeli haroratning reaksiya tezligiga ta’siri uchun fizik asos beradi. Bunga qarshi standart shakl bosim, ortiqcha sovish, gaz fugasiteligi yoki eksperimental transport sharoitlarini alohida o‘zgaruvchilar sifatida ochiq kiritmaydi.

Tadqiqotning asosiy yondashuvi Arrhenius modelini bekor qilish emas; bosim va ortiqcha sovish kabi qo‘shimcha o‘zgaruvchilarni o‘z ichiga olgan shaffof regressiya modeli xuddi shu ma’lumotlar to‘plamidagi empirik moslikni yaxshilaydimi yoki yo‘qmi, shuni sinashdir.

Ishlatilgan ma’lumotlar to‘plami

Ma’lumotlar Freer va hamkasblarining 2001 yildagi metan gidrati plyonkasi o‘sishi tajribalaridan olingan. Jami 30 ta kuzatuv mavjud. Har bir kuzatuvda quyidagi o‘zgaruvchilar bor:

• Bosim, P,
• Gidrat muvozanat harorati, T_eq,
• Hajm yoki tajriba muhiti harorati, T_b,
• Ortiqcha sovish, ΔT = T_eq − T_b,
• Eksperimental ma’lumotlardan hisoblangan k parametri.

Faqat 30 ta kuzatuv mavjudligi model murakkabligini cheklashni talab qiladi. Shu sabab tadqiqotchi sun’iy neyron tarmog‘i kabi ko‘p parametrli usullar o‘rniga ko‘p o‘zgaruvchili chiziqli regressiyani asosiy model sifatida tanlagan.

Kinetik parametr qanday hisoblangan?

Tadqiqotda avval harakatlanuvchi gidrat interfeysi uchun issiqlik balansi berilgan:

\[ \lambda_H\rho_H\frac{dX}{dt}=K(T_{\mathrm{eq}}-T_b) \]

Bu yerda:

• λ_H: Gidrat hosil bo‘lish yoki dissotsiatsiya issiqligi,
• ρ_H: Metan gidrati zichligi,
• dX/dt: Gidrat interfeysining siljish tezligi,
• K: Umumiy reaksiya yoki transport koeffitsienti,
• T_eq − T_b: Ortiqcha sovish.

Umumiy koeffitsientning kinetik va issiqlik uzatish qarshiliklari bilan bog‘liqligi quyidagicha ifodalangan:

\[ \frac{1}{K}=\frac{1}{k}+\frac{1}{h} \]

Bu yerda h eksperimental tizimning issiqlik uzatish koeffitsientidir. Ikki ifoda qayta tartiblanganda:

\[ \frac{1}{k}=\frac{T_{\mathrm{eq}}-T_b}{\lambda_H\rho_H\left(dX/dt\right)}-\frac{1}{h} \]

tenglamasi olinadi.

Hisoblashda ishlatilgan konstantalar:

Interfeys tezligi tajriba manbasida µm/s birlikda berilgan va hisoblashdan oldin m/s birlikka aylantirilgan. Issiqlik tenglamasida watt bo‘yicha birlik izchilligini saqlash uchun gidrat issiqligi ham kJ/kg dan J/kg ga o‘tkazilishi kerak; tadqiqot hisoblash kodini ulashmagani sabab bu konversiyaning amaldagi tafsilotlarini bevosita tekshirish mumkin emas.

13-sahifadagi sochilish grafiklari nimani ko‘rsatadi?

2-rasmda k parametri to‘rtta o‘zgaruvchiga nisbatan alohida chizilgan:

• Bosim: Bosim oshgani sari k qiymatida kuchli va taxminan monoton o‘sish kuzatiladi.
• Muvozanat harorati: Muvozanat harorati oshgani sari k odatda ortadi.
• Hajm harorati: Ma’lumotlar 1, 2 va 4 °C atrofida to‘planadi va aniq chiziqli tendensiya ko‘rsatmaydi.
• Ortiqcha sovish: ΔT oshgani sari k qiymati aniq ortishi ko‘rinadi.

Tadqiqotchi muhim chet qiymat aniqlamaganini va kichik ma’lumotlar to‘plami sabab hech bir kuzatuvni olib tashlamaganini bildirgan. Biroq grafiklarda ko‘rinadigan munosabatlar faqat shu eksperimental oraliqqa tegishli va gidrat o‘sishi fundamental jihatdan chiziqli degani emas.

Korrelyatsiya tahlili

Bosim, muvozanat harorati va ortiqcha sovish k bilan 0,9 dan yuqori korrelyatsiya ko‘rsatadi. Bunga qarshi hajm haroratining bevosita bog‘liqligi juda zaif.

Biroq kirishlarning o‘zaro korrelyatsiyasi ham juda yuqori. Ayniqsa P bilan T_eq orasidagi 0,987 darajadagi korrelyatsiya jiddiy multikollinearlikka ishora qiladi. Shu sabab yakuniy model koeffitsientlarini mustaqil fizik ta’sirlar yoki sababiy o‘zgarishlar sifatida talqin qilish kerak emas.

Hajm harorati nega yakuniy modelga kiritilmagan?

Hajm harorati k bilan zaif bevosita korrelyatsiya ko‘rsatgan. Bundan tashqari ortiqcha sovish allaqachon quyidagi munosabat bilan hisoblanadi:

\[ \Delta T=T_{\mathrm{eq}}-T_b \]

T_b, T_eq va ΔT ni bir vaqtning o‘zida modelga kiritish matematik jihatdan takroriy ma’lumot hosil qiladi. Shu sabab tadqiqotchi hajm haroratini nomzod model to‘plamidan chiqargan va bosim, muvozanat harorati hamda ortiqcha sovishga e’tibor qaratgan.

Biroq ΔT ishlatilishi T_b ma’lumoti butunlay yo‘qolgan degani emas. Hajm haroratining ta’siri ortiqcha sovish o‘zgaruvchisi ichida bilvosita mavjud.

Random forest muhimlik tahlili

Tadqiqotda Pearson korrelyatsiyasiga qo‘shimcha ravishda eksplorativ random forest o‘zgaruvchi muhimlik qiymatlari berilgan:

Regressor natijalarida ortiqcha sovish eng yuqori, bosim va muvozanat harorati bir-biriga yaqin, hajm harorati esa juda past muhimlikka ega. Bu tartib korrelyatsiya tahlilidagi umumiy tendensiyaga mos.

Bunga qarshi uzluksiz maqsad o‘zgaruvchisi uchun “Random Forest Classifier” qanday qo‘llangani, k qiymatlari sinflarga aylantirilganmi yoki yo‘qmi va qaysi giperparametrlar ishlatilgani tushuntirilmagan. Tadqiqotchi ham bu tahlilni model tanlashning asosiy asosi sifatida ishlatmagan, faqat eksplorativ taqqoslash sifatida taqdim etgan.

Qaysi regressiya modellari taqqoslangan?

Ko‘p o‘zgaruvchili chiziqli regressiya tuzilishi

Umumiy ko‘p o‘zgaruvchili chiziqli regressiya modeli quyidagicha berilgan:

\[ Y=\beta_0+\beta_1X_1+\beta_2X_2+\cdots+\beta_nX_n+\varepsilon \]

Bu yerda Y bashorat qilinayotgan chiqish, X o‘zgaruvchilari kirishlar, β₀ doimiy had, β koeffitsientlari kirishlarning model og‘irliklari va ε tushuntirilmagan xato hadidir.

Ushbu tadqiqot uchun yakuniy Model 7 quyidagi tenglamaga ega:

\[ k=-37157{,}28+15892{,}08P-5901{,}97T_{\mathrm{eq}}+2162{,}22\Delta T \]

P MPa, T_eq va ΔT °C birlikda; bashorat qilinadigan k birligi W/m²K.

Modelda bosim koeffitsienti musbat, muvozanat harorati koeffitsienti manfiy va ortiqcha sovish koeffitsienti musbat. Muvozanat harorati yakka holda k bilan musbat korrelyatsiya ko‘rsatganiga qaramay ko‘p o‘zgaruvchili modelda uning koeffitsienti manfiy chiqishi o‘zgaruvchilar o‘rtasidagi juda kuchli bog‘lanish natijasidir. Bu koeffitsientlarni alohida fizik mexanizmlarning kattaligi sifatida talqin qilmaslik kerak.

16-sahifadagi usul oqimi

3-rasm tadqiqot jarayonini quyidagicha ko‘rsatadi:

1. Ma’lumotlarni tayyorlash va oldindan ishlov berish,
2. Barcha nomzod xususiyatlarni aniqlash,
3. Dastlabki eksplorativ tahlil,
4. Korrelyatsiya va fizik asos bilan o‘zgaruvchi tanlash,
5. Tanlangan kirishlar bilan MLR modellarini qurish,
6. Eksperimental kinetik parametrni bashorat qilish,
7. Modellarni umumiy validatsiya mezonlari bilan taqqoslash.

Oqim sxemasida chiqish “intrinsik tezlik parametri” sifatida belgilangan ko‘rinadi. Bunga qarshi asosiy matn hisoblangan parametr sof intrinsik konstanta emasligini va tizim darajasidagi ta’sirlarni o‘z ichiga olishini aniq aytadi. Rasmdagi ifoda bilan matndagi ehtiyotkorroq ta’rif o‘rtasida terminologik nomuvofiqlik bor.

10 qatlamli kross-validatsiya qanday qo‘llangan?

Jami ma’lumotlar to‘plami 10 qismga ajratilgan. Har bir aylanishda to‘qqiz qism o‘qitish, qolgan bir qism test uchun ishlatilgan; jarayon har bir qism bir marta test ma’lumoti bo‘lguncha takrorlangan.

Jami atigi 30 kuzatuv bo‘lgani uchun har bir test qatlamida taxminan uchta kuzatuv mavjud. Bu ma’lumotlardan samarali foydalanishga yordam bersa-da, har qatlamdagi R² va xato mezonlari oz sonli kuzatuvlardan ta’sirlanishi mumkin.

Tadqiqot qatlamlar qaysi random seed bilan yaratilganini, ma’lumotlar aralashtirilgan yoki yo‘qligini, kross-validatsiya necha marta takrorlanganini yoki qatlamlar o‘rtasidagi samaradorlik standart og‘ishini bildirmagan. Shu sabab boshqa qatlam bo‘linishida natijalar qanchalik o‘zgarishi noma’lum.

Bundan tashqari yetti nomzod model bir xil kross-validatsiya natijalari asosida taqqoslanib, eng yaxshisi tanlangan. Alohida tashqi test to‘plami yoki ichma-ich kross-validatsiya ishlatilmagani sabab bildirilgan yakuniy samaradorlik model tanlashdan kelib chiqadigan ma’lum darajadagi optimizmni o‘z ichiga olishi mumkin.

Ishlatilgan validatsiya mezonlari

Determinatsiya koeffitsienti

\[ R^2=\frac{\sum_{i=1}^{n}(\hat{y}_i-\bar{y})^2}{\sum_{i=1}^{n}(y_i-\bar{y})^2} \]

R² bashoratlar kuzatilgan o‘zgaruvchanlikni qanchalik ifodalashini ko‘rsatadi. Qiymatning 1 ga yaqinlashishi kuchliroq moslikni bildiradi. Biroq yuqori R² munosabat sababiy yoki fizik jihatdan to‘g‘ri ekanini o‘zi isbotlamaydi.

Ildiz o‘rtacha kvadratik xato

\[ RMSE=\sqrt{\frac{1}{N}\sum_{i=1}^{N}(\hat{y}_i-y_i)^2} \]

RMSE katta xatolarni kvadratlash sabab ko‘proq jazolaydi. Ushbu tadqiqotda uning birligi k bilan bir xil, ya’ni W/m²K.

O‘rtacha mutlaq xato

\[ MAE=\frac{1}{N}\sum_{i=1}^{N}|y_i-\hat{y}_i| \]

MAE o‘rtacha mutlaq bashorat farqini ko‘rsatadi va katta yakka xatolarga RMSE ga qaraganda kamroq sezgir.

Yetti regressiya modelining natijalari

Bir o‘zgaruvchili modellar orasida bosim modeli M1 eng yaxshi natijani bergan. Bosimning o‘zi R² = 0,9207 ga erishib, ma’lumotlar to‘plamidagi k o‘zgarishining katta qismini tushuntirgan.

Bosimga muvozanat harorati yoki ortiqcha sovish qo‘shilganda samaradorlik oshgan. Ikki o‘zgaruvchili modellar orasida P va ΔT ni ishlatgan M5 R² = 0,9480 bilan eng muvaffaqiyatli natijani bergan.

Uch o‘zgaruvchili M7 taqqoslangan barcha regressiyalar ichida eng yuqori R² hamda eng past RMSE va MAE qiymatlarini hosil qilgan.

Bashorat qilingan va eksperimental qiymatlarni taqqoslash

22-sahifadagi 4-rasm Model 7 bashorat qilgan k qiymatlarini eksperimental ravishda hosil qilingan qiymatlar bilan taqqoslaydi. Nuqtalarning katta qismi mukammal moslikni ifodalovchi diagonal chiziq atrofida to‘plangan.

Past va o‘rta k qiymatlarida moslik umuman kuchli ko‘rinadi. Bunga qarshi taxminan 50.000–65.000 W/m²K oralig‘ida ayrim nuqtalar ideal chiziqdan sezilarliroq uzoqlashgan. Bu yuqori k hududida mutlaq xatolar kattalashishi mumkinligini ko‘rsatadi.

Qoldiq grafigi nimani ko‘rsatadi?

23-sahifadagi Qo‘shimcha S1-rasm eksperimental qiymat bilan bashorat o‘rtasidagi farqni bashorat qilingan k qiymatiga nisbatan ko‘rsatadi. Qoldiqlar nol chizig‘ining ham yuqorisida, ham pastida joylashgan va tadqiqot aniq sistematik naqsh yo‘qligini bildiradi.

Biroq yuqori bashorat oralig‘ida taxminan +10.000 va −5.000 W/m²K darajasiga yaqin qoldiqlar ko‘rinadi. Ma’lumotlar soni kam bo‘lgani sabab qoldiqlarning normal taqsimoti, doimiy dispersiya farazi yoki ta’sirchan kuzatuvlar haqida kuchli xulosa chiqarib bo‘lmaydi.

Tadqiqotda qoldiqlar normalligi, geteroskedastiklik, leverage qiymatlari, Cook masofasi yoki bashorat intervallari bildirilmagan.

Arrhenius modeli bilan taqqoslash

Regressiya modeli Arrhenius modeliga nisbatan R² qiymatini 0,0142 ga oshirgan; RMSE ni 619,5 W/m²K va MAE ni 537,9 W/m²K ga kamaytirgan.

Arrhenius modeli taxminan 173 kJ/mol aktivatsiya energiyasi bergani bildirilgan. Biroq Arrhenius old eksponensial koeffitsienti, har bir kross-validatsiya qatlamida parametrlar qay tarzda qayta moslashtirilgani va aktivatsiya energiyasining noaniqligi tushuntirilmagan.

24-sahifadagi 5-rasmda Arrhenius hamda MLR bashoratlari ideal moslik chizig‘i atrofida joylashadi. Ikki yondashuv natijalarining bir-biriga ancha yaqinligi haroratning ma’lumotlar to‘plamidagi ustun rolini qo‘llab-quvvatlaydi. MLR ning qo‘shimcha afzalligi bosim va ortiqcha sovish ma’lumotini ham alohida ishlata olishidir.

MLR modelining yaxshiroq statistik mosligi nimani anglatadi?

Pastroq xato MLR modelining o‘rganilgan 30 kuzatuv ichidagi o‘zgarishni Arrhenius modeliga qaraganda biroz yaxshiroq kuzatishini ko‘rsatadi. Bu natija quyidagilarni anglatmaydi:

• MLR gidrat o‘sish fizikasini Arrhenius modelidan aniqroq ifodalashi,
• Bosim yoki harorat koeffitsientlari sababiy kattaliklar ekani,
• Model tadqiqot oralig‘idan tashqarida ishonch bilan ishlatilishi,
• Model turli reaktorlar va gaz tarkiblari uchun amal qilishi,
• Intrinsik molekulyar reaksiya konstantasi bashorat qilinishi.

Arrhenius tenglamasi harorat va aktivatsiya energiyasi orqali fizik talqin beradi, MLR esa ma’lumotlar to‘plamidagi ko‘p o‘zgaruvchili naqshni moslashuvchanroq ifodalaydi. Tadqiqot bu ikki yondashuvni raqib emas, bir-birini to‘ldiruvchi vositalar sifatida ko‘radi.

Bosimning kuchli ko‘rinishining ehtimoliy ma’nosi

Bosim bir o‘zgaruvchili modellar ichida eng kuchli bashorat qiluvchi bo‘lgan. Tadqiqotda bu holat bosimning gaz fugasiteligi, gidrat hosil bo‘lish harakatlantiruvchi kuchi va interfeysda mavjud bo‘lishi mumkin bo‘lgan metan miqdoriga ta’siri bilan bog‘langan.

Biroq ishlatilgan ma’lumotlar to‘plamida bosim bilan muvozanat harorati orasida 0,987 kabi nihoyatda kuchli korrelyatsiya bor. Shu sabab bosim koeffitsienti faqat bosimga xos ta’sirni ifodalaydi deb aytib bo‘lmaydi. Bosim o‘zgarishi tajriba dizaynida muvozanat haroratidagi o‘zgarish bilan birga kechgan.

Tadqiqotning kuchli tomonlari

• Kichik ma’lumotlar to‘plamiga mos, sodda va tushuntiriladigan model tanlangan.
• Yetti xil kirish kombinatsiyasi bir xil validatsiya usuli bilan taqqoslangan.
• Bosim, harorat va ortiqcha sovish ta’sirlari alohida va birgalikda baholangan.
• Model faqat o‘qitish mosligi bilan emas, 10 qatlamli kross-validatsiya bilan baholangan.
• R² bilan birga RMSE va MAE ham bildirilgan.
• Arrhenius modeli xuddi shu ma’lumot va validatsiya yondashuvi bilan taqqoslangan.
• Tadqiqotchi MLR fizik mexanizm emasligini aniq aytgan.
• Multikollinearlik va ma’lumotlar to‘plamiga xos amal qilish manbada tan olingan.
• Bashorat-haqiqiy va qoldiq grafiklari berilgan.
• Kelajakda fizikaga asoslangan va ma’lumotga asoslangan usullarni birlashtirish taklif qilingan.

Tadqiqotning asosiy cheklovlari

• Tadqiqot peer review’dan o‘tmagan preprint.
• Ma’lumotlar to‘plami faqat 30 kuzatuvdan iborat.
• Barcha ma’lumotlar bitta eski eksperimental tadqiqotdan olingan.
• Turli laboratoriya, reaktor yoki tajriba qurilmasiga tegishli mustaqil tashqi validatsiya yo‘q.
• Bosim bilan muvozanat harorati orasida juda yuqori korrelyatsiya mavjud.
• Regressiya koeffitsientlari barqarorligini ko‘rsatadigan standart xato yoki ishonch intervali berilmagan.
• Variance inflation factor kabi multikollinearlik mezonlari hisoblanmagan.
• Kross-validatsiya qatlamlarining random seed va takror soni tushuntirilmagan.
• Qatlamlar orasidagi samaradorlik taqsimoti yoki standart og‘ish bildirilmagan.
• Alohida tashqi test to‘plami yoki ichma-ich kross-validatsiya ishlatilmagan.
• Ma’lumotlarning bir xil tajriba seriyalari ichidagi bog‘liqlik tuzilishi baholanmagan.
• Model faqat chiziqli munosabatni faraz qiladi.
• O‘zaro ta’sir hadlari yoki chiziqli bo‘lmagan hadlar o‘rganilmagan.
• Qoldiqlarning normalligi va doimiy dispersiyasi statistik tekshirilmagan.
• Random Forest Classifier uzluksiz maqsad uchun qanday ishlatilgani tushuntirilmagan.
• Regressiya va tahlil kodlari ulashilmagan.
• Xom 30 kuzatuv maqola ichida to‘liq jadval sifatida berilmagan.
• Arrhenius old eksponensial koeffitsienti va aktivatsiya energiyasi noaniqligi berilmagan.
• Gaz fugasiteligi, interfeys maydoni, diffuzivlik va massa uzatish koeffitsienti bevosita modelga kiritilmagan.
• Model turli gaz tarkiblari, tuzlar, ingibitorlar yoki oqimli tizimlarda sinovdan o‘tkazilmagan.

Manbadagi muhim hisobot muammolari

• 16-sahifadagi usul sxemasi chiqishni “intrinsik tezlik parametri” deb belgilaydi, asosiy matn esa parametr ko‘rinma va tizim ta’sirlarini o‘z ichiga oluvchi kattalik ekanini bildiradi.
• Random Forest Classifier uzluksiz k maqsadi uchun qanday qo‘llangani tushuntirilmagan.
• R² izohida “R-squared / adjusted R-squared” ifodasi ishlatilgan bo‘lsa-da, jadvalda adjusted R² alohida bildirilmagan.
• Kross-validatsiyaning tasodifiy bo‘linish tafsilotlari va takror soni berilmagan.
• Arrhenius moslashuvining old eksponensial ko‘paytuvchisi va noaniqliklari taqdim etilmagan.
• Hisoblarni qayta ishlab chiqish imkonini beradigan kod va to‘liq ma’lumotlar jadvali ulashilmagan.
• Model 7 ning juda kuchli bog‘langan kirishlariga qaramay koeffitsient noaniqliklari berilmagan.
• K parametrini hisoblashda ishlatilgan birlik konversiyalari kod darajasida tushuntirilmagan.

Tadqiqot qaysi natijalarni qo‘llab-quvvatlaydi?

• O‘rganilgan 30 kuzatuv ichida bosim, muvozanat harorati va ortiqcha sovish k parametri bilan kuchli statistik bog‘lanish ko‘rsatgan.
• Hajm harorati yakka holda k bilan zaif chiziqli bog‘lanish ko‘rsatgan.
• Bosim bir o‘zgaruvchili modellar ichida eng muvaffaqiyatli bashorat qiluvchi bo‘lgan.
• P, T_eq va ΔT ni birgalikda ishlatgan M7 sinovdan o‘tkazilgan regressiyalar ichida eng past kross-validatsiya xatosini bergan.
• Arrhenius modeli ham xuddi shu ma’lumotlar to‘plamida kuchli bashorat samaradorligini ko‘rsatgan.
• Qo‘shimcha operatsion o‘zgaruvchilardan foydalanish ma’lumotlar to‘plamiga xos empirik moslikni yaxshilagan.
• Sodda chiziqli modellar kichik gidrat kinetikasi ma’lumotlar to‘plamlarini o‘rganish uchun shaffof surrogat vositalar bo‘lishi mumkin.
• MLR va Arrhenius yondashuvlari turli maqsadlarda bir-birini to‘ldirishi mumkin.

Tadqiqot qaysi natijalarni isbotlamaydi?

• Model molekulyar miqyosdagi intrinsik gidrat o‘sish qonuni ekanini isbotlamaydi.
• Regressiya koeffitsientlari mustaqil sababiy ta’sirlar ekanini ko‘rsatmaydi.
• Model boshqa reaktorlar, quvurlar yoki gaz tarkiblarida amal qilishini ko‘rsatmaydi.
• Haqiqiy quvur liniyasining tiqilish vaqtini bashorat qilmaydi.
• Gidratga qarshi yoki dispersant kimyoviy moddaga ehtiyojni hisoblamaydi.
• Issiqlik va massa uzatish ta’sirlarini bir-biridan ajratmaydi.
• MLR fizik jihatdan Arrhenius modelidan ustun ekanini isbotlamaydi.
• Yuqori R² qiymati yangi va umumiy fizik qonun hosil qilganini ko‘rsatmaydi.
• Model ma’lumotlar oralig‘idan tashqaridagi bosim va haroratlarda xavfsiz ekstrapolyatsiya qila olishini ko‘rsatmaydi.

Tadqiqot usuli va topilmalari

Usul xulosasi

Yakuniy model va qo‘llash chegarasi

\[ k=-37157{,}28+15892{,}08P-5901{,}97T_{\mathrm{eq}}+2162{,}22\Delta T \]

Ushbu tenglama faqat tadqiqotdagi taxminiy oraliqlarda baholangan:

Tenglamaning ushbu oraliqlardan tashqarida ishlatilishi tadqiqot tomonidan tasdiqlanmagan. Chiziqli model fizik chegaralarni majburan qo‘llamagani uchun mos kelmaydigan kirishlarda manfiy yoki real bo‘lmagan bashoratlar hosil qilishi mumkin.

Asosiy samaradorlik taqqoslanishi

Model 7 ning o‘rtacha mutlaq xatosi ma’lumotlar to‘plamidagi o‘rtacha k qiymati bo‘lgan taxminan 29.696 W/m²K ning taxminan sakkiz foiziga to‘g‘ri keladi. Biroq bu nisbat kuzatuvlar orasidagi xato taqsimotini yoki yuqori k hududidagi kattaroq xatolarni yolg‘iz o‘zi ko‘rsatmaydi.

Tadqiqotchining kelajak ish takliflari

• Eksperimental ma’lumotlar bazasini kengaytirish,
• Kengroq bosim va harorat oraliqlarini o‘rganish,
• Turli gaz tarkiblari va reaktor konfiguratsiyalarini qo‘shish,
• Gaz fugasiteligini bevosita kirish sifatida ishlatish,
• Interfeys maydoni, diffuzivlik va massa uzatish koeffitsientini modelga qo‘shish,
• Yana ilg‘or mashinaviy o‘rganish modellarini ehtiyotkorlik bilan taqqoslash,
• Kichik ma’lumotlar to‘plamida izohlanuvchanlikni saqlash,
• Fizikaga asoslangan gidrat modellari bilan ma’lumotga asoslangan usullarni gibrid tarzda birlashtirish.

Statistik baholash chegarasi

Tadqiqotdagi samaradorlik qiymatlari bitta 30 kuzatuvli ma’lumotlar to‘plamida olingan. Mustaqil eksperimental tashqi test ma’lumoti yo‘qligi sabab R² = 0,9682 qiymatini umumiy foydalanish muvaffaqiyati sifatida talqin qilib bo‘lmaydi.

Bosim, muvozanat harorati va ortiqcha sovish o‘rtasidagi kuchli korrelyatsiya model bashoratini albatta yaroqsiz qilmaydi; biroq koeffitsientlarning o‘zgaruvchilar orasida sezgir taqsimlanishiga va yangi ma’lumotda beqarorlashishiga olib kelishi mumkin. Shu sabab modelning asosiy foydalanish sohasi mustaqil fizik talqin emas, balki o‘xshash sharoitdagi k qiymatlarini taxminiy bashorat qilishdir.

Manba va usul eslatmasi

Tadqiqotning to‘liq asl nomi: A Data-Driven Approach for Estimating an Experimentally Derived Kinetic Parameter in Methane Hydrate Growth

Muallif: Chinedu Charles Alamezie

Mualliflar tartibi: Tadqiqotda bitta muallif bor.

Teng birinchi muallif yoki teng hissa: Mavjud emas.

Mas’ul muallif: Chinedu Charles Alamezie

Muassasa aloqasi: School of Engineering, University of Aberdeen, Aberdeen AB24 3UE, Scotland, United Kingdom

Mas’ul muallif elektron pochtasi: ccalamezie@hotmail.com

ORCID: Yuklangan versiyada ORCID ma’lumoti mavjud emas.

DOI:10.2139/ssrn.6944011

Jurnal: Peer-reviewed jurnalda nashr etilmagan.

Asl jurnal nashriyoti: Mavjud emas.

Nashr platformasi: SSRN

Nashr sanasi: 15 June 2026

Sahifalar soni: 33

Manba turi: Ikkilamchi eksperimental ma’lumot tahlili va regressiya modellashtirishni o‘z ichiga olgan preprint tadqiqot

Peer-review holati: Tadqiqot peer review’dan o‘tmagan. Yuklangan versiyaning har bir sahifasida preprint peer-reviewed emasligi ko‘rsatilgan.

Rasmiy havola:SSRN rasmiy preprint sahifasi

Ushbu o‘zbekcha maqola yuklangan tadqiqotning matni, tenglamalari, to‘rtta asosiy rasmi, qo‘shimcha qoldiq grafigi, oltita jadvali, korrelyatsiya matritsasi, random forest muhimlik qiymatlari, yetti regressiya tenglamasi, kross-validatsiya natijalari, Arrhenius taqqoslanishi, xulosalari va tavsiyalari ko‘rib chiqilib tayyorlangan.

Tashqi manbalar faqat tadqiqot muallifi nomi, muassasasi, SSRN nashr sanasi, DOI, sahifa soni va rasmiy havolasini bibliografik tasdiqlash uchun ishlatilgan. Tadqiqotning ilmiy topilmalariga tashqi manbalardan yangi natija qo‘shilmagan.

Tadqiqotning asosiy cheklovlari: kichik va bitta manbali ma’lumotlar to‘plami, kuchli o‘zgaruvchi bog‘lanishlari, mustaqil tashqi validatsiya yo‘qligi, takroriy yoki ichma-ich kross-validatsiya qo‘llanmagani, regressiya koeffitsientlari noaniqliklari berilmagani va bashorat qilinayotgan parametr sof intrinsik kinetikadan ko‘ra termodinamik va transport ta’sirlarini birgalikda o‘z ichiga olishi. Model fizik gidrat o‘sish tenglamalarini almashtiradigan umumiy qonun sifatida emas, o‘xshash eksperimental sharoitlarda ishlatilishi mumkin bo‘lgan shaffof va ma’lumotlar to‘plamiga xos surrogat model sifatida baholanishi kerak.