Tecnología mejorada de absorción de capas totalmente adaptadas basada en depósitos de carbonato
Tian Kun, Li Zhenchun, Huang Jianping, Li Na, Kong Xue, Liu Yujin
(Facultad de Ciencias y Tecnología de la Tierra, Universidad China del Petróleo (Este de China), Qingdao, Shandong 266555)
Proyecto del fondo: Proyecto Nacional 973 (Nº 2011CB202402), Fondo de Innovación de la Universidad del Petróleo (Nº 27R1001046A) y (Y090104)
Acerca del autor: Tian Kun, hombre , es un candidato a doctorado y actualmente participa en la investigación de simulación directa de ondas sísmicas. Correo electrónico: tiankunwudi@163.com.
Resumen: Cuando se utiliza la ecuación de onda para la simulación numérica, dado que el rango de cálculo del medio es definitivamente limitado, esto provocará artificialmente límites de cálculo e inevitablemente producirá reflexiones en los límites. Por lo tanto, es necesario utilizar límites de absorción. condiciones para hacer La energía incidente en el límite se absorbe, lo que reduce el impacto de la reflexión del límite. Se ha demostrado que la tecnología de capa de absorción PML es una tecnología de absorción de límites muy eficaz, que tiene muy buenos efectos en la absorción de ondas corporales y ondas superficiales. Se ha utilizado ampliamente en simulaciones numéricas de ondas elásticas. Sin embargo, la tecnología PML tradicional todavía tiene ciertos problemas en el caso de incidencia de ángulo grande, es decir, en el caso de pastoreo, la atenuación no es suficiente, el coeficiente de reflexión es relativamente grande y después se producirán graves reflejos falsos. La discretización, que reduce el efecto de absorción y afecta la onda real, tiene un impacto grave, que aumenta con el aumento del ángulo de incidencia y el desplazamiento, y tendrá un impacto negativo en otros procesamientos, como las imágenes y la inversión en aplicaciones prácticas. El pastoreo es común, como en áreas de sección delgada, ubicaciones de fuente cercanas al borde del área de estudio, recepción con gran desplazamiento, etc. Basado en la ecuación de velocidad-estrés de primer orden, este artículo propone una tecnología de capa convolucional totalmente coincidente (CPML) y la implementa utilizando el método de diferencias finitas de cuadrícula escalonada. Los resultados muestran que la tecnología CPML mejora efectivamente la absorción en el caso de pastoreo. El efecto, especialmente en el caso de campos de ondas de energía débiles, como los campos de dispersión, es relativamente obvio. Además, no es necesario dividir variables durante el proceso de implementación, lo que hace que la aplicación sea más conveniente y simple, y fácil de implementar mediante programación. La convolución se calcula en forma de recursividad, lo que no aumenta la cantidad de cálculo y la cantidad. La cantidad de almacenamiento no cambia mucho.
Palabras clave: PML; condiciones de contorno de absorción; ondas elásticas; tecnología de capas convolucionales totalmente adaptadas
Una condición de contorno de absorción de capas perfectamente adaptadas para el yacimiento de carbonato karst
p>
Tian Kun, Li Zhenchun, Huang Jianping, Li Na, Kong Xue, Liu Yujin
(Escuela de Geociencias, Universidad del Petróleo de China (China Oriental), Qingdao 266555, China )
p>Resumen: Cuando simulamos a partir de la ecuación de onda, se pueden calcular los límites artificialmente y los reflejos de los límites artificiales inevitablemente, porque el rango de cálculo es limitado, por lo que es necesario absorber condiciones límite para generar la energía que la incidencia sobre el límite se absorbe y el impacto causado por las reflexiones de los límites es mínimo. La condición de límite absorbente de la capa perfectamente adaptada (PML) ha demostrado ser un método muy eficiente para el modelado numérico directo de la ecuación de onda elástica para absorber tanto las ondas corporales como las ondas superficiales. Pero en la incidencia rasante, el método PML discreto clásico todavía presenta algunos problemas, es decir, la atenuación no es suficiente, el coeficiente de reflexión es grande hasta cierto punto, esto puede producir una reflexión espuria grave después de la dispersión y reducir el efecto de absorción, lo que luego causa graves consecuencias. influencias en el campo real y estas influencias aumentan con el aumento del ángulo de incidencia y el desplazamiento, además, pueden afectar negativamente las imágenes y la inversión y otros tratamientos en la aplicación práctica. Pero la incidencia rasante es natural, por ejemplo.
ance en el caso de cortes de malla muy delgados, así como fuentes ubicadas cerca del borde de la malla, o en el caso de receptores ubicados en un desplazamiento muy grande. En este artículo, presentamos una condición de límite de absorción de PML mejorada en incidencia de pastoreo. que se denomina PML de convolución (CPML) y se basa en la ecuación de onda diferencial parcial de velocidad-estrés de primer orden y se implementa mediante un método de diferencias finitas de cuadrículas escalonadas. Los resultados muestran que la técnica CPML puede obtener un mejor resultado en el caso de. incidencia rasante, especialmente para campos de dispersión y similares como este campo de ondas de reflexión débil. Y su implementación es conveniente y fácil de usar y programar debido a la variable no dividida. El costo del cálculo no aumenta debido a la recursividad cuando la convolución. se calcula. El almacenamiento de memoria también es similar al del PML clásico
Palabras clave: PML; condición de contorno absorbente; modelado directo;
Debido al rápido desarrollo de la potencia informática, la investigación sobre métodos de simulación numérica de la propagación de ondas sísmicas en medios complejos se ha vuelto más extensa y profunda en las últimas décadas. Entre ellos, el más utilizado es el método de diferencias finitas [1 ~ 3]. En la simulación directa de diferencias finitas de un área limitada, ya sea resolviendo las ecuaciones de ondas electromagnéticas de Maxwell o las ecuaciones de dinámica de ondas elásticas, el método completo propuesto por Bérenger ( 1994) La condición límite de absorción de capa perfectamente adaptada (PML) es la condición límite de absorción más efectiva [4]. La ecuación de onda en la capa completamente emparejada puede considerarse como una extensión de la ecuación de onda convencional. Cuando la onda se propaga en ella, la fase cambia y la amplitud se atenúa exponencialmente los parámetros elásticos del PML y del área principal de estudio; pero los coeficientes de atenuación son diferentes y la impedancia de la onda coincide completamente. La teoría No se producirá ninguna reflexión.
Chew y Weedon (1994) introdujeron un sistema complejo de coordenadas de extensión para formular las condiciones límite de absorción de PML [5]; Rappapport (1995) demostró que el medio PML es equivalente a introducir medios anisotrópicos en la región límite de absorción [6]; de estudios [5, 7] muestran que la condición límite de absorción de PML tiene mejores propiedades que la condición límite de absorción de decaimiento exponencial [8, 9], la condición límite de absorción de Liao [10], la condición límite de absorción de Higdon [11] y el límite de absorción de aproximación paraxial condición [12]. Excelentes propiedades de absorción. Chen et al. (1997) y Wang y Oristaglio (2000) aplicaron con éxito las condiciones de contorno de absorción de PML a la solución de ecuaciones de ondas electromagnéticas [7, 13]. En los últimos años, las condiciones límite de absorción de PML también se han aplicado al modelado directo de diferencias finitas de ondas acústicas y ondas elásticas [13 ~ 16]. Además, Teixeir y Chew (1999) y Chew y Liu (1996) implementaron condiciones de contorno de absorción de PML en el sistema de coordenadas cilíndrico y el sistema de coordenadas esférico [17, 18].
Sin embargo, el PML tradicional también tiene ciertos defectos. El coeficiente de reflexión de la capa de PML discreta no es estrictamente 0, especialmente en ángulos de incidencia grandes. En este caso, una parte considerable de la energía se devuelve. El área principal de estudio en forma de ondas reflejadas, y son comunes situaciones con grandes ángulos de incidencia, como áreas de corte fino, ubicaciones de fuentes cercanas al borde del área de estudio, recepción con gran desplazamiento, etc. Para superar este problema, se propuso y aplicó a la solución de las ecuaciones de Maxwell (Roden y Gedney, 2000) un método PML que modifica la transformación de coordenadas complejas para mejorar el efecto de absorción en grandes ángulos incidentes después de la discretización y se denominó convolución Capa de coincidencia completa. (CPML) [19]. Este artículo utiliza este método para resolver la ecuación de onda elástica. Los resultados muestran que la tecnología CPML mejora efectivamente el efecto de absorción del pastoreo y no necesita dividir variables durante el proceso de implementación. La aplicación es más conveniente y simple. en forma de recursividad, que no cambia mucho a medida que aumenta la cantidad de cálculo.
1 Principios básicos de CPML
La tecnología PML esencialmente transforma la ecuación de onda en coordenadas complejas dentro de la capa PML. Para las coordenadas transformadas, la forma de la ecuación y su solución permanecen sin cambios. , pero la solución de coordenadas original está atenuada. La transformación de coordenadas en el dominio de frecuencia en PML tradicional es (tomando x como ejemplo):
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
En CPML, ecuación (3) Se ha ampliado para hacer su forma más general:
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
Donde αx≥0, κx≥ 1, se puede ver que la PML tradicional es un caso especial de CPML en el caso de αx=0, κx=1.
Transforme la relación de transformación de coordenadas al dominio del tiempo y utilice la transformada inversa de Fourier que representa 1/sx para obtener
Actas de la Exploración y Desarrollo Internacional de Petróleo y Gas No Convencionales (Qingdao) Conferencia
De acuerdo con la ecuación (4), podemos obtener
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
δ y 1 y e-atH(t) y 1/(a+iω) es un par de transformada de Fourier, por lo que
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
De esta manera, la relación de transformación de coordenadas en CPML se divide en dos partes. La primera parte es fácil de calcular. Simplemente calcule la derivada parcial espacial dividida por el coeficiente κ. Lo más importante es calcular la convolución de la segunda parte. .
Bajo la cuadrícula escalonada discreta, la convolución en el tiempo n se escribe como
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
En el caso de la cuadrícula escalonada, se puede escribir como
p>
Actas de la Conferencia Internacional de Exploración y Desarrollo de Gas y Petróleo No Convencional (Qingdao)
Entre ellos
Actas de la Conferencia Internacional de Exploración y Desarrollo de Petróleo y Gas No Convencionales (Qingdao)
En la fórmula:
Actas de la Conferencia Internacional de Exploración y Desarrollo de Petróleo y Gas No Convencionales (Qingdao)
Debido a que la fórmula (12) es una forma exponencial simple, la fórmula (11) se puede escribir como una fórmula La forma de empuje:
Actas de la Exploración y Desarrollo Internacional de Petróleo y Gas No Convencionales ( Qingdao) Conferencia
De esta manera, se realizan cálculos normales en el área de investigación principal y la transformación de coordenadas se realiza en la capa PML a través de la siguiente relación El cálculo final es suficiente:
Actas de la Conferencia Internacional de Exploración y Desarrollo de Petróleo y Gas No Convencionales (Qingdao)
Donde Ψx se puede obtener mediante la recursión de la ecuación (13).
Este es el principio básico de CPML y se puede extender fácilmente a las direcciones y y z. La situación de las esquinas es similar a la PML tradicional y se pueden considerar todas las direcciones simultáneamente.
2 Cálculo y comparación de simulación directa
Para verificar la efectividad y superioridad de CPML, a continuación se realiza un cálculo de prueba para medios isotrópicos y se compara con el cálculo de PML tradicional. Se compararon los resultados. La Figura 1 es un modelo medio en capas con un punto de muestreo de 301 × 301 y un espaciado de cuadrícula de 5 m × 5 m. Las velocidades de onda longitudinal y transversal de la capa superior son 3000 m/s y 1600 m/s respectivamente, y las ondas longitudinales y transversales de la capa inferior. las velocidades son 3300 m/s y 1900 m/s respectivamente, la densidad es 2800 kg/m3, las coordenadas de origen son (150, 5), el espacio entre canales es de 5 m, la velocidad de muestreo es de 0,5 ms y el tiempo de muestreo es de 1,5 s. 2 y 3 son la comparación de registros de pistola del componente z y instantáneas de campo de onda con la frecuencia principal de 20 Hz respectivamente, se puede ver claramente que PML produce una reflexión relativamente fuerte en la interfaz superior y tiene un impacto en el registro de pistola. mientras que el efecto de CPML es mejor.
La Figura 4 y la Figura 5 son, respectivamente, una comparación de registros de pistola de componente z e instantáneas de campo de onda con una frecuencia principal de 30 Hz. También se puede ver que CPML tiene un mejor efecto de absorción que PML. La Figura 6 es un diagrama de comparación de formas de onda de tres canales individuales. Las coordenadas de los canales individuales 1, 2 y 3 son (201, 5), (251, 5), (301, 5) respectivamente. Comparación de formas de onda de 0 a 1,2 s, el lado derecho es la comparación de formas de onda de 0,3 a 1,2 s. También se puede ver que CPML es mejor que PML y, a medida que aumenta el desplazamiento y el ángulo de incidencia, más obvio es el contraste. cuanto mayor sea el efecto de reflexión en la PML tradicional.
Figura 1 Modelo medio en capas
Figura 2 Comparación de registros de armas con una frecuencia principal de 20 Hz
Figura 3 Una instantánea del campo de ondas de 0,3 s con una frecuencia principal de 20 Hz (Comparación de superior) y ampliación parcial (inferior)
Figura 4 Comparación de registros de armas con una frecuencia principal de 30 Hz
La Figura 7 es un modelo mediano que contiene dispersores, en ( 225, 15) Hay un dispersor. Las velocidades de onda longitudinal y transversal del dispersor son 2700 m/s y 1300 m/s respectivamente. Los demás parámetros son los mismos que antes. Las Figuras 8 y 9 son, respectivamente, una comparación de registros de pistola de componente z e instantáneas de campo de onda con una frecuencia principal de 30 Hz. Se puede ver que la falsa reflexión de PML tiene un impacto considerable en el campo de dispersión, lo que afectará la imagen, la inversión y. otras aplicaciones prácticas.
Figura 5 Comparación de una instantánea del campo de onda de 0,3 s (arriba) y una ampliación parcial (abajo) con la frecuencia principal de 30 Hz
3 Conclusión
Este artículo se basa sobre no división El método de convolución generalizó la tecnología PML tradicional en la situación de diferencias finitas de cuadrícula escalonada, dedujo la fórmula recursiva de la relación de transformación de coordenadas en el dominio del tiempo de CPML, realizó un análisis comparativo a través de cálculos de prueba de modelos y obtuvo la siguiente conclusión: (1) El coeficiente de reflexión del PML tradicional no es 0 después de la discretización, especialmente en ángulos de incidencia grandes y grandes compensaciones, se producirán fuertes reflexiones falsas, lo que reducirá el efecto de absorción y afectará seriamente el campo de onda real, lo que aumentará. con el aumento del ángulo de incidencia y el desplazamiento (2) La tecnología CPML mejorada puede mejorar la atenuación del campo de onda en situaciones de pastoreo, reducir el coeficiente de reflexión y mejorar eficazmente el efecto de absorción de la capa coincidente (3) Para modelos medios complejos que contienen dispersores; , El PML tradicional también producirá reflejos falsos, que afectarán más seriamente los campos de ondas de energía débiles, como los campos de dispersión, y afectarán negativamente la migración, la inversión y otros procesos en aplicaciones prácticas. La tecnología CPML también puede mejorar el efecto de absorción y reducir los efectos adversos. efectos (4) la tecnología CPML utiliza el método de convolución sin dividir variables para el cálculo y calcula la convolución mediante recursividad, de modo que su capacidad de almacenamiento y cantidad de cálculo son comparables al PML tradicional. No hay muchos cambios en comparación con ambos, y el costo computacional; no se incrementará.
La tecnología CPML en este artículo puede mejorar efectivamente el efecto de absorción del pastoreo sin aumentar el costo computacional. Esto es muy bueno para muchas situaciones prácticas, como áreas de sección delgada y fuentes de terremotos cercanas al borde del estudio. Área y grandes perspectivas de aplicación de compensación. Esta también es una mejor opción para las condiciones límite de absorción en procesamientos como las imágenes de migración. Sin embargo, CPML no cambia la idea básica del PML tradicional, por lo que todavía tiene algunos problemas, como la inestabilidad inherente de algunos medios anisotrópicos y el coeficiente de reflexión de la capa coincidente después de la discretización no es estrictamente cero. Este artículo utiliza las condiciones de contorno de Dirichlet alrededor de la capa coincidente, en investigaciones posteriores, se pueden reemplazar por condiciones de contorno de absorción paraxial para mejorar aún más el efecto de absorción después de la discretización. Además, las condiciones de contorno de absorción CPML también para la ecuación de onda de desplazamiento de segundo orden. necesitan ser estudiados más a fondo. Una mayor investigación sobre la tecnología CPML y los métodos relacionados posteriores será beneficiosa para la simulación directa bajo el complejo modelo de velocidad cerca de la superficie en el área de exploración de carbonatos occidental, especialmente la investigación sobre la comprensión del mecanismo de las ondas dispersas cerca de la superficie en la exploración de carbonatos. área, que puede servir para la futura exploración y desarrollo de áreas de exploración de carbonatos.
Figura 6 Comparación de formas de onda CPML y PML de diferentes canales individuales
Figura 7 Modelo mediano con dispersores
Figura 8 Modelo de pistola de componente Z con dispersor Comparación de registros
Figura 9 Comparación de la instantánea del campo de onda de 0,3 s del componente z (arriba) y la ampliación parcial (abajo) del modelo que contiene dispersores
Referencias
[1]Alterman Z y Karal F C. Propagación de ondas elásticas en medios estratificados mediante métodos de diferencias finitas. Boletín de la Sociedad Sismológica de América, 1968, 58: 367-398.
[2]Madariaga, R..Dinámica de una falla circular en expansión. Boletín de la Sociedad Sismológica de América, 1976, 65: 163-182.
[3]Virieux, J.P-SV propagación de ondas en medios heterogéneos: método de diferencia finita de tensión-velocidad . Geophysics, 1986, 51: 889-901.
[4] Bérenger, J.P. Una capa perfectamente adaptada para la absorción de ondas electromagnéticas. Journal of Computational Physics, 1994, 114: 185 -200. p>
[5]Chew W.C. y W.H.Weedon.Un medio tridimensional perfectamente adaptado a partir de ecuaciones de Maxwell modificadas con coordenadas estiradas.Microwave and Optical Technology Letters, 1994, 7: 599-604.
[6]Rappaport, C.M.. Condiciones de contorno absorbentes perfectamente adaptadas basadas en mapeo anisotrópico con pérdida del espacio IEEE Microwave Guided Wave Lett., 1995, 5: 90-92.
[7 】Chen, Y.H., Chew, W.C., Oristaglio, M.L.. Aplicación de capas perfectamente adaptadas al modelado transitorio de problemas EM del subsuelo, 1997, 62: 1730-1736.
[8]Marfurt, K.J. Modelado en diferencias y elementos finitos de las ecuaciones de ondas escalares y elásticas.Geophysics, 1984, 49: 533-549.
[9]Shin, C..Spongeboundar
y condición para el modelado en el dominio de la frecuencia. Geophysics, 1995, 60: 1870-1874.
[10]Liao, Z.P., Wong, H.L., Yang, B.P., et al.. Un límite de transmisión para el análisis de ondas transitorias. .Scientia Sinica, 1984, 27 (10): 1063-1076.
[11] Higdon, R.L.Condición de frontera absorbente para ondas elásticas, 1991, 56: 231-241.
[12]Engquist, B., Majda, A..Condiciones de contorno absorbentes para la simulación numérica de ondas. Matemáticas de Computación, 1977, 31(9): 629-651.
[13]. Wang, T., Oristaglio, M.L.Simulación 3D de estudios GPR sobre tuberías en suelos dispersivos, 2000, 65: 1560-1568.
[14]Liu, Q.H., Tao, J.P..La combinación perfecta. capa para ondas acústicas en medios absorbentes.Acoust.Soc.Am., 1997, 102: 2072-2082.
[15]Zeng, Y., He, J., Liu, Q.H..La aplicación de la capa perfectamente adaptada en el modelado numérico de la propagación de ondas en medios poroelásticos Geophysics, 2001, 66: 1258-1266.
[16]Imhof, M.G..Calculando el efecto sísmico de estructuras subterráneas 3D y topografía con el. método de diferencias finitas.72nd Internat.Mtg., Soc.Expl.Geophys., Expanded Abstracts, 2002: 1939-1942.
[17]Teixeira, F.L., Chew, W.C..Sobre causalidad y estabilidad dinámica de capas perfectamente combinadas para simulaciones FDTD IEEE Transactions on Microwave Theory and Techniques, 1999, 47: 775-785.
[18]Chew, W.C., Liu, Q.H..Capas perfectamente combinadas para elastodinámica: una nueva. condición de contorno absorbente. Journal of Computationa
l Acoustics, 1996, 4: 341-359.
[19]Roden, J.A. y S.D.Gedney.Convolution PML (CPML): una implementación FDTD eficiente de CFS-PML para medios arbitrarios. Microondas y Cartas de tecnología óptica, 2000, 27: 334-339.