Wimbi la msukumo linalozalishwa na maporomoko ya ardhi ni wimbi la uso linalotokea wakati nyenzo ya maporomoko inayosonga kwa kasi inaingia katika maji na kuhamisha masa na momentamu kwenda kwenye maji. Utafiti huu unachunguza jinsi hali za awali za maporomoko ya viskoplastiki zinavyobadilika kuwa amplitudo ya wimbi na urefu wa wimbi, si kama uhusiano wa moja kwa moja tu kati ya mwanzo na matokeo, bali kama mnyororo wa mchakato wenye hatua tatu unaowakilisha wazi hali halisi ya kuingia majini iliyo kati yao. Jumla ya majaribio 290 ya maabara yalitumika; hali za awali za maporomoko, mienendo halisi ya kuingia majini na majibu ya mawimbi zilichanganuliwa kama tabaka tofauti za vigezo. Hoja kuu ni kwamba majibu ya mawimbi yanahusiana kwa karibu zaidi na masa ya nyenzo inayotokea kuingia majini na inersia ya kuingia kuliko na hali za awali. Hata hivyo, matokeo yanahusu nyenzo mahususi ya Carbopol, jiometri ya mfereji na masafa ya vigezo vya majaribio; njia za kitakwimu si uthibitisho wa usababisho, na uthibitishaji huru wa uwanjani unahitajika kabla ya kuhamisha matokeo moja kwa moja kwa maporomoko ya asili.
Katika utafiti, hali za awali zilipangwa kama \(X_0=(M_i,L,Bi)\), hali halisi ya kuingia majini kama \(X_e=(Fr,S,M)\), na jibu la wimbi kama \(Y=(A,H)\). Hapa \(M_i\) inaonyesha masa ya awali ya uwiano, \(L\) umbali wa kuteleza wa uwiano, \(Bi\) nambari ya Bingham, \(Fr\) nambari ya Froude ya kuingia, \(S\) unene wa kuingia wa uwiano, \(M\) masa tendaji ya kuingia ya uwiano, \(A\) amplitudo ya juu ya wimbi ya uwiano na \(H\) urefu wa juu wa wimbi wa uwiano. Hivyo, muundo mkuu wa kisayansi ni mnyororo wa mchakato \(X_0 \rightarrow X_e \rightarrow Y\).
Masa tendaji ya kuingia ya uwiano \(M\) ilionyesha uwiano wa Spearman wa 0,96 na \(A\) na 0,91 na \(H\). Kwa nambari ya Froude ya kuingia \(Fr\), uwiano husika ulikuwa 0,71 na 0,80. Katika uchanganuzi wa upatanishi, sehemu isiyo ya moja kwa moja ya masa ya awali ya uwiano \(M_i\) kupitia \(M\) ilikuwa 0,652 kwa \(A\) na 0,600 kwa \(H\); sehemu isiyo ya moja kwa moja ya nambari ya Bingham kupitia \(M\) ilikuwa mtawalia −0,669 na −0,615. Thamani hizi zinaonyesha jukumu kuu la njia ya kuingia majini inayohusiana na masa ndani ya seti ya data iliyochunguzwa na modeli iliyojengwa.
Mbinu ya mfululizo wa hatua inayotegemea Gradient Boosting ilianza kwa kujenga upya \(Fr\), \(S\) na \(M\) kutoka hali za awali, kisha ikatabiri jibu la wimbi. Hata hivyo, modeli ya mfululizo haikutoa ongezeko kubwa la usahihi ikilinganishwa na modeli ya moja kwa moja inayotumia hali za awali pekee. Kwa hiyo faida kuu ya mfumo wa mfululizo si usahihi wa utabiri; ni kufanya njia za kati za kifizikia na kitakwimu kati ya hali za awali, kuingia majini na uundaji wa wimbi zionekane wazi, na kuruhusu vyanzo vya makosa kutenganishwa kwa hatua.
Tatizo kuu la kisayansi
Wimbi linalozalishwa maporomoko yanapoingia majini halitegemei tu masa ya awali ya maporomoko au nafasi yake ya awali. Maporomoko yanaposonga kwenye mteremko yanaweza kuongezeka kasi, kutawanyika, kuwa membamba au mazito, na kupitia deformesheni ya ndani. Kwa hiyo masa, unene na hali ya reolojia iliyowekwa mwanzoni inaweza kuwa tofauti na hali halisi ya mwendo inapofika kwenye uso wa maji.
Tofauti hii ni muhimu hasa kwa nyenzo za viskoplastiki zinazoweza kubadilika umbo. Ingawa kizuizi kigumu kinaweza kuingia majini karibu kama kitu kimoja cha kijiometri, masa ya viskoplastiki hubadilisha umbo wakati wa kusonga kwenye mteremko. Kwa sababu hiyo, huenda si masa yote ya awali yanayoshiriki kwa kiwango sawa katika kuhamisha maji wakati muhimu wa uundaji wa wimbi.
Sehemu kubwa ya mbinu za awali zinazoendeshwa na data ziliunganisha vigezo vya awali moja kwa moja na sifa za mawimbi au zilitengeneza viwango vya umuhimu wa vigezo. Mbinu hii inaweza kuwa muhimu kwa utabiri; lakini ina kikomo katika kutenganisha ikiwa hali ya awali inaathiri wimbi moja kwa moja au ikiwa athari yake hutokea kwanza kwa kubadilisha masa ya kuingia majini, kasi ya kuingia au jiometri ya kuingia.
Mchango wa utafiti huu si kutengeneza algoriti mpya ya takwimu au kujifunza kwa mashine. Ubunifu wake ni kuunganisha mbinu zilizopo za uwiano, uchanganuzi wa njia, uchanganuzi wa upatanishi na Gradient Boosting katika mnyororo mmoja wa mchakato uliopangwa kulingana na mfuatano wa matukio ya kifizikia.
Kwa nini hali ya kuingia majini inachukuliwa kama hatua tofauti katika mawimbi ya maporomoko?
Kwa sababu nyenzo ya viskoplastiki ya maporomoko inaweza kubadilika umbo, kutawanyika, kuwa membamba au mazito na kupata kasi tofauti kati ya hali ya awali ilipoachiliwa na wakati inapofika majini; kwa hiyo masa, jiometri na hali ya inersia zinazobainisha wimbi moja kwa moja si lazima ziwe sawa na hali zilizowekwa mwanzoni.
Mfumo wa kifizikia wa majaribio
Mfumo wa maabara ni mfereji wa pande mbili wa maporomoko–wimbi. Vipengele vikuu ni tangi linaloshikilia nyenzo ya maporomoko, lango la nyumatiki linalofunguka, njia ya PVC yenye mteremko unaoweza kurekebishwa, mfereji mwembamba wa uwazi wa mawimbi na mfumo wa upigaji picha wa kasi kubwa.
| Kipengele | Sifa iliyotolewa katika chanzo |
|---|---|
| Urefu wa mfereji wa wimbi | 2,5 m |
| Upana wa mfereji wa wimbi \(B\) | 0,12 m |
| Mteremko unaoweza kurekebishwa wa njia | 30–50° |
| Mteremko uliotumika katika majaribio \(\theta\) | 40–45° |
| Kina kisichobadilika cha maji \(h_0\) | 0,20 m |
| Kamera | 400 fps |
| Azimio la picha | 600 × 800 pikseli |
| Jumla ya visa vya majaribio | 290 |
Lango la nyumatiki linapofunguliwa, nyenzo ya Carbopol husonga chini kwenye njia kwa mvuto na kuingia majini. Kuingia majini huhamisha maji yanayozunguka na kupeleka momentamu kwenye maji, hivyo kuunda usumbufu wa awali wa uso huru; usumbufu huu baadaye hubadilika kuwa wimbi la msukumo.
Eneo la upigaji picha lilijumuisha eneo la kuingia na uundaji wa wimbi katika uga wa karibu. Kabla ya majaribio, gridi ya kalibresheni yenye vipimo vinavyojulikana ilitumika kubadilisha kuratibu za pikseli kuwa umbali wa kifizikia. Nafasi ya kamera na uga wa kuona havikubadilishwa katika majaribio yote.
Amplitudo ya wimbi \(a_m\) ilifafanuliwa kama uhamisho wa juu kabisa wa uso huru kwenda juu ikilinganishwa na kiwango cha maji tulivu. Urefu wa juu wa wimbi \(h_m\) ni umbali wa juu kabisa wa wima kati ya kilele na bonde la wimbi linalolingana. Ingawa viwango hivi viwili vinaweza kuonekana kufanana, si kipimo kilekile cha kifizikia.
Masa tendaji ya kuingia ni nini?
Masa tendaji ya kuingia \(m_e\) ni masa ya nyenzo ya maporomoko iliyokuwa imeingia katika maji wakati urefu wa juu wa wimbi \(h_m\) ulipotokea. Hesabu inategemea kontua ya maporomoko iliyotolewa kutoka kwenye picha na msongamano unaojulikana wa nyenzo.
Ufafanuzi huu ni tofauti na masa ya awali. Masa ya awali \(m_i\) inaashiria jumla ya nyenzo iliyokuwapo kabla ya jaribio kuanza, ilhali \(m_e\) inawakilisha nyenzo iliyoshiriki kweli katika kuhamisha maji wakati muhimu wa uundaji wa wimbi kuu.
Kwa nini masa tendaji ya kuingia ni muhimu sana katika jibu la wimbi?
Kwa sababu masa tendaji ya kuingia ya uwiano \(M\) hubeba matokeo ya pamoja ya deformesheni, kutawanyika, upinzani wa reolojia na ugawaji upya wa maporomoko kabla hayajafika majini, inawakilisha moja kwa moja zaidi kuliko masa ya awali kiasi cha nyenzo kilichohamisha maji katika hatua kuu ya uundaji wa wimbi; katika utafiti huu \(M\) ilionyesha uhusiano wenye nguvu zaidi wa hali ya kuingia na amplitudo ya wimbi ya uwiano \(A\) pamoja na urefu wa wimbi wa uwiano \(H\).
Kwa nini Carbopol ilitumika?
Carbopol ni nyenzo ya viskoplastiki iliyorahisishwa yenye mkazo wa kuyield ulio wazi na tabia ya kupungua mnato inapokatwa. Tabia yake ya reolojia baada ya kuyield ilielezwa kwa modeli ya Herschel–Bulkley:
\[ \tau=\tau_c+K\dot{\gamma}^{\,n}, \qquad \tau>\tau_c \]
Hapa \(\tau\) ni mkazo wa kukata, \(\tau_c\) mkazo wa kuyield, \(K\) mgawo wa uthabiti, \(\dot{\gamma}\) kiwango cha kukata na \(n\) kiashiria cha tabia ya mtiririko. Ikiwa mkazo wa kukata unaotumika hauzidi \(\tau_c\), nyenzo haiingii katika mtiririko endelevu. Kwa hiyo \(\tau_c\) inawakilisha upinzani wa kizingiti unaopaswa kuzidiwa kabla ya deformesheni na mtiririko endelevu kuanza.
Utafiti ulitumia \(\tau_c\) kama kigezo kikuu kinachowakilisha tabia ya reolojia ya nyenzo ya maporomoko na ukachunguza uhusiano wake kutoka hali ya awali hadi inersia ya kuingia majini na masa tendaji ya kuingia kupitia nambari ya Bingham.
Masafa ya vigezo vya majaribio
| Kigezo | Maana | Masafa / thamani |
|---|---|---|
| \(\tau_c\) | Mkazo wa kuyield wa nyenzo | 60–90 Pa |
| \(m_i\) | Masa ya awali ya maporomoko | 2,0–5,6 kg |
| \(l_s\) | Umbali wa kuteleza | 0,55–0,75 m |
| \(\theta\) | Pembe ya mteremko | 40–45° |
| \(s_0\) | Unene wa kuingia kwenye uso wa maji | 0,026–0,050 m |
| \(u_0\) | Kasi ya kuingia | 1,12–2,36 m/s |
| \(h_0\) | Kina cha maji tulivu | 0,20 m |
Kutokuwa na uhakika wa vipimo kulitathminiwa kwa kutumia azimio la nafasi lililotokana na gridi ya kalibresheni na azimio la muda la kamera. Aidha, hali moja wakilishi ya jaribio ilirudiwa mara tatu kwa kujitegemea chini ya mipangilio ileile, na tofauti katika vipimo vilivyorudiwa ilitumika kutathmini kutokuwa na uhakika wa urudiwaji kati ya majaribio.
Lengo la mfumo wa vigezo visivyo na vipimo
Ili hali za majaribio zenye ukubwa tofauti wa kifizikia ziweze kulinganishwa katika mfumo mmoja wa kihisabati, vigezo vya kifizikia vilifanywa visivyo na vipimo.
Hali za awali:
\[ M_i=\frac{m_i}{\rho_w B h_0^2}, \qquad L=\frac{l_s}{h_0}, \qquad Bi=\frac{\tau_c}{\rho_s g s_i\sin\theta} \]
\(M_i\) ni masa ya awali iliyopimwa kwa msongamano wa maji, upana wa mfereji na kina cha maji. \(L\) hupima umbali wa kuteleza kwa kina cha maji. \(Bi\) inaonyesha uwiano wa upinzani wa kuyield wa nyenzo dhidi ya mkazo wa mvuto unaoendesha mwendo kwenye mteremko katika usanidi wa awali.
Hali halisi ya kuingia majini:
\[ Fr=\frac{u_0}{\sqrt{g h_0}}, \qquad S=\frac{s_0}{h_0}, \qquad M=\frac{m_e}{\rho_w B h_0^2} \]
\(Fr\) ni kipimo kisicho na vipimo cha kasi ya kuingia majini ikilinganishwa na kasi bainifu ya wimbi inayobainishwa na mvuto na kina cha maji. \(S\) hulinganisha unene wa kuingia na kina cha maji. \(M\) ni umbo lisilo na vipimo la masa inayoshiriki kwa ufanisi katika hatua kuu ya uundaji wa wimbi.
Majibu ya wimbi:
\[ A=\frac{a_m}{h_0}, \qquad H=\frac{h_m}{h_0} \]
\(A\) ni uwiano wa mkengeuko wa juu kabisa wa uso huru kwenda juu; \(H\) ni uwiano wa urefu wa wimbi kati ya kilele na bonde dhidi ya kina cha maji tulivu.
Mnyororo wa mchakato wa hatua tatu
| Hatua ya 1 | Hatua ya 2 | Hatua ya 3 | ||
|---|---|---|---|---|
| Hali za awali \(X_0=(M_i,L,Bi)\) | → | Hali halisi ya kuingia majini \(X_e=(Fr,S,M)\) | → | Jibu la wimbi \(Y=(A,H)\) |
Muundo huu unawezesha kuendelea kuchunguza mahusiano ya moja kwa moja kutoka hali za awali hadi jibu la wimbi huku pia ukitathmini kando ni sehemu gani ya uhusiano huo inaelezwa kupitia \(Fr\), \(S\) na \(M\).
Uwekaji sanifu na modeli ya njia
Kabla ya uchanganuzi wa njia, vigezo vyote viliwekwa sanifu kwa z-skor:
\[ Z_i=\frac{X_i-\mu_i}{\sigma_i} \]
Hapa \(X_i\) ni kigezo asilia, \(\mu_i\) wastani, \(\sigma_i\) mkengeuko sanifu na \(Z_i\) thamani iliyowekwa sanifu. Hatua hii huruhusu kulinganisha vigawo vya vigezo vilivyo katika mizani tofauti.
Njia ya hatua ya kwanza ilijengwa kutoka kwa vigezo vya awali kwenda kwa kila kigezo cha kuingia majini:
\[ Z_m=\alpha_{0m}+\sum_{k=1}^{3}a_{km}Z_k+\varepsilon_m \]
Hapa \(Z_k\) inawakilisha \(M_i\), \(L\) au \(Bi\); \(Z_m\) inawakilisha \(Fr\), \(S\) au \(M\); na \(a_{km}\) ni mgawo wa njia uliowekwa sanifu unaohusika.
Modeli ya hatua ya pili huunganisha vigezo vya awali na vigezo vya kuingia majini kwa pamoja na jibu la wimbi:
\[ Z_y=\beta_{0y}+\sum_{k=1}^{3}c'_{ky}Z_k+\sum_{m=1}^{3}b_{my}Z_m+\varepsilon_y \]
Hapa \(c'_{ky}\) inawakilisha uhusiano wa moja kwa moja wenye masharti uliobaki kutoka kwa kigezo cha awali hadi jibu la wimbi baada ya kuzingatia vigezo vya kuingia; \(b_{my}\) inawakilisha njia kutoka kwa kigezo cha kuingia hadi jibu la wimbi.
Sehemu isiyo ya moja kwa moja kupitia mpatanishi maalum:
\[ IE_{kmy}=a_{km}b_{my} \]
Jumla ya sehemu zisizo za moja kwa moja:
\[ TIE_{ky}=\sum_{m=1}^{3}IE_{kmy} \]
Na sehemu ya jumla:
\[ TE_{ky}=c'_{ky}+TIE_{ky} \]
Katika istilahi ya utafiti wa chanzo hizi huitwa “athari” za moja kwa moja, zisizo za moja kwa moja na za jumla; hata hivyo, ingawa zinaungwa mkono na mfuatano wa muda wa majaribio, ni mahusiano ya kitakwimu yenye masharti ndani ya modeli ya mstari iliyojengwa. Kwa zenyewe hazithibitishi usababisho.
Mantiki ya utabiri wa mfululizo
Ili kupima ikiwa hali halisi ya kuingia majini inaweza kutabiriwa kutokana na hali za awali pekee:
\[ \widehat{X}_e=f_1(X_0) \]
na kwa uwazi:
\[ (\widehat{Fr},\widehat{S},\widehat{M})=f_1(M_i,L,Bi) \]
modeli ya kutoka-hali-za-awali-hadi-kuingia ilijengwa.
Kisha mifumo mitatu ya utabiri wa mawimbi ililinganishwa:
| Modeli | Ingizo | Matumizi ya kisayansi |
|---|---|---|
| Modeli ya moja kwa moja ya hali za awali | \(X_0\) | Utabiri wa mbele kutoka taarifa zinazojulikana mwanzoni |
| Modeli ya hali ya kuingia iliyozingatiwa | \(X_e\) iliyopimwa | Ulinganisho wa maelezo / uchunguzi |
| Modeli ya mfululizo | \(X_0\rightarrow\widehat{X}_e\rightarrow\widehat{Y}\) | Utabiri wa mbele kutoka hali za awali + tafsiri ya mchakato wa kati |
Muhtasari wa kihisabati wa modeli ya mfululizo ni:
\[ \widehat{Y}=f_2[f_1(X_0)] \]
.
Mbinu na Matokeo ya Utafiti
Seti ya data na migawanyo ya msingi
| Kigezo | Wastani | Mkengeuko sanifu | Kiwango cha chini | Kiwango cha juu |
|---|---|---|---|---|
| \(M_i\) | 0,739 | 0,141 | 0,357 | 1,120 |
| \(L\) | 3,405 | 0,515 | 2,325 | 4,760 |
| \(Bi\) | 4,003 | 0,853 | 2,475 | 5,536 |
| \(Fr\) | 1,136 | 0,126 | 0,723 | 1,479 |
| \(S\) | 0,221 | 0,028 | 0,149 | 0,289 |
| \(M\) | 0,173 | 0,056 | 0,028 | 0,358 |
| \(A\) | 0,355 | 0,098 | 0,044 | 0,585 |
| \(H\) | 0,498 | 0,148 | 0,044 | 0,833 |
Uwiano kati ya hatua
Masa ya awali ya uwiano \(M_i\) ilionyesha uwiano chanya wa Spearman wa 0,61 na unene wa kuingia \(S\), na 0,36 na masa tendaji ya kuingia \(M\). Nambari ya Bingham \(Bi\) ilionyesha uwiano wa −0,79 na nambari ya Froude ya kuingia \(Fr\), na −0,29 na \(M\). Umbali wa kuteleza wa uwiano \(L\) ulionyesha uwiano chanya wa 0,33 na \(Fr\), na uwiano hasi wa −0,60 na −0,26 mtawalia na \(S\) na \(M\).
| Kigezo cha kuingia | Uwiano na \(A\) | Uwiano na \(H\) |
|---|---|---|
| \(Fr\) | 0,71 | 0,80 |
| \(S\) | 0,55 | 0,42 |
| \(M\) | 0,96 | 0,91 |
Katika ulinganisho huu, \(M\) inaonyesha uhusiano wenye nguvu zaidi na vigezo vyote viwili vya mawimbi. \(Fr\) ni kigezo cha pili kwa nguvu. Uwiano wa \(S\) ni chanya lakini ni mdogo zaidi.
Hali za awali zinapolinganishwa moja kwa moja na majibu ya mawimbi, mahusiano kwa ujumla ni dhaifu zaidi. Katika kundi hili, kigezo cha awali chenye nguvu zaidi ni \(Bi\), kikiwa na uwiano wa −0,44 kwa \(A\) na −0,56 kwa \(H\).
Nambari ya Bingham inaathirije uundaji wa wimbi?
Katika seti hii ya data, \(Bi\) ya juu inawakilisha upinzani mkubwa zaidi wa kuyield wa nyenzo ukilinganishwa na mkazo unaoendesha mwendo kwenye mteremko, na inaonyesha njia hasi zenye nguvu kwa nambari ya Froude ya kuingia \(Fr\) na masa tendaji ya kuingia \(M\); kwa hiyo upinzani mkubwa zaidi wa reolojia unahusishwa, ndani ya modeli, na inersia ndogo ya kuingia na nyenzo chache zinazoshiriki kwa ufanisi katika uzalishaji wa wimbi.
Vigawo vya njia vilivyowekwa sanifu
| Kigezo cha awali | \(Fr\) | \(S\) | \(M\) |
|---|---|---|---|
| \(M_i\) | 0,521 | 0,839 | 1,006 |
| \(L\) | 0,033; si muhimu | −0,741 | −0,622 |
| \(Bi\) | −1,061 | −0,261 | −1,031 |
Masa ya awali inapoongezeka, vigezo vyote vitatu vya kuingia majini vina mwelekeo wa kuongezeka. Kinyume chake, \(Bi\) ya juu inaonyesha mahusiano yenye masharti yanayopunguza kwa nguvu hasa \(Fr\) na \(M\). Njia huru ya umbali wa kuteleza kwenda \(Fr\) ni ndogo na si muhimu kitakwimu, ilhali njia zake kwenda \(S\) na \(M\) ni hasi kwa wazi.
| Mpatanishi | \(A\) | \(H\) |
|---|---|---|
| \(M\) | 0,649 | 0,596 |
| \(Fr\) | 0,316 | 0,458 |
| \(S\) | 0,072 | 0,100 |
\(M\) na \(Fr\) hubeba vigawo chanya vinavyotawala kuelekea majibu yote mawili ya mawimbi. Vigawo vya hatua ya pili vya \(S\) ni vidogo kiasi na havikufikia umuhimu wa kitakwimu katika uchanganuzi wa chanzo. Njia za moja kwa moja za vigezo vya awali baada ya kuzingatia hali ya kuingia pia kwa kiasi kikubwa ni ndogo; kati yao, ni njia ya \(L\rightarrow A\) pekee yenye −0,114 ndiyo muhimu kitakwimu.
Njia za upatanishi
| Kigezo cha awali | Njia | \(A\) | \(H\) |
|---|---|---|---|
| \(M_i\) | kupitia \(Fr\) | 0,165 | 0,239 |
| \(M_i\) | kupitia \(S\) | 0,060 | 0,084 |
| \(M_i\) | kupitia \(M\) | 0,652 | 0,600 |
| \(L\) | kupitia \(Fr\) | 0,010 | 0,015 |
| \(L\) | kupitia \(S\) | −0,053 | −0,074 |
| \(L\) | kupitia \(M\) | −0,403 | −0,371 |
| \(Bi\) | kupitia \(Fr\) | −0,335 | −0,486 |
| \(Bi\) | kupitia \(S\) | −0,019 | −0,026 |
| \(Bi\) | kupitia \(M\) | −0,669 | −0,615 |
Sehemu ya jumla iliyowekwa sanifu ya \(M_i\) juu ya \(A\) ni 0,945, na sehemu ya jumla juu ya \(H\) ni 0,855. Kwa \(Bi\), hizi ni −1,098 na −1,117; kwa \(L\), −0,560 na −0,471. Kwa vigezo vyote vitatu vya awali, sehemu nyingi kubwa zaidi zisizo za moja kwa moja zilitokea kupitia \(M\).
Gradient Boosting na uthibitishaji mtambuka
Gradient Boosting Regression ilitumika kujenga upya vigezo vya kuingia majini kutoka hali za awali na kutabiri majibu ya mawimbi. Katika tathmini ya mwisho, utafutaji wa hiperparameta unaotegemea data haukufanywa; muundo uleule uliowekwa mapema ulihifadhiwa bila kubadilishwa kwa aina zote za modeli, majibu na fold.
| Parameta | Thamani |
|---|---|
| n_estimators | 250 |
| learning_rate | 0,04 |
| max_depth | 2 |
| min_samples_leaf | 5 |
| subsample | 0,85 |
Modeli zote zilitathminiwa kwa uthibitishaji mtambuka wa fold tano. Katika modeli ya mfululizo, ili kuzuia utabiri wa hali ya kuingia ya kati kuvuja taarifa kwenye mafunzo, utaratibu wa ziada wa kujenga upya wa out-of-fold wa fold tano ulitumika ndani ya fold ya nje ya mafunzo.
Ujenzi upya wa hali ya kuingia majini
| Kigezo | \(R^2\) | MAE | RMSE |
|---|---|---|---|
| \(Fr\) | 0,843 | 0,038 | 0,050 |
| \(S\) | 0,865 | 0,008 | 0,010 |
| \(M\) | 0,824 | 0,018 | 0,024 |
Matokeo haya yanaonyesha kwamba, chini ya hali za majaribio zilizochunguzwa, hali halisi ya kuingia inaweza kujengwa upya kwa sehemu kwa utaratibu kutoka taarifa za awali. Hivyo \(X_e\) inaweza kutumika kama tabaka la kati lililo wazi katika modeli ya mfululizo.
Je, modeli ya mfululizo ni modeli sahihi zaidi ya utabiri?
Katika utafiti huu, hapana: usahihi wa modeli ya mfululizo uko karibu na kiwango kilekile cha modeli ya moja kwa moja inayotumia hali za awali pekee; faida kuu ya muundo wa mfululizo si usahihi wa juu zaidi, bali ni kuwakilisha hali ya kuingia majini kama tabaka la kati lililo wazi, kuimarisha tafsiri ya kifizikia na kutenganisha jinsi kosa la kujenga upya hali ya kuingia linavyoenezwa hadi kwenye utabiri wa mwisho wa wimbi.
Ulinganisho wa mifumo mitatu ya utabiri
| Modeli | Lengo | \(R^2\) | MAE | RMSE |
|---|---|---|---|---|
| Moja kwa moja kutoka hali za awali | \(A\) | 0,822 ± 0,042 | 0,030 ± 0,003 | 0,041 ± 0,009 |
| Moja kwa moja kutoka hali za awali | \(H\) | 0,794 ± 0,045 | 0,049 ± 0,006 | 0,066 ± 0,012 |
| Hali ya kuingia iliyozingatiwa | \(A\) | 0,945 ± 0,008 | 0,017 ± 0,001 | 0,023 ± 0,003 |
| Hali ya kuingia iliyozingatiwa | \(H\) | 0,924 ± 0,015 | 0,031 ± 0,001 | 0,040 ± 0,002 |
| Mfululizo | \(A\) | 0,831 ± 0,032 | 0,029 ± 0,003 | 0,040 ± 0,007 |
| Mfululizo | \(H\) | 0,794 ± 0,034 | 0,051 ± 0,005 | 0,066 ± 0,010 |
Modeli ya hali ya kuingia iliyozingatiwa ni bora wazi kwa nambari; hata hivyo, hutumia thamani halisi zilizopimwa za \(Fr\), \(S\) na \(M\). Kwa hiyo modeli hii haishughulikii tatizo lilelile la utabiri kama modeli inayofanya kazi kwa taarifa za awali pekee kwa tukio ambalo bado halijatokea. Kwa sababu hiyo watafiti huiangalia kama modeli ya ulinganisho wa maelezo au uchunguzi.
Ulinganisho halisi wa kutazama mbele ni kati ya modeli ya moja kwa moja inayotumia hali za awali na modeli ya mfululizo ambayo pia huanza tu kutoka hali za awali na kuunda yenyewe hali ya kati. Matokeo ya modeli hizi mbili yako karibu.
Ramani ya utaratibu
Utafiti uliweka \(L\) katika thamani ya wastani wa kati na kuchora ramani, katika tambarare ya \(M_i\)–\(Bi\), ya njia ya upatanishi yenye mchango mkubwa zaidi wa modeli kwa thamani kamili juu ya \(H\). Mchango wa ndani wa mpatanishi:
\[ C_m=b_{mH}\widehat{Z}_m \]
na mpatanishi unaotawala:
\[ m^*=\operatorname*{arg\,max}_{m\in\{Fr,S,M\}}|C_m| \]
vilifafanuliwa.
Katika modeli na nafasi ya vigezo iliyochunguzwa, katika takriban %79,4 ya seli, njia ya masa inayopatanishwa kupitia \(M\) ina mchango mkubwa zaidi kwa thamani kamili. Katika takriban %20,6, njia ya inersia kupitia \(Fr\) ndiyo inayotawala. Njia ya \(S\) haifanyi eneo kubwa la kutawala.
Maeneo haya si mipaka kamili ya tawala za kifizikia katika mazingira halisi. Yanawakilisha tu nguvu jamaa za njia katika modeli ya sasa na nafasi ya \(M_i\)–\(Bi\) iliyochunguzwa.
Uimara wa bootstrap
Uimara wa njia zisizo za moja kwa moja ulitathminiwa kwa sampuli upya 1000 za bootstrap. Kipindi cha kujiamini cha %95:
\[ CI_{95\%}=[Q_{2.5\%},Q_{97.5\%}] \]
kilifafanuliwa, na njia ambazo hazikujumuisha sifuri zilikubaliwa kuwa thabiti kitakwimu.
| Njia | Wastani | %2,5 | %97,5 |
|---|---|---|---|
| \(M_i\rightarrow M\rightarrow A\) | 0,635 | 0,444 | 0,848 |
| \(M_i\rightarrow M\rightarrow H\) | 0,584 | 0,396 | 0,772 |
| \(Bi\rightarrow M\rightarrow A\) | −0,653 | −0,888 | −0,459 |
| \(Bi\rightarrow M\rightarrow H\) | −0,600 | −0,805 | −0,412 |
| \(Bi\rightarrow Fr\rightarrow H\) | −0,486 | −0,583 | −0,394 |
Kutokuwapo kwa sifuri katika vipindi hivi vya kujiamini kunaunga mkono kwamba njia kuu za masa na inersia zilibaki thabiti kwa mwelekeo na ishara katika sampuli upya za bootstrap.
Kikomo cha multicollinearity
Ili kuchunguza uwezekano wa kutafsiri vigawo vya njia kwa kujitegemea, variance inflation factor ilitumika:
\[ VIF_j=\frac{1}{1-R_j^2} \]
| Hatua | Kigezo | VIF |
|---|---|---|
| 1 | \(M_i\) | 1,62 |
| 1 | \(L\) | 1,18 |
| 1 | \(Bi\) | 1,72 |
| 2 | \(M_i\) | 8,66 |
| 2 | \(L\) | 6,34 |
| 2 | \(Bi\) | 11,72 |
| 2 | \(Fr\) | 6,16 |
| 2 | \(S\) | 11,53 |
| 2 | \(M\) | 10,74 |
Katika hatua ya kwanza, thamani za VIF ni ndogo. Katika hatua ya pili, baadhi ya thamani kuzidi 10 kunaonyesha kwamba vigezo vya awali na vya kuingia vinawakilisha sifa zinazohusiana za mchakato uleule wa kifizikia, na hivyo vigawo vya mtu mmoja mmoja si huru kabisa kutoka kwa kila kimoja. Kwa hiyo watafiti wanatilia mkazo zaidi mwelekeo wa njia, muundo wa jumla na ulinganifu wa kifizikia kuliko tofauti ndogo za vigawo.
Uenezaji wa kosa katika mfululizo
Kosa la kujenga upya hali ya kuingia:
\[ e_m=|\widehat{X}_m-X_m| \]
na kosa la utabiri wa wimbi:
\[ e_y=|\widehat{Y}_y-Y_y| \]
vilifafanuliwa.
| Kosa la kujenga upya | Kosa la utabiri wa \(A\) | Kosa la utabiri wa \(H\) |
|---|---|---|
| \(Fr\) | 0,24 | 0,26 |
| \(S\) | 0,22 | 0,12 |
| \(M\) | 0,42 | 0,35 |
Kosa la kujenga upya \(M\) linaonyesha uhusiano wenye nguvu zaidi na kosa la mwisho la utabiri kwa \(A\) na \(H\). Matokeo haya ni muhimu hasa kwa sababu \(M\) pia ni mojawapo ya vigezo vyenye nguvu zaidi katika hatua ya pili. Kwa hiyo, katika modeli za mfululizo za baadaye, kuboresha utabiri wa masa tendaji ya kuingia ni mojawapo ya malengo ya kipaumbele ya kupunguza makosa.
Uchambuzi wa ziada wa uimara
Katika sehemu ya ziada ya makala, matokeo ya modeli ya njia ya mstari yalijaribiwa kwa kutumia udhibiti wa Ridge, kuondoa vigezo vya upatanishi kimoja baada ya kingine, vipengele vya daraja la pili na mwingiliano wa jozi, makosa sanifu imara ya HC3, uchunguzi wa mabaki na ushawishi, pamoja na uchanganuzi wa unyeti kwa umbali wa Cook.
Modeli ya Ridge ilihifadhi ishara za njia 20 kati ya 21 za OLS. Mabadiliko pekee ya ishara yalitokea katika njia ya moja kwa moja ya \(Bi\rightarrow H\), ambayo tayari ilikuwa dhaifu na si muhimu kitakwimu. Kuondoa \(M\) au \(Fr\) kutoka kwenye modeli kulipunguza utendaji wa uthibitishaji mtambuka wa milinganyo ya \(A\) na \(H\), ilhali kuondoa \(S\) kulisababisha mabadiliko madogo sana.
Uchunguzi wa mabaki ulionyesha hasa katika milinganyo ya \(A\) na \(H\) kwamba kulikuwa na heteroskedasticity, mabaki yasiyo ya kawaida na uchunguzi wenye ushawishi. Hata hivyo, uchunguzi uliotiwa alama kwa Cook ulipoondolewa kwa muda na modeli kujengwa upya, mwelekeo na utawala wa njia kuu za \(M\) na \(Fr\) kwa ujumla ulihifadhiwa. Kwa hiyo matokeo yanaunga mkono muundo mkuu wa mchakato, lakini yanahitaji ukubwa halisi wa nambari wa njia dhaifu kufasiriwa kwa tahadhari.
Matokeo yanayoungwa mkono na utafiti
Utafiti unaunga mkono kwamba katika seti ya data ya maabara iliyochunguzwa, hali halisi ya kuingia majini ina uhusiano wa karibu zaidi na uundaji wa wimbi kuliko hali za awali. Hasa, masa tendaji ya kuingia ya uwiano \(M\) iliibuka kama kigezo cha upatanishi kinachohusiana na masa chenye nguvu zaidi kwa amplitudo ya wimbi na urefu wa wimbi. Nambari ya Froude ya kuingia \(Fr\), hasa kwa urefu wa wimbi, ni njia ya pili muhimu ya inersia.
Utafiti pia unaonyesha kwamba \(Bi\) ya juu inahusiana vibaya na \(Fr\) na \(M\); umbali mrefu zaidi wa kuteleza wa uwiano \(L\) unaonyesha njia hasi na \(S\) na \(M\) chini ya hali za sasa za majaribio. Tafsiri ya kifizikia ya waandishi ni kwamba njia ndefu ya mwendo inaweza kusababisha kutawanyika, kupungua unene na ugawaji upya wa ndani katika masa ya viskoplastiki, na hivyo kupunguza masa tendaji ya kuingia. Tafsiri hii ni halali tu kwa nyenzo, jiometri na masafa ya vigezo yaliyopo.
Matokeo yasiyoungwa mkono na utafiti
Vigawo vya njia na upatanishi si uthibitisho wa usababisho kwa maana kali. Maeneo yanayopatanishwa na masa na inersia kwenye ramani ya utaratibu si mipaka isiyobadilika ya tawala katika mazingira halisi. Thamani kubwa ya \(R^2\) ya modeli ya hali ya kuingia iliyozingatiwa haimaanishi kwamba usahihi huohuo unaweza kupatikana katika utabiri wa kabla ya tukio kwa kutumia hali za awali pekee.
Vigawo vilivyotengenezwa kwa kutumia Carbopol katika jiometri isiyobadilika ya mfereji haviwezi kutumika moja kwa moja kwa maporomoko ya asili ya miamba au udongo. Katika maporomoko ya asili, utofauti wa nyenzo, kuvunjika vipande, kuingizwa kwa nyenzo kutoka mazingira, shinikizo la vinyweleo, usatishaji, kutawanyika kwa pande tatu, topografia changamano na athari za skeli lazima zishughulikiwe kando.
Mfumo wa sasa haujengi upya mabadiliko ya kina ya muda ya uso huru, sehemu za kasi za ndani, shinikizo za ndani au hidrodinamiki changamano ya pande tatu. Kwa upande huu hauchukui nafasi ya maelezo ya kifizikia ya CFD yenye usahihi mkubwa au suluhisho zinazotegemea chembe.
Utafiti ulilenga tu amplitudo ya juu ya wimbi ya uwiano \(A\) na urefu wa juu wa wimbi \(H\). Kipindi cha wimbi, upunguzaji, kupanda pwani, shinikizo la msukumo na ujazo wa kufurika si matokeo ya modeli ya sasa.
Maana ya kisayansi
Mchango mkuu wa utafiti ni kubadilisha swali “ni kigezo gani cha awali kilicho muhimu zaidi?” kuwa swali la kifizikia zaidi: Mabadiliko ya awali hubadilishaje hali ya kuingia majini, na hali hiyo ya kuingia iliyobadilika huunganishwaje na uundaji wa wimbi kupitia njia gani?
Tofauti hii ni ya thamani hasa kwa maporomoko yanayobadilika umbo. Maporomoko mawili yenye masa sawa ya awali yanaweza yasibebe masa tendaji na inersia sawa yanapofika majini kutokana na upinzani tofauti wa reolojia, kutawanyika tofauti njiani au kasi tofauti za kuingia. Kwa hiyo masa sawa ya awali haihakikishi jibu sawa la wimbi.
Picha ya jumla yenye nguvu zaidi inayojitokeza katika seti ya data ni kwamba uhamishaji wa masa na uhamishaji wa inersia ni njia mbili kuu zinazobeba jibu la wimbi katika hatua ya kuingia majini. Unene wa kuingia wa uwiano \(S\) hubeba taarifa; lakini \(M\) na \(Fr\) zinapoingizwa pamoja katika modeli, mchango wake huru hupungua kwa kiasi kikubwa.
Maelezo ya Chanzo na Mbinu
Kichwa asilia: Process-Chain Analysis of Viscoplastic Landslide-Generated Impulse Waves from Initial Slide Conditions Through Water-Entry Dynamics to Wave Responses
Waandishi: Zhenxia Yuan; Yadong Bian; Song Yin; Xian Liu; Xinchen Pan; Zhenzhu Meng; Bohan Cheng.
Mwandishi wa mawasiliano: Zhenxia Yuan.
Taasisi: School of Intelligent Construction and Civil Engineering, Zhongyuan University of Technology, Zhengzhou, China; Henan Environmental Geotechnical Engineering and Underground Engineering Disaster Control Engineering Research Center, Zhengzhou, China; School of Civil Engineering, Henan Polytechnic University, Jiaozuo, China; College of Civil Engineering, Sun Yat-Sen University, Zhuhai, China; School of Hydraulic Engineering, Zhejiang University of Water Resources and Electric Power, Hangzhou, China.
Jarida: Water.
Mchapishaji: MDPI.
Rekodi ya bibliografia: Water 2026, 18(17), 2189.
DOI: 10.3390/w18172189
Aina ya chanzo: Makala ya utafiti iliyopitiwa na wataalamu.
Tarehe ya kupokelewa: 3 Julai 2026.
Tarehe ya marekebisho: 5 Agosti 2026.
Tarehe ya kukubaliwa: 31 Agosti 2026.
Tarehe ya kuchapishwa: 3 Septemba 2026.
Leseni: Creative Commons Attribution (CC BY).
Ufadhili: Henan Provincial Key Research and Development Program, Grant No. 251111231500; Central Guidance Funds for Science and Technology Local Development Projects, Grant No. 2025ZY01091.
Upatikanaji wa data: Kwa sababu data zilizotumika katika utafiti zinaendelea kutumika katika utafiti wa ufuatiliaji, imeelezwa kuwa seti ya data inaweza kupatikana kutoka kwa mwandishi wa mawasiliano kwa ombi la msingi.
Mgongano wa maslahi: Waandishi wametangaza kuwa hakuna mgongano wa maslahi.
Michango ya waandishi: Zhenxia Yuan na Yadong Bian — uundaji wa dhana na mbinu; Song Yin, Xian Liu na Xinchen Pan — uthibitishaji; Zhenxia Yuan na Bohan Cheng — uchanganuzi rasmi; Zhenzhu Meng — uratibu wa data; Zhenxia Yuan — rasimu ya kwanza na mapitio/uhariri; Yadong Bian — usimamizi na usimamizi wa mradi; Zhenxia Yuan na Yadong Bian — upatikanaji wa ufadhili.
Mazingira ya ukokotoaji: Uchambuzi ulifanywa kwa Python 3.14.6, kwa kutumia NumPy 2.5.1, pandas 3.0.5, SciPy 1.18.0, scikit-learn 1.9.0, statsmodels 0.14.6 na Matplotlib 3.11.1. Modeli za njia zilijengwa kwa ordinary least squares baada ya uwekaji sanifu wa z-skor; katika uchanganuzi wa upatanishi, sampuli upya 1000 za bootstrap zilitumika kwa random seed 42.
Vikwazo vikuu vya mbinu: Majaribio yalifanywa kwa Carbopol, jiometri isiyobadilika ya mfereji na masafa machache ya vigezo katika skeli ya maabara. Vigawo vya modeli vinahitaji kukalibishwa upya kwa nyenzo, jiometri na skeli tofauti. Utofauti, kuvunjika vipande, entrainment, shinikizo la vinyweleo, usatishaji, topografia ya pande tatu na athari za skeli katika maporomoko ya asili hazikuundwa moja kwa moja katika modeli. Ufafanuzi wa masa tendaji ya kuingia unaotegemea wakati wa urefu wa juu wa wimbi unaondoa hadhi yake kama kigezo huru cha upatanishi kabla ya tukio kwa maana kali. Modeli ya njia ni ya mstari na ya kujumlisha; thamani za juu za VIF katika hatua ya pili zinaweka kikomo katika tafsiri ya vigawo binafsi. Utafiti unalenga \(A\) na \(H\) pekee.
Acha maoni
Anwani yako ya barua pepe haitachapishwa. Sehemu za lazima zimewekewa alama ya *