MX2013000646A - Sistema y metodos para crear superficie en espacio con falla. - Google Patents

Abstract

La presente invención se refiere a sistemas y métodos para crear una superficie en un espacio con falla, que incluye el uso de técnicas de interpolación.

Description

SISTEMAS Y METODOS PARA CREAR SUPERFICIE EN ESPACIO CON FALLA Campo de la Invención La presente invención generalmente se refiere a sistemas y métodos para crear una superficie en un espacio con falla. Más particularmente, la presente invención se refiere a la creación de una superficie por interpolación en un espacio con falla.

Antecedentes de la Invención Las técnicas de interpolación son frecuentemente usadas en la industria de petróleo y gas para crear una superficie, algunas veces referida como un horizonte, el cual se puede usar para ubicar hidrocarburos en una formación subsuperficial . La creación de una superficie desde los topes del pozo, los cuales son puntos en un pozo que representan un nivel de superficie en un horizonte en este punto, en un espacio con falla siempre es desafiante. También es posible usar información secundaria, tal como, pero no limitado a, datos sísmicos, para ayudar en la creación de la superficie interpolada. De manera ideal, la interpolación se debe hacer en el espacio sin falla y luego trasladarse de nuevo en el campo de falla. Esta técnica, la cual es comúnmente referida como eliminación de falla global, efectivamente restaura el espacio con falla al espacio sin falla en una base global antes de que ocurra la falla. Sin embargo, la eliminación de Ref. 238403 falla global es muy difícil, especialmente en la presencia de muchas fallas .

La rasterización de refinamiento es una técnica de interpolación que traduce una aproximación gruesa de datos conocidos en una aproximación más fina incrementando el número de nodos interpolados en una base etapa por etapa. La interpolación de rasterización de refinamiento por lo tanto, proporciona un procedimiento para crear una superficie en el espacio con falla directamente. Genera superficies muy lisas cuando no hay falla, y el desempeño es muy bueno, pero frecuentemente crea artefactos indeseables significativos 102 a lo largo de la falla y en sus bordes como se ilustra por la superficie interpolada 100 en la FIG. 1 El krigeado tradicional, el cual es sinónimo de la predicción óptima, es otra técnica de interpolación, que predice valores desconocidos de datos observados en ubicaciones conocidas. El krigeado usa variogramas para expresar la variación espacial, y minimiza el error de valores predichos, los cuales se estiman por la distribución espacial de los valores predichos . El krigeado tradicional en la presencia de una falla puede generar bordes más limpios a lo largo de la falla cuando los entornos de búsqueda son restringidos a un lado de la falla 202 u otro lado de la falla 204 como se ilustra por la superficie interpolada 200 en la FIG. 2. El krigeado tradicional, sin embargo, no siempre elimina los artefactos 302 cerca de las terminaciones de la falla como se ilustra por la superficie interpolada 300 en la FIG. 3.

Por lo tanto, existe una necesidad de crear superficies en un espacio con falla que reduzca la producción incidental de artefactos a lo largo de las fallas y cerca de las terminaciones de falla.

Breve Descripción de la Invención La presente invención cumple las necesidades anteriores y supera una o más deficiencias del arte previo proporcionando sistemas y métodos para crear una superficie en un espacio con falla usando técnicas de interpolación.

En una modalidad, la presente invención incluye un método para crear una superficie que tiene una o más fallas, el cual comprende: i) estimar un nivel de superficie en cada lado de cada falla, una diferencia de elevación ente los niveles de superficie estimados en cada lado de una falla respectiva que representa un valor de dislocación de falla respectivo; ii) ajustar cada valor de dislocación de falla usando un procesador de computadora; iii) ubicar cada tope de pozo local dentro de un entorno de búsqueda predefinido de un nodo de rejilla dentro de una rejilla en la superficie; iv) identificar cada falla que cruza un vector entre el nodo de rejilla y cada tope de pozo local respectivo, cada vector entre el nodo de rejilla y cada tope de pozo local respectivo que representa un vector de intersección; v) acumular cada valor de- dislocación de falla ajustado a lo largo de cada vector de intersección respectivo, él cual representa un valor de dislocación de falla acumulado para cada vector de intersección respectivo; vi) ajustar un valor para cada tope de pozo local con base en un valor de dislocación de falla acumulado respectivo; vii) realizar la interpolación en el nodo de rejilla usando el valor ajustado para cada tope de pozo local; y viii) repetir las etapas c) - g) para cada nodo de rejilla dentro de la rejilla.

En otra modalidad, la presente invención incluye un dispositivo portador de programa no transitorio que porta tangiblemente instrucciones ejecutables por computadora para crear una superficie que tiene una o más fallas. Las instrucciones son ejecutables para implementar: i) estimación de un nivel de superficie en cada lado de cada falla, una diferencia de elevación entre los niveles de superficie estimados en cada lado de una falla respectiva que representa un valor de dislocación de falla respectivo; ii) ajuste de cada valor de dislocación de falla usando un procesador de computadora; iii) ubicación de cada tope de pozo local dentro de un entorno de búsqueda predefinido desde un nodo de rejilla dentro de una rejilla en la superficie; iv) identificación de cada falla que cruza un vector entre el nodo de rejilla y cada tope de pozo local respectivo, cada vector entre el nodo de rejilla y cada tope de pozo local respectivo que representa un vector de intersección; v) acumulación de cada valor de dislocación de falla ajustado a lo largo de cada vector de intersección respectivo, el cual representa un valor de dislocación de falla acumulado para cada vector de intersección respectivo; vi) ajuste de un valor para cada tope de pozo local con base en un valor de dislocación de falla acumulado respectivo; vii) realización de la interpolación en el nodo de rejilla usando él valor ajustado para cada tope de pozo local; y viii) repetición de las etapas c) - g) para cada nodo de rejilla dentro de la rejilla.

Los aspectos, ventajas y modalidades adicionales de la invención llegarán a ser evidentes por aquellos expertos en el arte a partir de la siguiente descripción de las diversas modalidades y figuras relacionadas .

Breve Descripción de las Figuras La presente invención se describe a continuación con referencias a las figuras acompañantes en las cuales los elementos similares son referenciados con números de referencia similares, y en las cuales: La FIG. 1 ilustra una superficie y artefactos producidos por la interpolación de rasterización de refinamiento convencional usando datos reales.

La FIG. 2 ilustra una superficie producida por interpolación de krigeado tradicional y los mismos datos usados para producir la superficie interpolada en la FIG. 1.

La FIG. 3 ilustra una superficie y artefactos producidos por la interpolación de krigeado tradicional y datos artificiales.

La FIG. 4 es un diagrama de flujo que ilustra una modalidad de un método para implementar la presente invención.

La FIG. 5 ilustra una superficie creada por interpolación de acuerdo con la presente invención usando datos reales.

La FIG. 6 ilustra una superficie creada por interpolación de acuerdo cón la presente invención usando los mismos datos usados para producir la superficie interpolada en la FIG. 3.

La FIG. 7 es un diagrama de bloque que ilustra una modalidad de un sistema de computadora para implementar la presente invención.

Descripción Detallada de la Invención El tema de la presente invención se describe con especificidad, sin embargo, la descripción del mismo no se propone que limite el alcance de la invención. El tema por consiguiente, también puede ser incluido en otras maneras, para incluir diferentes etapas o combinaciones de etapas similares a las descritas en la presente, en conjunto con otras tecnologías presentes o futuras. Además, aunque el término "etapa" se puede usar en la presente para describir diferentes elementos de métodos empleados, el término no se debe interpretar como que implica cualquier orden particular entre o en medio de varias etapas en la presente descritas a menos que se limite expresamente de otra manera por la descripción a un orden particular. Mientras que la siguiente descripción se refiere a la industria del petróleo y gas, los sistemas y métodos de la presente invención no se limitan a esto y también se pueden aplicar a otras industrias para lograr resultados similares.

Descripción del Método Con referencia ahora a la FIG. 4, un diagrama de flujo ilustra una modalidad de un método 400 para implementar la presente invención. El método 400 crea una superficie en un espacio con falla usando técnicas de interpolación. El método 400 en este ejemplo usa solamente topes de pozo a través de fallas para crear la superficie interpolada. El método 400 no solamente genera bordes agudos alrededor de las fallas, sino también interpola suavemente lejos de ellas.

En la etapa 402, los datos de entrada se almacenan usando la interfaz de cliente y/o interfaz de vídeo descritas con referencia a la FIG. 7. Los datos de entrada pueden comprender topes de pozo y fallas que comprenden numerosos segmentos de falla. Cada tope de pozo preferiblemente se almacena en árbol KD, el .cual es una estructura de datos bien conocida para la recuperación rápida y eficiente de los datos de entrada. Cada segmento de falla preferiblemente se almacena en un árbol Quad, el cual también es una estructura de datos bien conocida usada para la recuperación rápida y eficiente.

En la etapa 404, se estima el nivel de superficie a lo largo de ambos lados de cada falla. La estimación preferiblemente se hace encontrando cada tope de pozo, que se almacena como datos de entrada en la etapa 402, en cada bloque de falla e interpolando el nivel de superficie a lo largo de ambos lados de cada falla usando uno o. más, pero no necesariamente todos, los topes de pozo encontrados en cada bloque de falla respectivo. Se piensa que la interpolación abarca cualquier técnica de interpolación conocida tal como, por ejemplo, krigeado, cokrigeado colocado, algoritmos de simulación y cosimulación colocada. Cada dislocación de falla, la cual se define como la cantidad de desplazamiento vertical (es decir, cambio de elevación) entre los niveles de superficie estimados, se almacena como un valor en la estructura de datos de árbol Quad descrita con referencia a la etapa 402 con los niveles de superficie estimados. Como se ilustra en la FIG. 5, la cual representa una superficie interpolada 500 usando datos de entrada reales en la etapa 402, un bloque de falla se define por segmentos de falla 502, 504, 506 y un borde 508 de la superficie interpolada 500. Otro bloque de falla se define por los segmentos de falla 510, 512, 514 y 516. Los datos reales usados para la superficie interpolada 500 incluyen 9,995 segmentos de falla y 329 topes de pozo, los cuales se ilustran dentro de una rejilla 540 de 500 por 500. Una falla 501 incluye el segmento de falla 502 en un extremo de la falla 501 y otro segmento dé falla 514 hacia otro extremo de la falla 501. Cada bloque de falla puede contener uno o más topes de pozo como se demuestra por el bloque de falla que contiene topes de pozo 518-526. El bloque de falla definido por un borde 528 para la superficie interpolada 500 y segmentos de falla 530, 532 y 534 también puede contener uno o más topes de pozo. El nivel de superficie 538 a lo largo de un lado de la falla 501 por lo tanto, se estimó i) hallando cada tope de pozo en el bloque de falla definido por segmentos de falla 502, 504, 506 y el borde 508 de la superficie interpolada 500; y ii) krigeando este nivel de superficie usando uno o más topes de pozo; en el bloque de falla. El nivel de superficie 536 a lo largo de otro lado de la falla 501 por lo tanto, se estimó i) hallando cada tope de pozo en el bloque de falla definido por los segmentos de falla 530, 532, 534 y el borde 528 de la superficie interpolada 500; y ii) krigeando este nivel de superficie usando uno o más topes de pozo en el bloque de falla. El nivel de superficie en cada lado de la falla 501 se puede estimar de esta manera (bloque de falla por bloque de falla) desde un extremo de la falla 501 al otro extremo de la falla 501. Alternativamente, esta etapa puede comenzar en cualquier parte entre cada extremo de la falla 501 y proceder de cualquier manera o secuencia hasta que se estime el nivel de superficie completo en cada lado de la falla 501.

En la etapa 406, cada valor de dislocación de falla se ajusta i) alisando los niveles de superficie estimados en la etapa 404 a lo largo de cada lado de cada falla que pertenece al mismo bloque de falla; y ii) ahusando cada valor de dislocación de falla a cero en cada extremo distal de cada falla. El ahusamiento no debe efectuar una falla que interseca otra falla y por consiguiente, necesariamente ocurre entre cada extremo distal de una falla y el punto más cercano en el cual otra falla interseca la falla.

En la etapa 408, los topes de pozo locales dentro de un entorno de búsqueda predefinido de un nodo de rejilla se ubican usando los topes de pozo almacenados como datos de entrada en la etapa 402, sin considerar si los topes de pozo locales cruzan los segmentos de falla y/o bloques de falla. Cada tope de pozo local tiene un valor de tope de pozo local que representa el nivel de superficie de este tope de pozo local. Como se ilustra por la rejilla 540 en la FIG. 5, los topes de pozo locales dentro de un entorno de búsqueda predefinido del nodo de rejilla 542 puede o no puede incluir topes de pozo 518-526 dependiendo del tamaño del entorno de búsqueda predefinido. El entorno de búsqueda predefinido puede incluir cualquier forma geométrica, sin embargo, preferiblemente es predefinido por un radio del nodo de rejilla que forma un entorno de búsqueda circular.

En la etapa 410, todas las fallas que cruzan un vector entre el nodo de rejilla y cada tope de pozo local respectivo ubicado en la etapa 408 se identifican usando las fallas almacenadas en la etapa 402. Cada vector entre el nodo de rejilla y cada tope de pozo local respectivo ubicado en la etapa 408 representa un vector de intersección.

En la etapa 412, los valores de dislocación de falla ajustados en la etapa 406 son acumulados (agregados juntos) a lo largo de cada vector de intersección respectivo y el resultado representa un valor de dislocación de falla acumulado para cada vector de intersección respectivo. Debido a que cada vector de intersección corresponde con (interseca) un tope de pozo local separado y el valor de dislocación de falla acumulado representa los valores de dislocación de falla ajustados a lo largo de cada vector de intersección, cada valor de dislocación de falla acumulado puede ser diferente para cada tope de pozo local respectivo.

En la etapa 414 , el valor de cada tope de pozo local encontrado en la etapa 408 se ajusta con base en el valor de dislocación de falla acumulado que corresponde con cada vector de intersección para un tope de pozo local respectivo. Cada valor de tope de pozo local ajustado por lo tanto, es igual al valor del tope de pozo local respectivo más el valor de dislocación de falla acumulado para el vector de intersección que interseca este tope de pozo local respectivo. Como un resultado, los valores de tope de pozo local ajustados temporalmente remueven la dislocación de falla local.

En la etapa 416, la interpolación se realiza en el nodo de rejilla usando los valores de tope de pozo local ajustados de la etapa 414. Se piensa que la interpolación abarca cualquier técnica de interpolación bien conocida tal como, por ejemplo, krigeado, cokrigeado colocado, algoritmos de simulación y cosimulación colocada. Opcionalmente, se pueden usar otras técnicas de interpolación bien conocidas. De esta manera, los nuevos topes de pozo locales ajustados se interpolan sobre la rejilla.

En la etapa 418, el método 400 determina si hay otro nodo de rejilla dentro de la rejilla. Si no hay otro nodo de rejilla dentro de la rejilla, entonces el método 400 finaliza. Si hay otro nodo de rejilla dentro de la rejilla, entonces el método 400 regresa a la etapa 408 y repite las etapas 408, 410, 412, 414 y 416 para cada nodo de rejilla dentro de la rejilla.

Aunque la red de falla que se sitúa dentro de la superficie interpolada 500 en la FIG. 5 es complicada, el método 400 produce una superficie interpolada 500 con pocos artefactos alrededor del extremo de cada falla en comparación con otros métodos convencionales . Como se ilustra por la comparación de la FIG. 3 y FIG. 6, la cual ilustra una superficie interpolada . de acuerdo con el método 400 usando los mismos datos usados para producir la superficie interpolada en la FIG. 3, el método 400 proporciona una mejoría significativa en la reducción de artefactos indeseable. Por ejemplo, los artefactos 302 ubicados cerca de la terminación de la falla 304 en la superficie interpolada 300 no están presentes cerca de la terminación de la falla 604 en la superficie interpolada 600 de la FIG. 6. Este procedimiento, por lo tanto, representa una mejoría sustancial para interpolar superficies que tienen . una o más fallas.

Descripción del Sistema La presente invención se puede implementar a través de un programa de instrucciones ejecutable por computadora, tales como módulos de programa, generalmente referido como aplicaciones de software o programas de aplicación ejecutados por una computadora. El software puede incluir, por ejemplo, rutinas, programas, objetos, componentes, y estructuras de datos que realizan tareas particulares o implementar tipos de datos de resumen particulares. El software formas una interfaz para permitir que una computadora reaccione de acuerdo con una fuente de entrada. DecisionSpaceTM, la cual es una aplicación de software comercial comercializada por Landmark Graphics Corporation, se puede usar como una aplicación de int.erfaz para implementar la presente invención. El software también puede cooperar con otros segmentos de código para iniciar una variedad de tareas en respuesta a los datos recibidos en conjunto con la fuente de los datos recibidos. El software se puede almacenar y/o portar en cualquier variedad de medios de memoria tales como CD-ROM, disco magnético, memoria de burbuja y memoria semiconductora (por ejemplo, varios tipos de RAM o ROM) . Además, el software y sus resultados se pueden transmitir sobre una variedad de medios portadores tal como fibra óptica, cable metálico y/o a través de cualquiera de una variedad de redes tal como la Internet .

Además, aquellos expertos en el arte apreciarán que la invención se puede practicar con una variedad de configuraciones de sistema de computadora, incluyendo dispositivo portátiles, sistemas multiprocesadores , aparatos electrónicos basados en microprocesador o programables por el consumidor, minicomputadoras, computadoras centrales, y similares. Cualquier número de redes de computadora y sistemas de computadora es aceptable para el uso con la presente invención. La invención se puede practicar en ambientes de computación distribuida donde las tareas se realizan por dispositivos de procesamiento remoto que se enlazan a través de una red de comunicaciones. En un ambiente de computación distribuida, los módulos de programa se pueden ubicar en medios de almacenamiento de computadora locales y remotos incluyendo dispositivos de almacenamiento de memoria. La presente invención puede por lo tanto, ser implementada en conexión con varios hardware, software o una combinación de los mismos, en un sistema de computadora u otro sistema de procesamiento.

Con referencia ahora a la FIG. 7, se ilustra un diagrama de bloque de un sistema para implementar la presente invención en una computadora. El sistema incluye una unidad de computación, algunas veces referida como un sistema de. computación, que contiene memoria, programas de aplicación, una interfaz de cliente, una interfaz de vídeo y una unidad de procesamiento. La unidad de computación solamente es un ejemplo de un ambiente de computación adecuado y no se propone para sugerir cualquier limitación en cuanto al alcance del uso o funcionalidad de la invención.

La memoria principalmente almacena los programas de aplicación, los cuales también se pueden describir · como módulos de programa que contienen instrucciones ejecutables por computadora, ejecutados por la unidad de computación para implementar la presente invención descrita en la presente e ilustrada en las FIGS. 4-6. La memoria por lo tanto, principalmente incluye un módulo de interpolación de superficie, el cual realiza las etapas 402-418 ilustradas en la FIG. 4. Aunque DecisionSpace™ se puede usar para interconectarse con el módulo de interpolación de superficie para realizar las etapas 402-418, otras aplicaciones de interfaz se pueden usar en lugar de DecisionSpace™ o el módulo de interpolación de superficie se puede usar como una aplicación independiente. DecisionSpace™ proporciona una interfaz común para el módulo de interpolación de superficie y otras aplicaciones y/o módulos ilustrados en la FIG. 7. Permite que el usuario acceda a datos, verlos en visualizadores ID, 2D y 3D y realizar varios tipos de computaciones e interpretación de datos. El módulo de interpolación de superficie se configura actualmente para ser implementado a través del módulo Dynamic Geomodeling™ en DecisionSpace™ como se ilustra en la FIG. 7.

Aunque la unidad de computación se muestra que tiene una memoria generalizada, la unidad de computación típicamente incluye una variedad de medios leíbles por computadora. Por vía de ejemplo, y no limitación, el medio leíble por computadora puede comprender medios de almacenamiento de computadora. La memoria de sistema de computación puede incluir medios de almacenamiento de computadora en la forma de memoria volátil y/o no volátil tal como una memoria de solo lectura (ROM, por sus siglas en inglés) y memoria de acceso aleatorio (RAM, por sus siglas en inglés) . Ün sistema básico de entrada/salida (BIOS, por sus siglas en inglés) , que contiene las rutinas básicas que ayudan a transferir información entre los elementos dentro de la unidad de computación, tal como durante el arranque, típicamente se almacena en ROM. La RAM típicamente contiene datos y/o módulos de programa que son accesibles inmediatamente y/o actualmente se operan por la unidad de procesamiento. Por vía de ejemplo, y no limitación, la unidad de computación incluye un sistema operativo, programas de aplicación, . otros módulos de programa, y datos de programa.

Los componentes mostrados en la memoria también se pueden incluir en otros medios de almacenamiento de computadora removibles/no removibles, volátiles/no volátiles o se pueden implementar en la unidad de computación a través de la interfaz de programa de aplicación ("API" , por sus siglas en inglés) , la cual puede residir en una unidad de computación separada conectada a través de una red o sistema de computadora. Por ejemplo, solamente una unidad de disco duro puede leer de o escribir en medios magnéticos no volátiles, no removibles, una unidad de disco magnético puede leer de o escribir en un disco magnético no volátil, removible, y una unidad de disco óptico puede leer de o escribir en un disco óptico no volátil, removible tal como un CD ROM u otro medio óptico. Otros medios de almacenamiento de computadora removibles/no removióles, volátiles/no volátiles que se pueden usar en el ambiente de operación ejemplar pueden incluir, pero no se limitan a, casetes de cinta magnética, tarjetas de memoria flash, discos versátiles digitales, cinta de vídeo digital, RAM de estado sólido, ROM de estado sólido, y similares. Las unidades y sus medios de almacenamiento de computadora asociados discutidos anteriormente proporcionan almacenamiento de instrucciones leíbles por computadora, estructuras de · datos, módulos de programa y otros datos para la unidad de computación.

Un cliente puede introducir comandos e información en la unidad de computación a través de la interfaz de cliente, que pueden ser dispositivos de entrada tales como un dispositivo de señalización y teclado, comúnmente referido como un ratón, bola de seguimiento o almohadilla táctil. Los dispositivos de entrada pueden incluir un micrófono, palanca de mando, plato satelital, escáner, o similar. Estos y otros dispositivos de entrada son frecuentemente conectados a la unidad de procesamiento a través de un bus de sistema, pero se pueden conectar por otras estructuras de bus e interfaz, tal como un puerto paralelo o un bus serial universal (USB, por sus siglas en inglés) .

Un monitor u otro tipo de dispositivo visualizador se puede conectar al bus de sistema vía una interfaz, tal como una interfaz de vídeo. Una interfaz de usuario gráfica ("GUI" , por sus siglas en inglés) también se puede usar con la intérfaz de vídeo para recibir instrucciones de la interfaz de cliente y transmitir instrucciones a la unidad de procesamiento. Además del monitor, las computadoras pueden incluir otros dispositivos de salida periféricos tales como altavoces e impresora, los cuales se pueden conectar a través de una interfaz periférica de salida.

Aunque muchos otros componentes internos de la unidad de computación no se muestran, aquellos de experiencia ordinaria en el arte apreciarán que tales componentes y su interconexión son bien conocidos.

Se hace constar que con relación a esta fecha, el mejor método conocido por la solicitante para llevar a la práctica la citada invención, es el que resulta claro de la presente descripción de la invención.

Claims (20)

REIVINDICACIONES Habiéndose descrito la invención como antecede se reclama como propiedad lo contenido en las siguientes reivindicaciones :

1. Un método para crear una superficie que tiene una o más fallas, caracterizado porque comprende: a) estimar un nivel de superficie en cada lado de cada falla, una diferencia de elevación ente los niveles de superficie estimados en cada lado de una falla respectiva que representa un valor de dislocación de falla respectivo; b) ajustar cada valor de dislocación de falla usando un procesador de computadora; c) ubicar cada tope de pozo local dentro de un entorno de búsqueda predefinido de un nodo de rejilla dentro de una rejilla en la superficie; d) identificar cada falla que cruza un vector entre el nodo de rejilla y cada tope de pozo local respectivo, cada vector entre el nodo de rejilla y cada tope de pozo local respectivo que representa un vector de intersección; e) acumular cada valor de dislocación de falla ajustado a lo largo de cada vector de intersección respectivo, el cual representa un valor de dislocación de falla acumulado para cada vector de intersección respectivo; f) ajustar un valor para cada tope de pozo local con base en un valor de dislocación de falla acumulado respectivo; g) realizar la interpolación en el nodo de rejilla usando el valor ajustado para cada tope de pozo local; y h) repetir las etapas c) - g) para cada nodo de rejilla dentro de la rejilla.

2. El método de conformidad con la reivindicación 1, caracterizado porque adicionalmente comprende almacenar datos de entrada que comprenden topes de pozo y fallas, cada falla comprende uno o más segmentos de falla, cada tope de pozo se almacena en una estructura de datos de árbol KD y cada segmento de falla se almacena en una estructura de datos de árbol Quad.

3. El método de conformidad con la reivindicación 1, caracterizado porque el nivel de superficie en cada lado de cada falla se estima hallando cada tope de pozo en cada bloque de falla e interpolando el nivel de superficie a lo largo de cada lado de cada falla usando uno o más de los topes de pozo encontrados en cada bloque de falla respectivo.

4. El método de conformidad con la reivindicación 3, caracterizado porque la interpolación se realiza por krigeado, cokrigeado colocado, simulación o cosimulación colocada.

5. El método de conformidad con la reivindicación 1, caracterizado porque. cada valor de dislocación de falla se ajusta alisando los niveles de superficie estimados a lo largo de cada lado de cada falla que pertenece al mismo bloque de falla y ahúsando cada valor de dislocación de falla á cero en cada extremo distal de cada falla.

6. El método de conformidad con la reivindicación 1, caracterizado porque el entorno de búsqueda predefinido forma un entorno de búsqueda circular con un radio del nodo de rejilla.

7. El método de conformidad con la reivindicación 1, caracterizado porque cada valor de dislocación de falla ajustado se acumula agregando cada valor de dislocación de falla a lo largo de cada vector de intersección respectivo.

8. El método de conformidad con la reivindicación 1, caracterizado porque el valor para cada tope de pozo local se ajusta agregando el valor de cada tope de pozo local al valor de dislocación de falla acumulado respectivo.

9. El método de conformidad con la reivindicación 8, caracterizado porque el valor de dislocación de falla acumulado respectivo es el valor de dislocación de falla acumulado para el vector de intersección que interseca el tope de pozo respectivo.

10. El método de conformidad con la reivindicación 1, caracterizado porque la realización de la interpolación en el nodo de rejilla comprende krigeado, cokrigeado colocado, simulación o cosimulación colocada en el nodo de rejilla usando el valor ajustado para cada tope de pozo local.

11. dispositivo portador de programa transitorio caracterizado porque porta tangiblemente instrucciones ejecutables por computadora para crear una superficie que tiene una o más fallas, las instrucciones son ejecutables para implementar.- a) estimación de un nivel de superficie en cada lado de cada falla, una diferencia de elevación entre los niveles de superficie estimados en cada lado de una falla respectiva que representa un valor de dislocación de falla respectivo; b) ajuste de cada valor de dislocación de falla usando un procesador de computadora; c) ubicación de cada tope de pozo local dentro de un entorno de búsqueda predefinido desde un nodo de rejilla dentro de una rejilla en la superficie; d) identificación de cada falla que cruza un vector entre el nodo de rejilla y cada tope de pozo local respectivo, cada vector entre el nodo de rejilla y cada tope de pozo local respectivo que representa un vector de intersección; e) acumulación de cada valor de dislocación de falla ajustado a lo largo de cada vector de intersección respectivo, el cual representa un valor de dislocación de falla acumulado para cada vector de intersección respectivo; f) ajuste de un valor para cada tope de pozo local con base en un valor de dislocación de falla acumulado respectivo; g) realización de la interpolación en el nodo de rejilla usando el valor ajustado para cada tope de pozo local; y h) repetición de las etapas c) - g) para cada nodo de rejilla dentro de la rejilla.

12. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque adicionalmente comprende almacenar datos de entrada que comprenden topes de pozo y fallas, cada falla comprende uno o más segmentos de falla, cada tope de pozo se almacena en una estructura de datos de árbol KD y cada segmento de falla se almacena en una estructura de datos de árbol Quad.

13. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque el nivel de superficie en cada lado de cada falla se estima hallando cada tope de pozo en cada bloque de falla e interpolando el nivel de superficie a lo largo de cada lado de cada falla usando, uno o más de los topes de falla encontrados en cada bloque de falla respectivo.

14. El dispositivo portador de programa de conformidad con la reivindicación 13, caracterizado porque la interpolación se realiza por krigeado, cokrigeado colocado, simulación o cosimulación colocada.

15. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque cada valor de dislocación de falla se ajusta alisando los niveles de superficie estimados a lo largo de cada lado de cada falla que pertenece al mismo bloque de falla y ahusando cada valor de. dislocación de falla a cero en cada extremo distal de cada falla.

16. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque el entorno de búsqueda predefinido forma un entorno de búsqueda circular con un radio del nodo de rejilla.

17. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque cada valor de dislocación de falla ajustado se acumula agregando cada valor de dislocación de falla a lo largo de cada vector de intersección respectivo.

18. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque el valor para cada tope de pozo local se ajusta agregando el valor de cada tope de pozo local al valor de dislocación de falla acumulado respectivo.

19. El dispositivo portador de programa de conformidad con la reivindicación 18, caracterizado porque el valor de dislocación de falla acumulado respectivo es el valor de dislocación de falla acumulado para el vector de intersección que interseca el tope de pozo respectivo.

20. El dispositivo portador de programa de conformidad con la reivindicación 11, caracterizado porque la realización de la interpolación en el nodo de rejilla comprende krigeado, cokrigeado colocado, simulación o cosimulación colocada en el nodo de rejilla usando el valor ajustado para cada tope de pozo local.