Көчкү жараткан импульстук толкун — ылдам кыймылдаган көчкү материалы суу массасына кирип, сууга масса жана импульс өткөргөндө пайда болгон беттик толкун. Бул изилдөө вископластикалык көчкүнүн баштапкы шарттары толкун амплитудасына жана толкун бийиктигине кантип айланарын башталыш менен натыйжанын ортосундагы түз байланыш катары гана эмес, алардын ортосундагы чыныгы сууга кирүү абалын ачык чагылдырган үч баскычтуу процесстер чынжыры катары карайт. Жалпысынан 290 лабораториялык эксперимент колдонулуп, баштапкы көчкү шарттары, иш жүзүндө болгон сууга кирүү динамикасы жана толкун жооптору өзгөрмөлөрдүн өзүнчө катмарлары катары талданган. Негизги жыйынтык — толкун жооптору баштапкы шарттарга караганда сууга чындап кирген материалдын массасы жана кирүү инерциясы менен тыгызыраак байланыштуу. Бирок жыйынтыктар белгилүү Carbopol материалына, канал геометриясына жана эксперименттик параметрлер диапазонуна тиешелүү; статистикалык жолдор себептүүлүктүн далили эмес жана табигый көчкүлөргө түз колдонуу үчүн көз карандысыз талаалык текшерүү талап кылынат.
Изилдөөдө баштапкы шарттар \(X_0=(M_i,L,Bi)\), иш жүзүндө болгон сууга кирүү абалы \(X_e=(Fr,S,M)\) жана толкун жообу \(Y=(A,H)\) түрүндө уюштурулган. Бул жерде \(M_i\) салыштырмалуу баштапкы массаны, \(L\) салыштырмалуу жылышуу аралыкты, \(Bi\) Bingham санын, \(Fr\) кирүү Froude санын, \(S\) салыштырмалуу кирүү калыңдыгын, \(M\) салыштырмалуу натыйжалуу кирүү массасын, \(A\) салыштырмалуу максималдуу толкун амплитудасын жана \(H\) салыштырмалуу максималдуу толкун бийиктигин билдирет. Ошентип, борбордук илимий түзүлүш \(X_0 \rightarrow X_e \rightarrow Y\) процесстер чынжыры болуп саналат.
Салыштырмалуу натыйжалуу кирүү массасы \(M\), \(A\) менен 0,96 жана \(H\) менен 0,91 Spearman корреляциясын көрсөткөн. Кирүү Froude саны \(Fr\) үчүн тиешелүү корреляциялар 0,71 жана 0,80 болгон. Медиатордук талдоодо баштапкы салыштырмалуу масса \(M_i\) нин \(M\) аркылуу кыйыр компоненти \(A\) үчүн 0,652, \(H\) үчүн 0,600; Bingham санынын \(M\) аркылуу кыйыр компоненти болсо тиешелүүлүгүнө жараша −0,669 жана −0,615 деп табылган. Бул маанилер каралган маалымат топтому жана курулган модель чегинде масса менен байланышкан сууга кирүү жолунун үстөм ролун көрсөтөт.
Gradient Boosting негизиндеги каскаддык ыкма баштапкы шарттардан алгач \(Fr\), \(S\) жана \(M\) ни кайра түзүп, андан кийин толкун жообун божомолдогон. Бирок каскаддык модель баштапкы шарттарды гана колдонгон түз моделге салыштырмалуу чоң тактык өсүшүн берген эмес. Демек, каскаддык алкактын негизги артыкчылыгы божомол тактыгы эмес; баштапкы шарттар, сууга кирүү жана толкун пайда болуусунун ортосундагы физикалык жана статистикалык аралык жолдорду көрүнүктүү кылуу жана ката булактарын баскычтарга бөлүп кароо мүмкүнчүлүгү.
Негизги илимий маселе
Көчкү сууга киргенде пайда болгон толкун көчкүнүн баштапкы массасына же баштапкы абалына гана көз каранды эмес. Көчкү жантайыңкы бет боюнча кыймылдаганда ылдамдашы, жайылышы, жукарышы же калыңдашы жана ички деформацияга учурашы мүмкүн. Ошондуктан башында белгиленген масса, калыңдык жана реологиялык абал суу бетине жеткен учурдагы чыныгы кыймыл абалынан айырмаланышы мүмкүн.
Айрыкча деформациялана алган вископластикалык материалдарда бул айырмачылык маанилүү. Катуу блок сууга болжол менен бир бүтүн геометриялык дене катары кириши мүмкүн болсо, вископластикалык масса жантайыңкы бет боюнча кыймыл учурунда формасын өзгөртөт. Натыйжада баштапкы массанын баары толкун пайда болуусунун чечүүчү учурунда сууну жылдырууга бирдей деңгээлде катышпашы мүмкүн.
Мурунку маалыматка негизделген ыкмалардын маанилүү бөлүгү баштапкы өзгөрмөлөрдү түздөн-түз толкун касиеттери менен байланыштырып же өзгөрмөлөрдүн маанилүүлүк рейтингин түзгөн. Мындай ыкма божомол үчүн пайдалуу болушу мүмкүн; бирок баштапкы шарт толкунга түз таасир этеби же анын таасири адегенде сууга кирүү массасын, кирүү ылдамдыгын же кирүү геометриясын өзгөртүү аркылуу пайда болобу — муну ажыратууда чектелүү.
Бул изилдөөнүн салымы жаңы статистикалык же машиналык үйрөнүү алгоритмин иштеп чыгуу эмес. Изилдөөнүн жаңылыгы — учурдагы корреляция, жол талдоосу, медиатордук талдоо жана Gradient Boosting ыкмаларын физикалык окуялардын иретине ылайык уюштурулган бирдиктүү процесстер чынжырында бириктирүүсү.
Көчкү толкундарында сууга кирүү абалы эмне үчүн өзүнчө баскыч катары каралат?
Анткени вископластикалык көчкү материалы бошотулган баштапкы абалдан сууга жеткенге чейин деформацияланып, жайылып, жукарып же калыңдап жана башка ылдамдыкка ээ болушу мүмкүн; демек, толкунду түз аныктаган масса, геометрия жана инерция абалы башында берилген шарттар менен сөзсүз бирдей болбойт.
Эксперименттик физикалык система
Лабораториялык система эки өлчөмдүү көчкү–толкун каналы. Негизги компоненттер — көчкү материалын кармап турган резервуар, пневматикалык ачылма дарбаза, жантайышы жөнгө салынуучу PVC нуа, тар тунук толкун каналы жана жогорку ылдамдыктагы сүрөткө тартуу системасы.
| Элемент | Булакта берилген мүнөздөмө |
|---|---|
| Толкун каналынын узундугу | 2,5 m |
| Толкун каналынын туурасы \(B\) | 0,12 m |
| Нуанын жөнгө салынуучу жантайышы | 30–50° |
| Эксперименттерде колдонулган жантайыш \(\theta\) | 40–45° |
| Туруктуу суунун тереңдиги \(h_0\) | 0,20 m |
| Камера | 400 fps |
| Сүрөттүн чечилиши | 600 × 800 пиксел |
| Жалпы эксперименттик учурлар | 290 |
Пневматикалык дарбаза ачылганда Carbopol материалы тартылуу күчүнүн таасири менен нуа боюнча ылдый жылып, сууга кирет. Сууга кирүү айланадагы сууну жылдырып жана сууга импульс өткөрүп, алгачкы эркин бет бузулуусун түзөт; бул бузулуу кийин импульстук толкунга айланат.
Сүрөткө тартуу аймагы кирүү зонасын да, жакынкы аймактагы толкун пайда болуусун да камтыган. Экспериментке чейин өлчөмдөрү белгилүү калибрлөө тору колдонулуп, пиксел координаттары физикалык аралыктарга айландырылган. Камеранын орду жана көрүү талаасы эксперименттер бою өзгөртүлгөн эмес.
Толкун амплитудасы \(a_m\) эркин беттин тынч суунун деңгээлине салыштырмалуу максималдуу өйдө жылышы катары аныкталган. Максималдуу толкун бийиктиги \(h_m\) болсо тиешелүү толкун чокусу менен ойдуңунун ортосундагы максималдуу тик аралык. Бул эки чоңдук окшош көрүнгөнү менен бирдей физикалык өлчөө эмес.
Натыйжалуу кирүү массасы деген эмне?
Натыйжалуу кирүү массасы \(m_e\) — максималдуу толкун бийиктиги \(h_m\) пайда болгон учурда суу массасына кирип калган көчкү материалынын массасы. Эсептөө сүрөттөрдөн алынган көчкү контуруна жана материалдын белгилүү тыгыздыгына негизделет.
Бул аныктама баштапкы массадан айырмаланат. Баштапкы масса \(m_i\) эксперимент башталганга чейин болгон жалпы материалды билдирсе, \(m_e\) негизги толкун пайда болуусунун чечүүчү учурунда сууну жылдырууга чындап катышкан материалды билдирет.
Натыйжалуу кирүү массасы эмне үчүн толкун жообунда мынчалык маанилүү?
Салыштырмалуу натыйжалуу кирүү массасы \(M\) көчкүнүн сууга жеткенге чейинки деформациясынын, жайылышынын, реологиялык каршылыгынын жана убакыттагы кайра бөлүштүрүлүшүнүн бириккен натыйжасын алып жүргөндүктөн, баштапкы массага караганда негизги толкун пайда болуу баскычында сууну канча материал жылдырганын түзүрөөк чагылдырат; бул изилдөөдө \(M\) салыштырмалуу толкун амплитудасы \(A\) жана салыштырмалуу толкун бийиктиги \(H\) менен эң күчтүү кирүү-абалы байланышын көрсөткөн.
Carbopol эмне үчүн колдонулган?
Carbopol — ачык агуу чыңалуусуна жана кесүүдө суюлуу жүрүм-турумуна ээ идеалдаштырылган вископластикалык материал. Агуудан кийинки реологиялык жүрүм-туруму Herschel–Bulkley модели менен берилген:
\[ \tau=\tau_c+K\dot{\gamma}^{\,n}, \qquad \tau>\tau_c \]
Бул жерде \(\tau\) кесүү чыңалуусу, \(\tau_c\) агуу чыңалуусу, \(K\) консистенция коэффициенти, \(\dot{\gamma}\) кесүү ылдамдыгы жана \(n\) агым жүрүм-турумунун индекси. Колдонулган кесүү чыңалуусу \(\tau_c\) ден ашпаса, материал үзгүлтүксүз агымга өтпөйт. Ошондуктан \(\tau_c\) үзгүлтүксүз деформация жана агым баштала электе ашып өтүлүшү керек болгон босого каршылыкты билдирет.
Изилдөө \(\tau_c\) ни көчкү материалынын реологиялык жүрүм-турумун билдирген негизги өзгөрмө катары колдонуп, анын баштапкы абалдан сууга кирүү инерциясына жана натыйжалуу кирүү массасына чейинки байланыштарын Bingham саны аркылуу изилдеген.
Эксперименттик параметрлердин диапазондору
| Өзгөрмө | Мааниси | Диапазон / маани |
|---|---|---|
| \(\tau_c\) | Материалдын агуу чыңалуусу | 60–90 Pa |
| \(m_i\) | Баштапкы көчкү массасы | 2,0–5,6 kg |
| \(l_s\) | Жылышуу аралыгы | 0,55–0,75 m |
| \(\theta\) | Жантайыш бурчу | 40–45° |
| \(s_0\) | Суу бетиндеги кирүү калыңдыгы | 0,026–0,050 m |
| \(u_0\) | Кирүү ылдамдыгы | 1,12–2,36 m/s |
| \(h_0\) | Тынч суунун тереңдиги | 0,20 m |
Өлчөө белгисиздиги калибрлөө торунан алынган мейкиндик чечилиши жана камеранын убакыттык чечилиши менен бааланган. Мындан тышкары, өкүлчүлүктүү бир эксперименттик шарт ошол эле жөндөөлөрдө үч жолу көз карандысыз кайталанып, кайталанган өлчөөлөрдөгү өзгөрмөлүүлүк эксперименттер аралык кайталануучулук белгисиздигин баалоо үчүн колдонулган.
Өлчөмсүз өзгөрмөлөр системасынын максаты
Ар башка физикалык чоңдуктардагы эксперименттик шарттарды жалпы математикалык алкакта салыштыруу үчүн физикалык өзгөрмөлөр өлчөмсүздөштүрүлгөн.
Баштапкы шарттар:
\[ 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\) — баштапкы массанын суунун тыгыздыгына, каналдын туурасына жана суунун тереңдигине карата масштабдалган түрү. \(L\) жылышуу аралыкты суунун тереңдиги менен масштабдайт. \(Bi\) болсо материалдын агуу каршылыгынын баштапкы конфигурациядагы жантайыш боюнча гравитациялык айдоочу чыңалууга катышын билдирет.
Иш жүзүндө болгон сууга кирүү абалы:
\[ 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\) — сууга кирүү ылдамдыгынын тартылуу күчү жана суунун тереңдиги менен аныкталган мүнөздүү толкун ылдамдыгына карата өлчөмсүз өлчөмү. \(S\) кирүү калыңдыгын суунун тереңдигине салыштырат. \(M\) — толкундун негизги пайда болуу баскычына натыйжалуу катышкан массанын өлчөмсүз түрү.
Толкун жооптору:
\[ A=\frac{a_m}{h_0}, \qquad H=\frac{h_m}{h_0} \]
\(A\) эркин беттин максималдуу өйдө четтөөсүнүн; \(H\) болсо чоку–ойдуң толкун бийиктигинин тынч суунун тереңдигине катышы.
Үч баскычтуу процесстер чынжыры
| 1-баскыч | 2-баскыч | 3-баскыч | ||
|---|---|---|---|---|
| Баштапкы шарттар \(X_0=(M_i,L,Bi)\) | → | Иш жүзүндө болгон сууга кирүү абалы \(X_e=(Fr,S,M)\) | → | Толкун жообу \(Y=(A,H)\) |
Бул түзүлүш баштапкы шарттардан толкун жообуна түз байланыштарды изилдөөнү улантуу менен бирге, ошол байланыштын канча бөлүгү \(Fr\), \(S\) жана \(M\) аркылуу берилерин өз-өзүнчө баалоого мүмкүндүк берет.
Стандартташтыруу жана жол модели
Жол талдоосуна чейин бардык өзгөрмөлөр z-упай менен стандартташтырылган:
\[ Z_i=\frac{X_i-\mu_i}{\sigma_i} \]
Бул жерде \(X_i\) баштапкы өзгөрмө, \(\mu_i\) орточо маани, \(\sigma_i\) стандарттык четтөө жана \(Z_i\) стандартташтырылган маани. Бул операция ар башка масштабдагы өзгөрмөлөрдүн коэффициенттерин салыштырууга мүмкүндүк берет.
Биринчи баскычтагы жол баштапкы өзгөрмөлөрдөн ар бир сууга кирүү өзгөрмөсүнө курулган:
\[ Z_m=\alpha_{0m}+\sum_{k=1}^{3}a_{km}Z_k+\varepsilon_m \]
Бул жерде \(Z_k\), \(M_i\), \(L\) же \(Bi\) ни; \(Z_m\), \(Fr\), \(S\) же \(M\) ни; \(a_{km}\) тиешелүү стандартташтырылган жол коэффициентин билдирет.
Экинчи баскыч модели баштапкы өзгөрмөлөр менен сууга кирүү өзгөрмөлөрүн бирге толкун жообу менен байланыштырат:
\[ Z_y=\beta_{0y}+\sum_{k=1}^{3}c'_{ky}Z_k+\sum_{m=1}^{3}b_{my}Z_m+\varepsilon_y \]
Бул жердеги \(c'_{ky}\) кирүү өзгөрмөлөрү эсепке алынгандан кийин баштапкы өзгөрмөдөн толкун жообуна калган түз шарттуу байланышты; \(b_{my}\) болсо кирүү өзгөрмөсүнөн толкун жообуна барган жолду билдирет.
Белгилүү бир медиатор аркылуу кыйыр компонент:
\[ IE_{kmy}=a_{km}b_{my} \]
Жалпы кыйыр компонент:
\[ TIE_{ky}=\sum_{m=1}^{3}IE_{kmy} \]
Жалпы компонент болсо:
\[ TE_{ky}=c'_{ky}+TIE_{ky} \]
Булак изилдөөнүн терминологиясында булар түз, кыйыр жана жалпы “таасирлер” деп аталат; бирок эксперименттик убакыт ирети менен колдоого алынганы менен, курулган сызыктуу модель ичиндеги шарттуу статистикалык байланыштар. Өздөрү эле себептүүлүктү далилдебейт.
Каскаддык божомол логикасы
Иш жүзүндө болгон сууга кирүү абалын баштапкы шарттардан гана божомолдоого болор-болбосун текшерүү үчүн:
\[ \widehat{X}_e=f_1(X_0) \]
жана ачык түрдө:
\[ (\widehat{Fr},\widehat{S},\widehat{M})=f_1(M_i,L,Bi) \]
түрүндө баштапкыдан-кирүүгө модели курулган.
Андан кийин үч толкун божомолдоо алкагы салыштырылган:
| Модель | Киргизүү | Илимий колдонуу |
|---|---|---|
| Түз баштапкы-шарт модели | \(X_0\) | Башында белгилүү маалыматтардан алдыга карай божомол |
| Байкалган кирүү-абалы модели | Өлчөнгөн \(X_e\) | Түшүндүрмө / диагностикалык салыштыруу |
| Каскаддык модель | \(X_0\rightarrow\widehat{X}_e\rightarrow\widehat{Y}\) | Баштапкы шарттан алдыга карай божомол + аралык процесстин чечмелениши |
Каскаддык моделдин математикалык жыйынтыгы:
\[ \widehat{Y}=f_2[f_1(X_0)] \]
түрүндө.
Изилдөөнүн Ыкмасы жана Жыйынтыктары
Маалымат топтому жана негизги бөлүштүрүүлөр
| Өзгөрмө | Орточо | Стандарттык четтөө | Минимум | Максимум |
|---|---|---|---|---|
| \(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 |
Баскычтар аралык корреляциялар
Баштапкы салыштырмалуу масса \(M_i\), кирүү калыңдыгы \(S\) менен 0,61 жана натыйжалуу кирүү массасы \(M\) менен 0,36 оң Spearman корреляциясын көрсөткөн. Bingham саны \(Bi\), кирүү Froude саны \(Fr\) менен −0,79 жана \(M\) менен −0,29 корреляция көрсөткөн. Салыштырмалуу жылышуу аралыгы \(L\), \(Fr\) менен 0,33 оң; \(S\) жана \(M\) менен тиешелүүлүгүнө жараша −0,60 жана −0,26 терс корреляция көрсөткөн.
| Кирүү өзгөрмөсү | \(A\) менен корреляция | \(H\) менен корреляция |
|---|---|---|
| \(Fr\) | 0,71 | 0,80 |
| \(S\) | 0,55 | 0,42 |
| \(M\) | 0,96 | 0,91 |
Бул салыштырууда \(M\) эки толкун өзгөрмөсү менен тең эң күчтүү байланыш көрсөткөн. \(Fr\) экинчи күчтүү өзгөрмө. \(S\) корреляциясы оң болгону менен төмөнүрөөк.
Баштапкы шарттар толкун жооптору менен түз салыштырылганда байланыштар жалпысынан алсызыраак. Бул топто эң күчтүү баштапкы өзгөрмө \(Bi\) болуп, \(A\) үчүн −0,44 жана \(H\) үчүн −0,56 корреляция көрсөткөн.
Bingham саны толкун пайда болуусуна кандай таасир этет?
Бул маалымат топтомунда жогорураак \(Bi\) материалдын агуу каршылыгынын жантайыш боюнча айдоочу чыңалууга салыштырмалуу чоңураак экенин билдирет жана кирүү Froude саны \(Fr\) жана натыйжалуу кирүү массасы \(M\) менен күчтүү терс жолдорду көрсөтөт; демек, жогорку реологиялык каршылык модель ичинде төмөнүрөөк кирүү инерциясы жана толкун жаратууга натыйжалуу катышкан азыраак материал менен байланыштуу.
Стандартташтырылган жол коэффициенттери
| Баштапкы өзгөрмө | \(Fr\) | \(S\) | \(M\) |
|---|---|---|---|
| \(M_i\) | 0,521 | 0,839 | 1,006 |
| \(L\) | 0,033; маанилүү эмес | −0,741 | −0,622 |
| \(Bi\) | −1,061 | −0,261 | −1,031 |
Баштапкы масса көбөйгөн сайын үч сууга кирүү өзгөрмөсү тең өсүү тенденциясын көрсөтөт. Тескерисинче, жогорку \(Bi\), айрыкча \(Fr\) жана \(M\) ни күчтүү азайтуучу шарттуу байланыштарды көрсөтөт. Жылышуу аралыгынын \(Fr\) менен көз карандысыз жолу кичине жана статистикалык жактан маанилүү эмес, ал эми \(S\) жана \(M\) менен жолдору ачык терс.
| Медиатор | \(A\) | \(H\) |
|---|---|---|
| \(M\) | 0,649 | 0,596 |
| \(Fr\) | 0,316 | 0,458 |
| \(S\) | 0,072 | 0,100 |
\(M\) жана \(Fr\) эки толкун жообуна барган үстөм оң коэффициенттерди алып жүрөт. \(S\) нин экинчи баскычтагы коэффициенттери салыштырмалуу кичине жана булак талдоосунда статистикалык маанилүүлүккө жеткен эмес. Кирүү абалы эсепке алынгандан кийин баштапкы өзгөрмөлөрдүн түз жолдору да көбүнчө кичине; алардын ичинен \(L\rightarrow A\) жолу гана −0,114 менен статистикалык жактан маанилүү.
Медиатордук жолдор
| Баштапкы өзгөрмө | Жол | \(A\) | \(H\) |
|---|---|---|---|
| \(M_i\) | \(Fr\) аркылуу | 0,165 | 0,239 |
| \(M_i\) | \(S\) аркылуу | 0,060 | 0,084 |
| \(M_i\) | \(M\) аркылуу | 0,652 | 0,600 |
| \(L\) | \(Fr\) аркылуу | 0,010 | 0,015 |
| \(L\) | \(S\) аркылуу | −0,053 | −0,074 |
| \(L\) | \(M\) аркылуу | −0,403 | −0,371 |
| \(Bi\) | \(Fr\) аркылуу | −0,335 | −0,486 |
| \(Bi\) | \(S\) аркылуу | −0,019 | −0,026 |
| \(Bi\) | \(M\) аркылуу | −0,669 | −0,615 |
\(M_i\) нин \(A\) га жалпы стандартташтырылган компоненти 0,945, \(H\) га жалпы компоненти 0,855. \(Bi\) үчүн булар −1,098 жана −1,117; \(L\) үчүн −0,560 жана −0,471. Үч баштапкы өзгөрмөнүн тең эң чоң кыйыр компоненттеринин көбү \(M\) аркылуу пайда болгон.
Gradient Boosting жана кайчылаш текшерүү
Сууга кирүү өзгөрмөлөрүн баштапкы шарттардан кайра түзүү жана толкун жоопторун божомолдоо үчүн Gradient Boosting Regression колдонулган. Акыркы баалоодо маалыматка негизделген гиперпараметр издөө жүргүзүлгөн эмес; алдын ала белгиленген бир эле түзүлүш бардык модель түрлөрү, жооптор жана катмарлар үчүн туруктуу сакталган.
| Параметр | Маани |
|---|---|
| n_estimators | 250 |
| learning_rate | 0,04 |
| max_depth | 2 |
| min_samples_leaf | 5 |
| subsample | 0,85 |
Бардык моделдер беш катмарлуу кайчылаш текшерүү менен бааланган. Каскаддык моделде аралык кирүү-абалы божомолдору окутууга маалымат агып кетпеши үчүн тышкы окутуу катмарынын ичинде кошумча беш катмарлуу out-of-fold кайра түзүү процедурасы колдонулган.
Сууга кирүү абалын кайра түзүү
| Өзгөрмө | \(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 |
Бул жыйынтыктар каралган эксперименттик шарттарда иш жүзүндө болгон кирүү абалын баштапкы маалыматтардан системалуу түрдө жарым-жартылай кайра түзүүгө болорун көрсөтөт. Ошентип, \(X_e\) каскаддык моделде ачык аралык катмар катары колдонулушу мүмкүн.
Каскаддык модель тагыраак божомол модели болобу?
Бул изилдөөдө жок: каскаддык моделдин тактыгы баштапкы шарттарды гана колдонгон түз модель менен болжол менен бир деңгээлде; каскаддык түзүлүштүн негизги пайдасы жогорку тактык эмес, сууга кирүү абалын ачык аралык катмар катары көрсөтүү, физикалык чечмелөөнү күчөтүү жана кирүү-абалын кайра түзүү катасы акыркы толкун божомолуна кантип өтөрүн ажыратуу мүмкүнчүлүгү.
Үч божомолдоо алкагын салыштыруу
| Модель | Максат | \(R^2\) | MAE | RMSE |
|---|---|---|---|---|
| Түз баштапкы-шарт | \(A\) | 0,822 ± 0,042 | 0,030 ± 0,003 | 0,041 ± 0,009 |
| Түз баштапкы-шарт | \(H\) | 0,794 ± 0,045 | 0,049 ± 0,006 | 0,066 ± 0,012 |
| Байкалган кирүү-абалы | \(A\) | 0,945 ± 0,008 | 0,017 ± 0,001 | 0,023 ± 0,003 |
| Байкалган кирүү-абалы | \(H\) | 0,924 ± 0,015 | 0,031 ± 0,001 | 0,040 ± 0,002 |
| Каскад | \(A\) | 0,831 ± 0,032 | 0,029 ± 0,003 | 0,040 ± 0,007 |
| Каскад | \(H\) | 0,794 ± 0,034 | 0,051 ± 0,005 | 0,066 ± 0,010 |
Байкалган кирүү-абалы модели сандык жактан ачык эле ийгиликтүүрөөк; бирок ал \(Fr\), \(S\) жана \(M\) нин чыныгы өлчөнгөн маанилерин колдонот. Демек, бул модель боло элек окуя үчүн баштапкы маалыматтар менен гана иштеген модель менен бир эле божомол маселесинде эмес. Ошондуктан изилдөөчүлөр аны түшүндүрмө же диагностикалык салыштыруу модели катары баалашат.
Чыныгы алдыга карай салыштыруу — баштапкы шарттарды колдонгон түз модель менен кайра эле баштапкы шарттардан баштап аралык абалды өзү түзгөн каскаддык модель ортосунда. Бул эки моделдин жыйынтыктары бири-бирине жакын.
Механизм картасы
Изилдөө \(L\) ни медианалык мааниде туруктуу кармап, \(M_i\)–\(Bi\) тегиздигинде кайсы медиатордук жолдун \(H\) га абсолюттук модель салымы эң чоң экенин карталаган. Жергиликтүү медиатор салымы:
\[ C_m=b_{mH}\widehat{Z}_m \]
жана үстөм медиатор:
\[ m^*=\operatorname*{arg\,max}_{m\in\{Fr,S,M\}}|C_m| \]
деп аныкталган.
Каралган модель жана параметр мейкиндигинде уячалардын болжол менен %79,4 ында \(M\) аркылуу масса-медиатордук жол эң чоң абсолюттук салымга ээ. Болжол менен %20,6 ында \(Fr\) аркылуу инерция жолу үстөм. \(S\) жолу кең үстөм аймак түзбөйт.
Бул аймактар чыныгы табияттагы так физикалык режим чек аралары эмес. Алар болгону учурдагы моделдин жана каралган \(M_i\)–\(Bi\) мейкиндигинин салыштырмалуу жол күчтөрүн билдирет.
Bootstrap туруктуулугу
Кыйыр жолдордун туруктуулугу 1000 bootstrap кайра үлгүлөө менен бааланган. %95 ишеним аралыгы:
\[ CI_{95\%}=[Q_{2.5\%},Q_{97.5\%}] \]
деп аныкталып, нөлдү камтыбаган жолдор статистикалык жактан туруктуу деп кабыл алынган.
| Жол | Орточо | %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 |
Бул ишеним аралыктарынын нөлдү камтыбашы негизги масса жана инерция жолдору bootstrap кайра үлгүлөөлөрүндө багыты жана белгиси боюнча туруктуу калганын колдойт.
Көп коллинеардуулук чеги
Жол коэффициенттеринин өз алдынча чечмеленишин изилдөө үчүн дисперсияны инфляциялоо фактору колдонулган:
\[ VIF_j=\frac{1}{1-R_j^2} \]
| Баскыч | Өзгөрмө | 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 |
Биринчи баскычта VIF маанилери төмөн. Экинчи баскычта айрым маанилердин 10дон ашышы баштапкы жана кирүү өзгөрмөлөрү ошол эле физикалык процесстин байланышкан касиеттерин билдирерин, демек жеке коэффициенттер бири-биринен толук көз карандысыз эместигин көрсөтөт. Ошондуктан изилдөөчүлөр майда коэффициент айырмачылыктарына караганда жолдордун багытына, жалпы түзүлүшүнө жана физикалык шайкештигине көбүрөөк көңүл бурушат.
Каскаддык катанын таралышы
Кирүү абалын кайра түзүү катасы:
\[ e_m=|\widehat{X}_m-X_m| \]
жана толкун божомол катасы:
\[ e_y=|\widehat{Y}_y-Y_y| \]
деп аныкталган.
| Кайра түзүү катасы | \(A\) божомол катасы | \(H\) божомол катасы |
|---|---|---|
| \(Fr\) | 0,24 | 0,26 |
| \(S\) | 0,22 | 0,12 |
| \(M\) | 0,42 | 0,35 |
\(M\) ни кайра түзүү катасы \(A\) жана \(H\) үчүн акыркы божомол катасы менен эң күчтүү байланыш көрсөткөн. Бул жыйынтык өзгөчө маанилүү, анткени \(M\) ошол эле учурда экинчи баскычтагы эң күчтүү өзгөрмөлөрдүн бири. Ошондуктан келечектеги каскаддык моделдерде натыйжалуу кирүү массасын жакшыраак божомолдоо катаны азайтуу үчүн артыкчылыктуу максаттардын бири.
Кошумча туруктуулук талдоолору
Макаланын кошумча бөлүгүндө сызыктуу жол моделинин жыйынтыктары Ridge регуляризациясы, медиатор өзгөрмөлөрдү бирден алып салуу, экинчи даражадагы жана жуптук өз ара аракеттенүү мүчөлөрү, HC3 туруктуу стандарттык каталары, калдык жана таасир диагностикасы жана Cook аралыкка сезгичтик талдоолору менен текшерилген.
Ridge модели OLS жолдорунун 21инин 20сынын белгисин сактаган. Жалгыз белги өзгөрүшү ансыз да алсыз жана статистикалык жактан маанилүү эмес \(Bi\rightarrow H\) түз жолунда болгон. \(M\) же \(Fr\) моделден чыгарылганда \(A\) жана \(H\) теңдемелеринин кайчылаш текшерүү көрсөткүчтөрү төмөндөгөн, ал эми \(S\) ни чыгаруу абдан аз өзгөрүү жараткан.
Калдык диагностикасы айрыкча \(A\) жана \(H\) теңдемелеринде гетероскедастикалуулук, нормалдуу эмес калдыктар жана таасирдүү байкоолор бар экенин көрсөткөн. Ошого карабастан, Cook менен белгиленген байкоолор убактылуу алынып, моделдер кайра курулганда негизги \(M\) жана \(Fr\) жолдорунун багыты жана үстөмдүгү жалпысынан сакталган. Демек, жыйынтыктар негизги процесс түзүлүшүн колдойт, бирок алсыз жолдордун так сандык чоңдуктарын этият чечмелөөнү талап кылат.
Изилдөө колдогон жыйынтыктар
Изилдөө каралган лабораториялык маалымат топтомунда иш жүзүндө болгон сууга кирүү абалы толкун пайда болуусу менен баштапкы шарттарга караганда тыгызыраак байланыштуу экенин колдойт. Айрыкча салыштырмалуу натыйжалуу кирүү массасы \(M\) толкун амплитудасы жана толкун бийиктиги боюнча массага байланышкан эң күчтүү медиатор өзгөрмө болуп чыккан. Кирүү Froude саны \(Fr\) болсо айрыкча толкун бийиктиги үчүн маанилүү экинчи инерция жолу.
Изилдөө ошондой эле жогорку \(Bi\) \(Fr\) жана \(M\) менен терс байланышта экенин; узунураак салыштырмалуу жылышуу аралыгы \(L\) учурдагы эксперименттик шарттарда \(S\) жана \(M\) менен терс жолдорду көрсөткөнүн аныктаган. Авторлордун физикалык түшүндүрмөсү боюнча, узунураак кыймыл жолу вископластикалык массада жайылуу, жукаруу жана ички кайра бөлүштүрүүгө алып келип, натыйжалуу кирүү массасын азайтышы мүмкүн. Бул түшүндүрмө учурдагы материал, геометрия жана параметр диапазону үчүн гана жарактуу.
Изилдөө колдобогон жыйынтыктар
Жол жана медиатордук коэффициенттер катуу мааниде себептүүлүктүн далили эмес. Механизм картасындагы масса-медиатордук жана инерция-медиатордук аймактар чыныгы табияттагы өзгөрбөс режим чек аралары эмес. Байкалган кирүү-абалы моделинин жогорку \(R^2\) мааниси баштапкы шарттардан окуяга чейинки божомолдо ошол эле тактык алынат дегенди билдирбейт.
Carbopol колдонулуп, туруктуу канал геометриясында иштелип чыккан коэффициенттерди табигый тоо тектүү же топурактуу көчкүлөргө түз колдонууга болбойт. Табигый көчкүлөрдө материалдын гетерогендүүлүгү, бөлүнүп-жарылуу, айлана-чөйрөдөн материал кошулушу, тешикче басымы, каныгуу, үч өлчөмдүү жайылуу, татаал топография жана масштаб таасирлери өзүнчө каралышы керек.
Учурдагы алкак эркин беттин майда-чүйдө убакыт эволюциясын, жергиликтүү ылдамдык талааларын, жергиликтүү басымдарды же татаал үч өлчөмдүү гидродинамиканы кайра түзбөйт. Бул жагынан ал жогорку тактыктагы CFD же бөлүкчө негизиндеги чечимдердин физикалык деталын алмаштырбайт.
Изилдөөдө салыштырмалуу максималдуу толкун амплитудасы \(A\) жана максималдуу толкун бийиктиги \(H\) гана максатталган. Толкун периоду, басаңдоо, жээкке көтөрүлүү, сокку басымы жана ашып төгүлүү көлөмү учурдагы моделдин чыгыштары эмес.
Илимий мааниси
Изилдөөнүн негизги салымы “кайсы баштапкы өзгөрмө эң маанилүү?” деген суроону физикалык жактан мазмундуураак суроого айландыруусу: Баштапкы өзгөрүү сууга кирүү абалын кантип өзгөртөт жана бул өзгөргөн кирүү абалы толкун пайда болуусуна кайсы жол аркылуу байланышат?
Бул айырмачылык айрыкча деформациялануучу көчкүлөрдө баалуу. Баштапкы массасы бирдей эки көчкү ар башка реологиялык каршылык, жол боюндагы ар башка жайылуу же ар башка кирүү ылдамдыгынан улам сууга жеткенде бирдей натыйжалуу масса жана инерция алып келбеши мүмкүн. Демек, бирдей баштапкы масса бирдей толкун жообуна кепилдик бербейт.
Маалымат топтомунда пайда болгон эң күчтүү жалпы көрүнүш — масса өткөрүлүшү жана инерция өткөрүлүшү сууга кирүү баскычында толкун жообуна өтүүчү эки негизги жол экени. Салыштырмалуу кирүү калыңдыгы \(S\) маалымат алып жүрөт; бирок \(M\) жана \(Fr\) модельге чогуу киргенде анын өз алдынча салымы кыйла азаят.
Булак жана Ыкма Жөнүндө Эскертүү
Түпнуска аталышы: Process-Chain Analysis of Viscoplastic Landslide-Generated Impulse Waves from Initial Slide Conditions Through Water-Entry Dynamics to Wave Responses
Авторлор: Zhenxia Yuan; Yadong Bian; Song Yin; Xian Liu; Xinchen Pan; Zhenzhu Meng; Bohan Cheng.
Жооптуу автор: Zhenxia Yuan.
Мекемелер: 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.
Журнал: Water.
Басмакана: MDPI.
Библиографиялык жазуу: Water 2026, 18(17), 2189.
DOI: 10.3390/w18172189
Булак түрү: Рецензияланган изилдөө макаласы.
Алынган күнү: 3 июль 2026.
Түзөтүү күнү: 5 август 2026.
Кабыл алынган күнү: 31 август 2026.
Жарыяланган күнү: 3 сентябрь 2026.
Лицензия: Creative Commons Attribution (CC BY).
Каржылоо: Henan Provincial Key Research and Development Program, Grant No. 251111231500; Central Guidance Funds for Science and Technology Local Development Projects, Grant No. 2025ZY01091.
Маалыматтардын жеткиликтүүлүгү: Изилдөөдө колдонулган маалыматтар улантылып жаткан кийинки изилдөөлөрдө пайдаланылып жатканына байланыштуу, маалымат топтому негиздүү өтүнүч боюнча жооптуу автордон алынышы мүмкүн экени билдирилген.
Кызыкчылыктардын кагылышы: Авторлор кызыкчылыктардын кагылышы жок экенин билдиришкен.
Авторлордун салымы: Zhenxia Yuan жана Yadong Bian концептуалдаштыруу жана ыкма; Song Yin, Xian Liu жана Xinchen Pan текшерүү; Zhenxia Yuan жана Bohan Cheng формалдуу талдоо; Zhenzhu Meng маалыматтарды курациялоо; Zhenxia Yuan алгачкы долбоор жана кароо/редакциялоо; Yadong Bian көзөмөл жана долбоорду башкаруу; Zhenxia Yuan жана Yadong Bian каржылоону камсыз кылуу милдеттерин аткарышкан.
Эсептөө чөйрөсү: Талдоолор Python 3.14.6 менен; NumPy 2.5.1, pandas 3.0.5, SciPy 1.18.0, scikit-learn 1.9.0, statsmodels 0.14.6 жана Matplotlib 3.11.1 колдонулуп жүргүзүлгөн. Жол моделдери z-упай стандартташтыруудан кийин ordinary least squares менен курулуп; медиатордук талдоодо random seed 42 менен 1000 bootstrap кайра үлгүлөө колдонулган.
Негизги методологиялык чектөөлөр: Эксперименттер Carbopol, туруктуу канал геометриясы жана чектелген өзгөрмө диапазондору менен лабораториялык масштабда жүргүзүлгөн. Модель коэффициенттерин башка материал, геометрия жана масштабдарда кайра калибрлөө керек. Табигый көчкүлөрдөгү гетерогендүүлүк, бөлүнүп-жарылуу, entrainment, тешикче басымы, каныгуу, үч өлчөмдүү топография жана масштаб таасирлери түз моделденген эмес. Натыйжалуу кирүү массасынын максималдуу толкун бийиктигинин учуруна байланышкан аныктамасы аны катуу мааниде окуяга чейинки көз карандысыз медиатор өзгөрмө болуудан чыгарат. Жол модели сызыктуу жана кошумчалуу; экинчи баскычтагы жогорку VIF маанилери жеке коэффициенттерди чечмелөөнү чектейт. Изилдөө \(A\) жана \(H\) ни гана максат кылат.
Пикир калтырыңыз
E-mail дарегиңиз жарыяланбайт. Милдеттүү талаалар * менен белгиленген