In chaotic systems the typical situation is that at a low level trajectories of points are wild, but overall there is a nice statistical description of the system.
For example, consider the trajectories of a system of differential equations. Suppose there is some sort of chaos like in the Lorenz system of differential equations. Each trajectory is very wild and sensitive to the initial condition. However, on the global scale there is a strange attractor that attracts all initial conditions.
It seems that this is the typical situation, that there are finitely many attractors that attract all initial conditions, and on each one of them there is a measure which describes the frequency of visits of each trajectory to some region of the attractor. That is, one can say on average how many times a point visits a certain region near the attractor. Moreover, this structure is stable.
I would like to understand why chaotic, irregular behavior of trajectories lead to a sort of regular behavior when one considers all trajectories together and takes a global point of view?
In order to understand the situation a bit better, I like to understand what is in the complement. That is, are there chaotic systems where there is no nice, global picture as described above? How wild should the trajectories of a system (say in $\mathbb{R}^3$) be so that the global picture is also ``chaotic''? Is this even possible?
What does it even mean for the global picture to be chaotic? Infinitely many strange attractors? Or the other extreme: lack of a strange attractor and any specific pattern?
What is the most chaotic system one can have?
There's a lot here, so I'll try to address "why chaotic, irregular behavior of trajectories leads to a sort of regular behavior".
The most classical version of this theory is that of $C^2$ diffeomorphisms $f$ with compact hyperbolic attractors ($C^2$ is needed so that the derivative $df$ is bounded in operator norm). In this situation, the diffeomorphism may behave however 'orderly' or 'chaotically' you can imagine away from the attractor, but in a neighborhood of the attractor, you can apply a very nice theory (that of SRB measures) which indicates that the statistics of orbits (of a Lebesgue-full set of initial conditions) in that neighborhood are governed by a nice measure, which is supported on the attractor itself. Precisely, there is a probability measure $\mu$ such that, given a continuous observable $f$ and a generic point $x \in U$ where $U \supset \Lambda$ is attracted to the attractor $\Lambda$, we have $$ \frac{1}{n} \sum_{k = 0}^{n-1} f^k x \rightarrow \int_{\Lambda} f d\mu $$ Moreover, the measure $\mu$ is in the best possible sense "compatible" with the Lebesgue measure.
The reason "why" is interesting as well: it turns out that in the above situation, there is a way to symbolically encode the dynamics on $\Lambda$ as a subshift of finite type (aka a topological markov shift). This encoding gives rise to a special measure $\mu$ on $\Lambda$ with all the nice properties I've mentioned.
I should mention that hyperbolicity is important: it means more-or-less that the asymptotic dynamics is governed by the derivative term at a point (technically, it means that the eigenvalues of $df$ are bounded away from the unit circle).
This stuff is the province of differentiable ergodic theory. For reference, you should google the authors Ruelle, Bowen, and Lanford.