Vegetación superficial
En el modelo geológico tridimensional existen relaciones cruzadas entre estratos y entre estratos y fallas. Existen varias relaciones de intersección entre estratos (como pellizco, intrusión, lente, etc.). Hay cortes de fallas normales y cortes de fallas inversas entre estratos y fallas. Por lo tanto, después de construir todas las interfaces estratigráficas y las interfaces de falla, es necesario realizar un procesamiento de corte y ajuste entre estratos y estratos y entre estratos y niveles de falla.
El problema de corte entre capas implica primero el cálculo de la intersección entre TIN y TIN, y finalmente se reduce al juicio entre unidades triangulares. Lindenbeck et al. (2002) desarrollaron un algoritmo TRICUT basado en árboles OBB (Oriented Bounding Box) para lograr cortes entre los planos de falla y el plano de tierra. Möller (1997) propuso un algoritmo de prueba de intersección de triángulos rápido, cuyo núcleo es el "método de superposición de intervalos", que se considera un algoritmo más rápido que el árbol OBB.
Hay dos triángulos espaciales T1 y T2, que están compuestos por puntos , , y puntos , , , respectivamente. Los planos espaciales donde se ubican los dos triángulos se registran como π1 y π2 respectivamente. Primero, es necesario aclarar la relación espacial entre las unidades triangulares: los planos π1 y π2 pueden ser paralelos a ****, pueden ser paralelos a diferentes planos o pueden no ser paralelos. Al mismo tiempo, T1 y T2 pueden ser paralelos. se cruzan o no (Fig. 5.6).
Figura 5.6 Configuración de intersección triángulo-triángulo
a π1||π2, π1≠π2; b. T2 Disjoint
Zhang Fang (2005) mejoró el método de Möller y consideró las dos situaciones en las que π1 y π2 tienen intersecciones (Figura 5.6c, d) sin considerar la superficie *** Hay intersecciones (Figura 5.6b) , que es principalmente para mejorar el algoritmo.
(1) Primero, escribe la ecuación normal puntual del plano π2 (o π1): donde N2= ×; d2=-N2-.
(2) Calcular la distancia direccional (distancia con signo) desde los tres puntos del triángulo T1 (o T2) al plano π2 (o π1) (Figura 5.7):
Tres geología -dimensional Investigación sobre la aplicación del modelado y la ingeniería digital del espacio subterráneo
Según la distancia direccional positiva, negativa o cero, se puede juzgar si los triángulos son paralelos, **** o pueden cruzarse .
Si d≠0 (i=0, 1, 2; es decir, ningún punto está en el plano π2), y los signos de d son iguales, entonces T1 está completamente a un lado del plano π2 (Figura 5.6a), es decir, es imposible intersecar.
Si d = 0 (i = 0, 1, 2), significa que T1 está ubicado en un lado del plano π2*** (Figura 5.6b), entonces la prueba de bidimensional se ingresa la intersección del triángulo; de lo contrario, T1 y T2 pueden cruzarse, y hay dos casos de intersección. En la Figura 5.6c, es una intersección real y en la Figura 5.6d, es una no intersección.
(3) Unifique las dos posibles situaciones de intersección, primero genere la línea de intersección P1P2 del triángulo T1 y el plano π2, y luego genere la línea de intersección Q1Q2 del triángulo T2 y el plano π1 (Figura 5.6 cd ). 5.6c, d), obviamente ambos están ubicados en la intersección L y son parte de la intersección L y la ecuación paramétrica de L es L=O tD, donde D=N1×N2, N1 y N2 son π1 y π2 respectivamente. línea normal exterior.
Figura 5.7 Diagrama de distancia direccional
Tome P1P2 como ejemplo (Figura 5.6c) para ilustrar el proceso de cálculo. Como se muestra en la Figura 5.7, sin pérdida de generalidad, sea y ubicado en un lado del plano π2, y 11 ubicado en el otro lado del plano π2.
Proyecte y respectivamente en la intersección del plano L:
Digitalización del espacio subterráneo e ingeniería, modelado geológico tridimensional e investigación de aplicaciones
Al calcular el punto de intersección B de 10V11 y L, solo necesita Para encontrar un parámetro t1 arriba, B= ∩L=O t1D.
(4) Determine la relación entre P1P2 y Q1Q2. Si hay superposición entre los dos, extraiga los segmentos de línea superpuestos. De lo contrario, significa que los dos triángulos están separados.
Después de calcular la línea de intersección, la superficie TIN se reconstruye con restricciones.