We study the following quasilinear pursuit-evasion model:

$$\{\begin{array}{ll}{u}_t=\mathrm{\Delta }u-\chi \mathrm{\nabla }\cdot (u(u+1{)}^{\alpha }\mathrm{\nabla }w)+u({\lambda }_1-{\mu }_1{u}^{r_1-1}+av),& x\in \mathrm{\Omega },t>0,\\ v_t=\mathrm{\Delta }v+\xi \mathrm{\nabla }\cdot (v(v+1{)}^{\beta }\mathrm{\nabla }z)+v({\lambda }_2-{\mu }_2{v}^{r_2-1}-bu),& x\in \mathrm{\Omega },t>0,\\ 0=\mathrm{\Delta }w-w+v,& x\in \mathrm{\Omega },t>0,\\ 0=\mathrm{\Delta }z-z+u,& x\in \mathrm{\Omega },t>0,\end{array}$$

in a smooth and bounded domain $\mathrm{\Omega }\subset {\mathbb{R}}^n(n\ge 1),$ where $a,b,\chi ,\xi ,{\lambda }_1,{\lambda }_2,{\mu }_1,{\mu }_2>0,$ $\alpha ,\beta \in \mathbb{R},$ and $r_1,r_2>1.$ When $r_1>max\{1,1+\alpha \},r_2>max\{1,1+\beta \},$ it has been proved that if $min\{(r_1-1)(r_2-\beta -1),(r_1-\alpha -1)(r_2-\beta -1)\}>\frac{(n-2{)}_{+}}{n},$ then for some suitable nonnegative initial data $u_0$ and $v_0,$ the system admits a unique globally classical solution which is bounded in $\mathrm{\Omega }\times (0,\mathrm{\infty })$.

The use of unmanned aerial vehicles (UAVs) for photogrammetry allows the rapid acquisition of high-resolution images of geological masses in complex landforms. However, effective analysis of the acquired image information remains a key research issue. At K158 + 837 on the Chongqing-Huaihua Railway, Baima jointed rock masses were reconstructed with high accuracy using UAV close-range photogrammetry technology, and rock discontinuities were extracted from the projected image. The proposed modeling algorithm for jointed rock masses enables the preprocessing of two-dimensional jointed rock mass slopes. Numerical simulations using the disk-based discontinuous deformation analysis method show that the discontinuity network formed by initial cutting significantly affects the subsequent crack development. Meanwhile, simulation results under different scenarios indicate the importance of the pre-reinforcement measures applied to unstable rock masses. The workflow developed based on these results can serve as a reference for the comprehensive acquisition, recognition and numerical modeling analysis of similar jointed rock masses.

In this manuscript, the following chemotaxis system has been considered:

$$\{\begin{array}{ll}{v}_t=\mathrm{\nabla }\cdot (\varphi (v)\mathrm{\nabla }v-\phi (v)\mathrm{\nabla }{w}_1+\psi (v)\mathrm{\nabla }{w}_2)+av-b{v}^{\kappa },& x\in \mathrm{\Omega },t>0,\\ 0=\mathrm{\Delta }{w}_1+\alpha v^{{\gamma }_1}-\beta w_1,& x\in \mathrm{\Omega },t>0,\\ 0=\mathrm{\Delta }{w}_2+\gamma v^{{\gamma }_2}-\delta w_2,& x\in \mathrm{\Omega },t>0,\end{array}$$

where $\mathrm{\Omega }$ is a bounded smooth domain of ${\mathbb{R}}^n(n\ge 1),$ the parameters $a,b,\alpha ,\beta ,\gamma ,\delta ,{\gamma }_1,{\gamma }_2>0,\kappa >1,$ and nonnegative functions $\varphi (\varrho )=(\varrho +1{)}^m,$ $\phi (\varrho )=\chi \varrho (\varrho +1{)}^{\theta -1}$ and $\psi (\varrho )=\xi \varrho (\varrho +1{)}^{l-1}$ for $\varrho \ge 0$ with $m,\theta ,l\in \mathbb{R}$ and $\chi ,\xi >0.$ In the present work, we improve the boundedness criteria established in previous work and further show that under the corresponding critical cases, namely, assume that $\theta +{\gamma }_1=max\{l+{\gamma }_2,\kappa \}\ge m+\frac{2}{n}+1$ with $m>-\frac{2}{n},n\ge 3,$ if one of the following conditions holds:

(a) when $\theta +{\gamma }_1=l+{\gamma }_2=\kappa ,$ if $\theta \ge l\ge 1$ and $\frac{[(\kappa -1-m)n-2](2\alpha \chi -\gamma \xi )}{2(l-1)+(\kappa -1-m)n}=b,$ or $l\ge \theta \ge 1$ and $\frac{2\alpha \chi [(\kappa -1-m)n-2]}{2(\theta -1)+(\kappa -1-m)n}=b;$

(b) when $\theta +{\gamma }_1=\kappa >l+{\gamma }_2,$ if $\theta \ge 1$ and $\frac{2\alpha \chi [(\kappa -1-m)n-2]}{2(\theta -1)+(\kappa -1-m)n}=b,$

then the system still possesses at least a global classical solution, which is bounded in $\mathrm{\Omega }\times (0,\mathrm{\infty })$. Additionally, we have also explored the long time behavior of the classical solution mentioned above.