Res2dinv-2.en.es (1) TRADUCIDO

Ver RES2DINV. 3.59 para Windows XP / Vista / 7

Rápido 2-D Resistividad e IP inversión utilizando el método de mínimos cuadrados

Wenner ( α, β, gamma), dipolo-dipolo, en línea polo polos, poledipole, ecuatorial dipolo-dipolo, compensado polo-dipolo,

Wenner-Schlumberger, gradiente y no convencionales arrays En tierra, el agua y encuestas de corte del pozo

Geoelectrico de imágenes 2D y 3D

SOFTWARE GEOTOMO Malasia www.geoelectrical.com de febrero de 2010

Este software es distribuido por:

Geotomo Software, 115 Cangkat Minden Jalan 5, Minden Heights, 11700 Gelugor, Penang, Malasia

Tel: +60 4 6574525, Fax: +60 4 6588437 correo electrónico: electrónico: [email protected], [email protected], [email protected] I nternet: www.geoelectrical.com

Tabla de contenido

Tema 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 dieciséis

17 18 19

Introducción

Página No.

operaciones de archivo de datos y formato de datos

1 3 6 7 9 10 11

Edición de los datos  datos 

dieciséis

Cambio de los ajustes del programa

19 25 33 37 42 43 43 44 44 46 50 51 52 55 58 60 62 63 76 80 81 82 92 96 98 104 106 109 112 116 119 127

Requisitos del sistema informático y la instalación protección de copia

Teoría contenido del CD

Usando el programa

opciones de inversión

modelado topográfico Mostrar los resultados de la inversión La configuración del programa  Ayuda Impresión de los pseudosecciones

Otros programas útiles Errores en 2-D inversión  Algunos ejemplos de campo Solución de problemas Expresiones Expresiones de gratitud

referencias  Apéndice A: encuestas encuestas dipolo-dipolo-dipolo y el polo

 Apéndice B: La matriz Wenner-Schlumberger  Apéndice C: Las resistividades resisti vidades de fijación

D Apéndice: La matriz Ecuatorial dipolo-dipolo  Apéndice E: estudios submarinos  Apéndice F: inversión inversión de datos IP  Apéndice G: Nuevas opciones de modo modo de proceso por lotes  Apéndice H: Velocidad de visualización visualizació n de gráficos gráfi cos

 Apéndice I: inversión inversión de datos cruzada pozo

J Apéndice: electrodos remotos K Apéndice: datos sólidos y el modelo de inversión  Apéndice L: L : matrices no convencion c onvencionales ales o generales general es

 Apéndice M: tipos tipos de matriz en RES2DINV  Apéndice N: guardar guardar los resultados en formato SURFER SURFER

O Anexo: Modelo refinamiento  Apéndice P: Combinación C ombinación de d e archivos 2-D

 Apéndice Q: detección de bordes Be drock

 Apéndice I: Time-lapse resistividad inversión  Apéndice S: La incorporación de las estimaciones estimaciones de ruido de datos

 Apéndice T: La incorporación de límites de capas conocidos 129 Apéndice Apéndice T: información topografía especial

Tabla de contenido

Tema 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 dieciséis

17 18 19

Introducción

Página No.

operaciones de archivo de datos y formato de datos

1 3 6 7 9 10 11

Edición de los datos  datos 

dieciséis

Cambio de los ajustes del programa

19 25 33 37 42 43 43 44 44 46 50 51 52 55 58 60 62 63 76 80 81 82 92 96 98 104 106 109 112 116 119 127

Requisitos del sistema informático y la instalación protección de copia

Teoría contenido del CD

Usando el programa

opciones de inversión

modelado topográfico Mostrar los resultados de la inversión La configuración del programa  Ayuda Impresión de los pseudosecciones

Otros programas útiles Errores en 2-D inversión  Algunos ejemplos de campo Solución de problemas Expresiones Expresiones de gratitud

referencias  Apéndice A: encuestas encuestas dipolo-dipolo-dipolo y el polo

 Apéndice B: La matriz Wenner-Schlumberger  Apéndice C: Las resistividades resisti vidades de fijación

D Apéndice: La matriz Ecuatorial dipolo-dipolo  Apéndice E: estudios submarinos  Apéndice F: inversión inversión de datos IP  Apéndice G: Nuevas opciones de modo modo de proceso por lotes  Apéndice H: Velocidad de visualización visualizació n de gráficos gráfi cos

 Apéndice I: inversión inversión de datos cruzada pozo

J Apéndice: electrodos remotos K Apéndice: datos sólidos y el modelo de inversión  Apéndice L: L : matrices no convencion c onvencionales ales o generales general es

 Apéndice M: tipos tipos de matriz en RES2DINV  Apéndice N: guardar guardar los resultados en formato SURFER SURFER

O Anexo: Modelo refinamiento  Apéndice P: Combinación C ombinación de d e archivos 2-D

 Apéndice Q: detección de bordes Be drock

 Apéndice I: Time-lapse resistividad inversión  Apéndice S: La incorporación de las estimaciones estimaciones de ruido de datos

 Apéndice T: La incorporación de límites de capas conocidos 129 Apéndice Apéndice T: información topografía especial

para la exploración de minerales

V Apéndice: Trazado de métodos m étodos para la matriz de gradiente W Apéndice: La matriz polo-dipolo de compensación

 Apéndice X: coordenadas coordenadas globales o GPS Nota especial OhmMapper

Renuncia  Apoyo Qué hay de nuevo

132

133 138 140 143 144 144 145

 Aviso de derechos de autor RES2DINV los derechos de autor • Geotomo Software, 2000-2010. Todos los derechos reservados.

Ninguna parte de este documento puede ser reproducida sin el permiso por escrito de G eotomo software.

1

1

Introducción

RES2DINV es un programa informático que va a determinar automáticamente un modelo de resistividad de dos dimensiones (2-D) para el subsuelo para los datos obtenidos a partir de encuestas de imagen eléctricos (Griffiths y Barker 1993). Dado que es un programa basado en Windows, todas las tarjetas de gráficos e impresoras compatibles con Windows son compatibles de forma automática. Se ha probado con los modos de pantalla de vídeo de hasta 1600 x 1200 píxeles y 16 millones de colores.

La Figura 1 muestra un ejemplo de la disposición de electrodos y la secuencia de medición que se puede utilizar para una encuesta 2-D de imagen eléctrica. Este programa está diseñado para invertir grandes conjuntos de datos (con alrededor de 200 a 21.000 puntos de datos) recogidos con un sistema con un gran número de electrodos (alrededor de 25 a 16.000 electrodos).

El modelo 2-D utilizado por el programa de inversión, que consiste en un número de bloques rectangulares, se muestra en la Figura 2. La disposición de los bloques está vagamente ligada a la distribución de los puntos de datos en el pseudosección. La distribución y el tamaño de los bloques se genera automáticamente por el programa utilizando la distribución de los puntos de datos como una guía aproximada. La profundidad de la fila inferior de bloques se establece para que sea aproximadamente igual a la profundidad equivalente de investigación (Edwards 1977) de los puntos de datos con el mayor espaciamiento de los electrodos. La encuesta se lleva a cabo por lo general con un sistema en el que los electrodos están dispuestos a lo largo de una línea con una separación constante entre los electrodos adyacentes. Sin embargo, el programa también puede manejar los conjuntos de datos con una separación entre electrodos no uniforme.

Una subrutina modelado directo se utiliza para calcular los valores de resisti vidad aparente, y una técnica de optimización no lineales de mínimos cuadrados se usa para la rutina de inversión (deGroot-Hedlin y Constable 1990, Loke y 1996a Barker). El programa es compatible con las técnicas de modelado con visión de elementos finitos y diferencias finitas. Este programa puede ser utilizado para las encuestas utilizando el Wenner, polo polos, dipolo-dipolo, polo-dipolo, matrices ecuatoriales dipolo-dipolo (rectangulares) Wenner-Schlumberger y. Además de estas matrices comunes, el programa soporta incluso las matrices no convencionales con un número casi ilimitado de posibles configuraciones de electrodos! Puede procesar pseudosecciones con hasta 16000 electrodos y 21000 puntos de datos en una sola vez en un equipo con 1 GB de RAM. La mayor distancia entre los electrodos puede ser de hasta 36 veces la separación más pequeña utilizada en un solo conjunto de datos. Los límites de datos de programa se extenderá en el futuro a medida que se encontraron con grandes conjuntos de datos de campo.

 Además de las encuestas normales llevados a cabo con los electrodos en la superficie del terreno, el programa también es compatible con estudios submarinos y transversal del pozo!

2

Figura 1. Secuencia de mediciones para construir una pseudosección utilizando una configuración de múltiples electrodos encuesta controlado por ordenador.

Figura 2. Disposición de los bloques utilizados en un modelo junto con los puntos de datos en la pseudosección.

3

2

Requisitos del sistema informático y la instalación

Este programa está diseñado para ejecutarse en IBM PC microcomputadoras compatibles a través de los sistemas operativos XP / Vista / 7 de Windows. Debido a la potencia de cálculo necesaria para el modelado hacia adelante y mínimos cuadrados subrutinas 2-D, está escrito para los sistemas informáticos de 32 bits con al menos 256 megabytes de RAM. En un equipo Pentium 4, se tardará menos de un minuto para invertir el conjunto de datos obtenidos a partir de una encuesta c on cerca de 100 electrodos.

Requisitos mínimos: disco duro con al menos 256 megabytes de espacio libre, y la placa base con al menos 256 megabytes de RAM. sistema de gráficos en color SVGA con una resolución de pantalla de al menos 800 por 600 píxeles y 256 colores. / Sistema operativo Windows XP Vista / 7. Pentium III o Pentium 4 (o compatible) CPU. Usted debe ser capaz de procesar el mayor co njunto de datos de ejemplo, DUFUYA.DAT que tiene 297 electrodos, con el sistema informático anteriormente. Sin embargo, se r ecomienda que utilice un sistema con 512 (o más) MB de memoria y espacio de disco duro l ibre al menos 1 GB para procesar grandes co njuntos de datos con más de 400 electrodos y 3000 puntos de datos. Este es un programa de Windows de 32 bits que pueden acceder hasta 4 gigabytes de memoria. Cuanta más memoria RAM y el espacio del disco duro que tiene, cuanto mayor sea el conjunto de datos que se pueden procesar con este programa.

También se recomienda el uso de los 1024 por 768 (para 17 monitores de pulgada) o el 1280 por 1024 (para los monitores de 19 pulgadas) con 256 (o superior) Modo de gráficos colores XVGA. Se puede cambiar el modo gráfico haciendo clic en 'Inicio' y luego 'Configuración' segui do de 'Control Panel' y luego 'Display'.

Si usted había comprado el software con un CD, el siguiente gestor de configuración debe aparecer automáticamente al insertar el CD.

Haga clic en el elemento adecuado con el fin de instalar el programa RES2DINV; seguido de los programas RES3DINV, RES2DMOD y RES3DMOD. Si usted tiene la

4

Keylok dongle USB, haga clic en el 'Keylok USB Driver' para instalar el controlador para el dispositivo de seguridad.

Si usted había descargado los programas de Internet, tendrá que instalar los programas m anualmente. El paquete RES2DINV viene en un solo archivo de instalación SETUP .EXE comprimido que es un programa de instalación basado en Windows. Para instalar el programa, haga clic en Inicio, y luego Configuración seguido por el Panel de control y luego en Agregar / quitar programas. Después de instalar el programa RES2DINV.EXE, debe hacer clic en la o pción “Ejecutar JACOBWIN” undret él menú

'Archivo' para crear varios archivos de soporte. Sólo tiene que ejecutar esta opción una vez después de instalar el programa RES2DINV.

Para la versión completa de este software, es necesario instalar un controlador de sistema para el dongle. Si es necesario instalar el controlador manualmente, siga los pasos a continuación. un).

Para los usuarios con la tecla Keylok (puerto paralelo o USB)

Por favor, ejecute el programa INSTALL.EXE en la carpeta KEYLOK_SYSTEM_DRIVER para los pilotos Keylok. Este programa instalará el controlador apropiado para su sistema operativo en su disco duro del ordenador .. Si está utilizando el Keylok USB llave, por favor siga los pasos a continuación.

1). Si el RES2DINV y el software RES3DINV ya están instalados en su ordenador, tendrá que reemplazar los archivos RES2DINV.EXE y RES3DINV.EXE con las últimas versiones en el CD. Es necesario copiar los archivos RES2DINV.EXE y RES3DINV.EXE de la carpeta Win_Updates a las carpetas correspondientes en su ordenador (normalmente C: \ RES2DINV y C: \ RES3DINV). 2). Si aún no se han instalado el software de RES2DINV y RES3DINV, ejecute los archivos SETUP.EXE que se encuentran en las carpetas RES2DINV y RES3DINV en el DISCOS COMPACTOS.

3). En primer lugar, asegúrese de que el dispositivo USB está no y conectado al equipo. Ahora tendrá que instalar el controlador de software para el dispositivo de seguridad Keylok ejecutando el programa INSTALL.EXE (mediante el comando install.EXE / B) en la carpeta KEYLOK_SYSTEM_DRIVER en el CD. Tenga en cuenta que en Windows NT / 2000 / XP, debe tener privilegios de administrador para instalar el controlador. 4). Después de instalar el controlador, apague el ordenador y conectar la llave USB. Se recomienda conectar el dispositivo USB a la primero puerto USB del ordenador. Algunos sistemas informáticos mayores tienen dificultades para detectar que en el segundo puerto USB. Ahora encienda el ordenador. Windows mostrará un mensaje que indica que ha encontrado un nuevo hardware y que a continuación, debe copiar automáticamente el controlador para el dispositivo de seguridad. 5). Intente ejecutar el programa RES2DINV y RES3DINV, y que debe ser capaz de detectar el dongle. Si los programas todavía no puede detectar el dispositivo de seguridad, tendrá que instalar el controlador manualmente. Haga clic en el Panel d e control en Configuración y, a continuación, seleccione la opción“Agregar  nuevo hardware”.  A continuación , Windows bu sca

5

nuevos dispositivos Plug and Play. Haga clic en el “Sin el dispositivo no está en la lista” y luego el “No deseo seleccionar el hardware de una lista” opciones. Cerca de la parte inferior cerca de la lista de hardware, probablemente se mostrará el “Dongle USB”. Haga clic en él, y luego la opción “Utilizar disco” en el cuadro de diálogo siguiente. A continuación, haga clic en “Examinar” y cambiar la carpeta a la carpeta KEYLOK_SYSTEM_DRIVER en el CD. A continuación, se muestra una lista de archivos con la extensión “inf”. Para Windows 98 / Me seleccione el archivo “usbkey98.inf”, para Windows NT / 2000 / XP, seleccione el archivo “Usbkey2K.inf”. Después de que haga clic en “Aceptar”, y una caja de diálogo “Copia de archivos” se debe mostrar

y Windows debe copiar los controladores.

Después de copiar los controladores, es necesario reiniciar el equipo y los programas RES2DINV / RES3DINV continuación, debería ser capaz de detectar el dispositivo de seguridad. Más información acerca de los controladores del sistema se puede obtener en la página web siguiente: www.keylok.com

Vista: Para Windows Vista, es necesario utilizar el controlador que se encuentra en el archivo en la carpeta InstallVistaUSB.zip “KEYLOK_SYSTEM_DRIVER \ Vista Driver” en el CD. Las instrucciones se dan en el archivo “Para instalar Dongle USB en Vista.pdf”.

segundo). Para los usuarios con la tecla Rainbow Technologies Sentinel SuperPro

Esta clave se distribuye con versiones muy antiguas del software, y ya no s e distribuye pero aún está soportado en las actualizaciones de los programas. Ejecutar el RainbowSSD5.39.2.exe archivo

en

la RAINBOW_SYSTEM_DRIVER

subdirectorio en el CD-ROM. Este programa instalará automáticamente el controlador del sistema p ara el dispositivo de seguridad. Después de instalar el controlador, reinicie el equipo para que el conductor tenga efecto. En caso de problemas, utilice el programa SuperProMedic.exe para comprobar el estado del dispositivo de protección. Más información acerca de los controladores del sistema se puede obtener en el siguiente sitio web: www.rainbow.com.

Prácticamente todos los equipos tienen un sistema de administración de energía incorporada que se ralentiza y finalmente se apaga el sistema informático si no se tiene acceso al teclado / ratón o el disco duro después de un cierto límite de tiempo. Esto puede interferir con el funcionamiento del programa de RES2DINV si va a invertir una gran cantidad de datos ya que este programa pasa la mayor parte del tiempo de cálculo utilizando sólo la CPU. Windows también tiene una función de protector de pantalla incorporado que reemplaza el contenido de la pantalla con un protector de pantalla. Antes de ejecutar el programa de RES2DINV, tendrá que desactivar tanto la administración de energía y protector de pantalla programas. Para desactivarlos, es necesario seguir los siguientes pasos. 1). Poco después de que el ordenador está encendido, por lo general se ejecuta una comprobación de la memoria. Durante este tiempo, es necesario ir al programa de configuración de CMOS en la BIOS del ordenador. Para la mayoría de las computadoras, esto se hace pulsando la tecla "Supr". Próximamente se muestra la pantalla de configuración del CMOS. Vaya a la sección que se refiere a la adminis tración de energía, y desactivar las opciones de administración de energía para el equipo.

6

2). Después se inicia Windows, vaya al panel de control. Haga clic en el icono Administración de energía o Opciones de energía. Aquí probablemente mostrará el número de minutos que deben de Windows se apagará automáticamente el equipo si no se utiliza el ratón o el teclado. Desactivar esta característica. 3). A continuación, todavía dentro del Panel de control, haga clic en el icono Pantalla. En el programa de visualización, haga clic en la opción de protector de pantalla y desactivar el protector de pantalla. Esto asegurará que Windows no se ejecuta automáticamente un protector de pantalla. Después de desactivar la administración de energía y protector de pantalla de opciones, puede cambiar físicamente apagar el monitor si el equipo se deja en marcha durante la noche. Hay preferiblemente no debería haber demasiados programas en segundo plano que se ejecutan antes de ejecutar este programa. Los programas en segundo plano activos generalmente se muestran en la barra d e 'Inicio' en la parte inferior de la pantalla. Por el cierre de los otros programas, más memoria estará disponible para este programa. Esto reducirá la memoria de intercambio de discos duro q ue ralentiza el programa. Antes de llevar a cabo la inversión de un gran conjunto de datos con más de un millar de puntos de datos, por favor c errar otros programas activos. Además, el uso de la sub-opción “Establecer la prioridad del programa” en la opción “File” del menú principal, establecer la prioridad a programas de alta para asegurarse de que no se interrumpe por otros programas durante la inversión del conjunto de datos.

En Windows, un programa de la papelera de reciclaje se utiliza a veces para guardar archivos borrados recientemente. Esta característica podría llegar a ser una molestia cuando se invierte grandes conjuntos de datos como el programa utiliza un número de archivos temporales. Aunque el programa elimina los archivos, que podría terminar en la papelera de reciclaje, y rápidamente llenar su disco duro! El programa RES2DINV, por defecto, utilice la unidad de disco duro con la mayor cantidad de espacio libre como la unidad de memoria intermedia. Se recomienda que desactive la papelera de reciclaje para esta unidad. Esto se puede hacer haciendo clic en el icono de la papelera de reciclaje con el botón derecho del ratón, seguido de la opción Propiedades. A continuación, puede configurar los ajustes para cada unidad de disco duro.

Se recomienda que quite el disco CD, si está presente, desde la unidad de CD- ROM cuando se utiliza este programa. Debido a la naturaleza multi-tarea de estos sistemas, a veces el programa en el disco podría ser lanzado mientras RES2DINV está ocupado invirtiendo un gran conjunto de datos! Para todas las versiones de Windows, se recomienda ejecutar el programa SCANDISK vez en cuando para comprobar el estado del disco duro, y también el programa de desfragmentación para reducir l a fragmentación de archivos.

3

protección de copia

Una llave de hardware USB (dongle) se suministra con la versión registrada del programa. Después de instalar el controlador, conecte esta clave a un puerto USB del ordenador. La versión de demostración del programa no requiere una llave de hardware, pero algunas de las características importantes del programa no se puede utilizar.

7

Las versiones anteriores de este software utilizan un dispositivo de seguridad que se une al puerto de impresora paralelo del ordenador. Se puede conectar el conector de 25 pines de la impresora en el otro extremo de la llave. El programa debe ser capaz de detectar el dispositivo de seguridad si el puerto paralelo de impresora se ajusta a las especificaciones de IBM PC. Si no es capaz de encontrar la llave, el puerto de la impresora no completamente compatible con IBM PC. En este caso, después de conectar el dongle al ordenador, conectar una impresora al otro extremo del dongle. Encienda la impresora e imprimir algo para confirmar que el puerto de la impresora está funcionando. Deja la impresora y, a continuación, intentar ejecutar el programa RES2DINV.EXE nuevo. Si usted todavía no es capaz de resolver el problema, póngase en contacto con la empresa Rainbow Technologies o Keylok por correo electrónico.

Muchos de los nuevos ordenadores portátiles no tienen un puerto paralelo. Si usted tiene el tipo de puerto paralelo del dispositivo de seguridad, y la necesidad de utiliz ar el programa en un ordenador portátil tal, se recomienda que el intercambio con un dongle USB. Por favor, póngase en contacto con [email protected] sobre el procedimiento para el intercambio de la llave.

Teoría 4 La rutina de inversión utilizados por el programa se basa en el m étodo de mínimos cuadrados smoothnessconstrained (deGroot-Hedlin y Constable 1990, Sasaki 1992). El método de mínimos cuadrados suavidad con restricci ones se basa en la siguiente ecuación

( J T J + u F) d = J T gramo

(1)

dónde F = f x F x T + F z F z T F x = filtro de planeidad horizontal F z = filtro de planeidad vertical,

J = matriz de derivadas parciales u = factor de amortiguamiento d = modelo de vector de perturbación g = vector de discrepancia Una ventaja de este método es que los filtros de factores y de planeidad de amortiguación se pueden ajustar para adaptarse a diferentes tipos de datos. Una descripción detallada de las distintas variaciones del método de mínimos cuadrados suavidad limitados se puede encontrar en las notas clase particular libre por Loke (2001).

El programa es compatible con una nueva implementación del método de mínimos cuadrados basado en una técnica cuasi-Newton optimización (Loke y Barker, 1996a). Esta técnica es significativamente más rápido que el método de mínimos cuadrados convencionales para

8

grandes conjuntos de datos y requiere menos memoria. También puede utilizar el método de Gauss-Newton convencional en este programa. Es mucho más que el método cuasi-Newton, pero en las zonas con grandes contrastes de resistividad de mayor que 10: 1, da resultados ligeramente mejores. Una tercera opción en este programa es usar el método GaussNewton para los primeros 2 o 3 i teraciones, después de lo cual se utiliza el método cuasi-Newton. En muchos casos, esto proporciona el mejor compromiso (Loke y Dahlin 2002). Debido a las mejoras en el código del programa y los PC se r ecomienda que la opción de utilizar el método de Gauss-Newton debe ser el método por defecto, en especial para el modelo interpretación final, como para la mayoría de los conjuntos de datos, la inversión de datos tomar sólo unos minutos en el PC moderna Los sistemas.

El modelo 2-D utilizado por este programa divide el subsuelo en un número de bloques rectangulares (Figura 2). El propósito de este programa es determinar la resistividad de los bloques rectangulares que producirán un pseudosección resistividad aparente que está de acuerdo con las medidas reales. Para las matrices de Wenner y Schlumberger, el espesor de la primera capa de bloques se establece en 0,5 veces la separación entre electrodos. Para el polo polos, dipolo-dipolo y arrays polo-dipolo, el espesor se ajusta a aproximadamente 0,9, 0,3 y 0,6 veces la distancia entre los electrodos respectivamente. El espesor de cada capa más profunda posterior normalmente se aumentó en un 10% (o 2 5%). Las profundidades de las capas también se pueden cambiar m anualmente por el usuario.

El método de optimización intenta básicamente para reducir la diferencia entre los valores de resistividad aparente calculados y medidos mediante el ajuste de la resistividad de los bloques de modelo. Una medida de esta diferencia está dada por el error de la raíz-meansquared (RMS). Sin embargo, el modelo con el error RMS bajo posible a veces puede mostrar variaciones grandes y poco realistas en los valores de resistividad modelo y no siempre puede ser la "mejor" modelo a partir de un punto de vista geológico. En general, el enfoque más prudente es elegir el modelo en la iteración después de lo cual el error RMS no cambia significativamente. Esto ocurre generalmente entre la tercera y iteraciones 5º.

Más información sobre el método de inversión se puede encontrar en los libres “Notas Tutorial” en la imagen eléctrica (Tutorial: encuestas de imagen eléctrica 2-D y 3-D) que está disponible en el sitio web www.geoelectrical.com.

9

5

contenido del CD

La carpeta RES2DINV en el CD contiene el principal archivo de instalación setup.exe. Después de instalar el archivo de instalación, usted debe tener los siguientes archivos en el subdirectorio que contiene el programa RES2DINV.

RES2DINV.EXE JACOBWIN.EXE REGISTER.TXT

programa de inversión principal

LANDFILL.DAT

Un archivo de datos de ejemplo de campos para la matriz Wenner obtenido por una

Programa para crear algunos archivos de soporte. archivo de registro encuesta con 50 electrodos. Este ejemplo muestra una aplicación interesante de imágenes 2-D eléctrica para trazar una pluma de contaminación de un vertedero (Niederleithinger,

1994).

GRUNDFOR.DAT

Otro ejemplo de campo para la matriz Wenner, pero con una distribución más irregular de los puntos de datos.

ODARSLOV.DAT

Un conjunto de datos bastante grande recogió sobre un di que de alta resistividad (array Wenner).

ROMO.DAT

Otro conjunto de datos bastante grande (array Wenner).

DUFUYA.DAT

Este es un gran conjunto de datos con cerca de 300 el ectrodos y más de 1200 puntos de datos (matriz Wenner). En un 266 Mhz Pentium II que tomó cerca de 23 minutos para invertir este conjunto de datos, donde la matriz jacobiana se calcula de nuevo después de cada iteración.

GLADOE2.DAT BLOCKWEN.DAT

Un ejemplo de archivo de datos con información topográfica. archivo de datos de prueba de entrada de gama Wenner con unos pocos puntos de datos malos.

BLOCKDIP.DAT BLOCKTWO.DAT

Ejemplo archivo de datos de entrada para array dipolo-dipolo.

RATHCRO.DAT

archivo de datos de matriz Wenner de un estudio arqueológico que

Ejemplo archivo de datos de entrada para array polo polos.

también contiene información topográfica.

PIPESCHL.DAT

datos de ejemplo de campos establecidos para la matriz Wenner-Schlumberger.

WATER.DAT LAKELELIA.DAT

Ejemplo de archivo de datos para una encuesta bajo el agua.

MODEL101.DAT DIPOLEN5.DAT

Un moderadamente amplio conjunto de datos de prueba.

BLUERIDGE.DAT

datos Ejemplo dipolo- dipolo fijados con diversas “a” y

Ejemplo de un conjunto de datos de estudio de campo bajo el agua.

Ejemplo archivo de datos dipolo-dipolo con valores no enteros para el “ norte ”Factor.

combinaciones “n”.

WENSCHN5.DAT

datos Ejemplo Wenner-Schlumberger establecidos con valores no enteros para el “ norte ”Factor.

PDIPREV.DAT

Ejemplo de archivo de polo-dipolo de datos con el "adelante" y la disposición "inversa" de los electrodos.

10

POLDPIN5.DAT

datos Ejemplo polo-dipolo establecidos con valores no enteros para el “ norte ”Factor.

datos de campo OHMMAPPER.DAT set de ejemplo de un sistema de topografía móvil KNIVSAS.DAT Ejemplo conjunto de datos utilizando la matriz de gradiente

ITALY_CAVE.DAT

Datos de ejemplo establecido en formato de datos de l a matriz general, pero con la topografía en una lista separada

IPMODEL.DAT

IP con el conjunto de datos con los valores de cargabilidad

IPSHAN.DAT

El campo de datos IP encuesta con los valores establecidos PFE

IPMAGUSI.DAT IPKENN.DAT

datos de la encuesta de campo IP establecen con valores del factor de metal.

BOREHOLE.DAT BOREHIP.DAT

Ejemplo conjunto de datos a través del pozo de sondeo.

BORELANC.DAT

conjunto de datos de campo a través del pozo.

BOREDIFF.DAT

Ejemplo datos transversal del pozo set donde los electrodos en los dos poz o de

datos IP de estudio de campo ocupados con valores de ángulo de fase Ejemplo datos transversal del pozo establecidos con las mediciones de IP.

sondeo no están en las mismas profundidades.

BORERES.DAT

Ejemplo datos transversal del pozo establecidos con medidas dadas como valores de resistencia.

BOREHOLE_TOPO. DAT

Ejemplo datos de perforación establecen con topografía

TIMELAPS.DAT MIXED.DAT

conjunto de datos Ejemplo de lapso de tiempo.

CROMER02.DAT RES2DINV.HLP RES2DINV.CNT

conjunto de datos matriz de gradiente.

RES2DINV.INI

archivo de configuración para el programa RES2DINV.

RES2DINV.IN2

archivo de configuración alternativa.

LEAME.TXT

La información más reciente acerca de la versión actual del programa de inversión

datos set de ejemplo en formato general matriz con mezcla de mediciones con Wenner-Schlumberger y arrays dipolo -dipolo.

archivo de ayuda para este programa. Un archivo de soporte para el archivo de ayuda.

de 2-D.

RESIS.BTH

6

 Archivo de ejemplo para la opción de modo de procesamiento por lotes

Usando el programa

Después de ejecutar el programa de instalación SETUP.EXE, y el programa JACOBWIN.EXE,

la principal RES2DINV.EXE programa y los archivos de soporte (GRADWEN, GRADTWO y GRADDIP) y archivos de datos de ejemplo deben estar en el m ismo subdirectorio en el disco duro. Por favor, compruebe el estado del disco duro m ediante los programas de CHKDSK y ScanDisk o programas similares en paquetes de utilidades de software como Norton Utilities.

Para ejecutar el programa de inversión de resistividad, haga clic en el icono de RES2DINV. El programa comprobará primero el sistema informático para garantizar que cuenta con los recursos necesarios que este programa requiere. Se comprobará la memoria disponible y

11

espacio de disco duro. Si el programa muestra un aviso, debe dejar de hacerlo desde el programa y hacer los cambios necesarios. Después de comprobar la configuración del equipo, el programa mostrará entonces la siguiente barra de menú principal en la parte superior de la pantalla.

Puede seleccionar una opción haciendo clic con el cursor del ratón. Si esta es la primera vez que utiliza el programa, trata de leer en el archivo de datos de resistividad LANDFILL.DAT proporcionado con el programa. A continuación, seleccione la opción de inversión para llevar a cabo una inversión del conjunto de datos. La Figura 3 muestra un ejemplo de los resultados de la inversión de este conjunto de datos. Después de leer un conjunto de datos, usted debe comprobar en general, para los malos puntos de datos utilizando el “ mi datos dit”opción. En las siguientes secciones, se da una descripción más detallada de cada elemento en el menú anterior.

7

operaciones de archivo de datos y formato de datos

Cuando se selecciona la opción “Archivo”, se muestran las siguientes opciones de menú.

12

R archivo de datos EAD - Cuando se selecciona esta opción, la lista de archivos en el directorio actual que tiene se mostrará una extensión de DAT. Se puede utilizar el ratón o el teclado para seleccionar el archivo correspondiente, o para cambiar la carpeta. Los valores de resistividad aparente primero se deben escribir en un archivo de texto. Se puede utilizar cualquier editor de texto de uso general, tales como el programa Bloc de notas de Windows. Los datos s e disponen de una manera  ASCII delimitado donde una coma o un espacio en blanco o un LF / CR se utiliza para separar los diferentes elementos de datos numéricos. El programa requiere que los datos que se van dispuestos en un formato determinado. Si hay un problema en el funcionamiento de este programa, una posible causa es que los datos de entrada se dispone en un formato incorrecto. Un ejemplo de un conjunto de datos de entrada se da en el LANDFILL.DAT archivo. La sección inicial de este archivo con una descripción del formato de datos se da en el recuadro de observaciones a continuación.

archivo LANDFILL.DAT

comentarios

ESTUDIO DE VERTEDERO

| Línea 1- Nombre de la línea de levantamiento

3.0 1

| Línea 2 - Más pequeño espacio entre electrodos

334 1

| Línea 4 - Número total de puntos de datos

| Línea 3 - tipo Array (Wenner = 1, Palo polos = 2, | dipolo-dipolo = 3, polo-dipolo = 6, Schlumberger = 7 | Ecuatorial dipolo dipolo - 8)

| Línea 5 - Tipo de p osición en x de puntos de datos. | Introduzca 0 si la ubicación del primer electrodo en la matriz utilizada para | indicar la posición en x. Introduzca 1 si el punto medio de la | a rray (es decir, localización de punto de datos en pseudosección) | se utiliza.

0

| Línea 6 - Bandera de datos IP (introduzca 0 para los datos de resistividad | solamente)

4,50 3,0 84,9

| Línea 7 - La posición en x, distancia entre electrodos y | valor de resistividad aparente medido para la primera | punto de datos.

7,50 3,0 62,8

| Línea 8 - La posición en x, distancia entre electrodos y | valor de resistividad aparente medido para la segunda | punto de datos.

Las líneas siguientes dan información similar para los otros puntos de datos. Esto es seguido por cuatro 0 de banderas que están reservados para otras funciones. los dipolo-dipolo, polo-dipolo y Schlumberger conjuntos de datos de la matriz tienen un formato ligeramente diferente, ya que un parámetro adicional, el factor de separación dipolo, es necesario. Por favor refiérase a Apéndices A y segundo para una descripción de los formatos de datos para estas matrices. El BLUERIDGE.DAT archivo es un ejemplo de datos establecidos para la matriz dipolo-dipolo.

13

Figura 3: El pseudosección resistividad aparente observados y calculados para los datos LANDFILL.DAT establecidos junto con un modelo obtenido por el programa de inversión.

Tenga en cuenta que el programa asume que la posición en x de los puntos de datos se incrementa desde la izquierda a la parte derecha de la pseudosección. Si los X-localizaciones en el conjunto de datos están dispuestos en la otra dirección, aparecerá la pseudosección haber sufrido una inversión de izquierda a derecha cuando se visualiza en la pantalla!

El GLADOE2.DAT archivo es un ejemplo con datos topográficos. La sección final de este archivo c on una descripción del formato de los datos topográficos es la siguiente: -

archivo GLADOE2.DAT

comentarios

237 2 39,207 203 2 14,546 227 2 31,793 233 2 30,285 2 26

| Últimos cuatro puntos de datos

- 100 33 - 40 34,5

| coordenadas horizontales y verticales de primera,

| con x-ubicación del punto de datos, | separación de los electrodos y

| valores de resistividad aparente medidos | marca de datos de topografía. Si no hay datos de topografía, coloque 0 aquí. | Número de puntos de datos de topografía

| Topografía segundo punto de datos

14

- 20 35,0

| Esto es seguido por los datos similar para

0 35.209 10 34.904 20 33.678

| los puntos de datos restantes topografía

30 32.264 40 31.191 50 30.560 60 29.905 70 29.200 80 28.774 90 28.569 100 28.473 110 28.553 120 28,462 130 28,486 140 28.601 150 28.652 160 28,825 170 28,996 180 29,284 190 29.643 200 29.563 270 30,0

|

300 33

| Última topografía punto de datos

2 0 0

| El número de puntos de datos topografía con el primer electrodo

0 0

| | | | | | | | | | | | | | | | | | | |

| Introduzca cuatro 0s en el

| final del archivo | |

Tenga en cuenta que los datos topográficos se coloca inmediatamente después de los puntos de datos de resistividad aparente. El primer elemento es una bandera para indicar si el archivo contiene datos de topografía. Si no hay datos de topografía, su valor es 0. Introduzca 1 o 2, si los datos topográfica está presente. En el caso en que se dan las c oordenadas horizontales y verticales reales de puntos a lo largo de la línea de sondeo, introduzca 1. Incluso si las distancias horizontales reales se dan en la sección de datos de topografía, todavía se debe utilizar el x-distancia a lo largo de la su perficie del suelo en el sección de datos de resistividad aparente. En la mayoría de las encuestas de las distancias de los puntos de la superficie del suelo, y no es cierto distancias horizontales, son en realidad mide con una cinta. En este caso, introduzca un valor de 2 para el indicador de datos de topografía. Esto es seguido por el número de puntos de datos topográficos.

No es necesario medir la elevación para cada electrodo. Por ejemplo, los datos en el archivo GLADOE2.DAT implica 161 electrodos, pero sólo las

15

se dan elevaciones en 26 puntos. El número máximo de puntos de datos topográficos que puede tener es 500. Para cada punto de datos, la ubicación horizontal y la elevación se introduce en el archivo de datos. Después de que el último punto de datos topográfica, se da el número del punto de datos topográfica donde se encuentra el primer electrodo. En el ejemplo anterior, el primer electrodo se encuentra a -40 metros, lo que corresponde a la segunda punto de datos topográfica. Tenga en cuenta que se requiere la elevación del primer electrodo. Si esto no se midió en el campo, se puede estimar que desde las elevaciones de los puntos de datos vecinos. Para conjuntos de datos con las mediciones de IP, así, por favor refiérase al  Apéndice F.

Un problema común encontrado en el uso de este programa son los errores e n el archivo de datos de entrada. Si el programa se detiene con un mensaje de error q ue podría ser debido a un error e n el archivo de datos de entrada. En primer lugar comprobar que los datos en el arc hivo están dispuestos según el formato descrito anteriormente. El programa intentará detectar errores comunes, tales cero o negativos valores de resistividad aparente, pero puede haber algunos errores nuevos exóticos que la versión actual de este programa no puede detectar. Otros posibles problemas se discuten en el "trampas en 2-D de inversión" de este manual.

importación de datos: Este programa también le permite ejecutar otros programas para convertir los datos dispuestos en un formato propietario en el formato requerido por este programa. Este programa se suministra normalmente por el fabricante del sistema de metro de la resistencia de múltiples electrodos que está utilizando.

Recopilar datos en formato RES3DINV - Por favor consulte el Apéndice P para los detalles.

correr J programa ACOBWIN - Este es un programa de utilidad que genera un número de archivo de soporte requerida por el programa RES2DINV. Debe ejecutar este programa después de instalar el programa RES2DINV. Sólo tiene que ejecutar el programa JACOBWIN vez.

Cambio segundo Uffer conducir - El programa utilizará automáticamente la unidad de disco duro con la mayor cantidad de espacio libre en disco como la unidad de memoria intermedia para almacenar los archivos temporales que se utiliza durante la inversión de un conjunto de datos. Esta opción le permite cambiar la unidad de memoria intermedia utilizada por el programa.

dieciséis

8

Edición de los datos

Esta opción le permite realizar algunos cam bios en los datos que ha leído en el archivo de datos de entrada se describe en la opción anterior. Se le permitirá eliminar puntos de datos malos, y para seleccionar una parte del conjunto de datos a invertir para grandes conjuntos de datos. Al seleccionar esta opci ón, se mostrará el siguiente submenú:

mi xterminate malos puntos de referencia - En esta opción, los valores de los datos de resistividad aparente se muestran en forma de perfiles para cada nivel de datos. Puede utilizar el ratón para eliminar cualquier punto de datos mal. El objetivo principal de esta opción es eliminar los puntos de datos que tienen valores de resistividad que son claramente erróneas. Tales puntos de datos malos podrían ser debido a la insuficiencia de los relés e n uno de los electrodos, un mal contacto de tierra del electrodo debido al suelo seco, o el cortocircuito a través de los cables debido a las condiciones del suelo muy húmedas. por lo general estos puntos de datos malos tienen valores de resistividad aparente que son obviamente demasiado grande o demasiado pequeño en comparación con los puntos de datos vecinos. La mejor manera de manejar este tipo de cosas malas es dejar caer ellos para que no influyen en el modelo obtenido. La Figura 4 muestra un ejemplo de un conjunto de datos con unos pocos puntos malos. Para eli minar un mal punto de datos, mover el cursor en forma de cruz con el ratón al punto de datos y haga clic en el botón izquierdo del ratón. El color del punto de datos debe cambiar de negro a púrpura. Si hace clic en el mismo punto de datos nuevo, no va a ser retirado del conjunto de datos. Para salir de esta opción, sólo tiene que pulsar el Q tecla o ha ga clic en la opción Salir.

S Plice gran conjunto de datos - Esta opción le permite elegir una sección del conjunto de datos completa (que es demasiado grande para ser procesada en una sola vez) para invertir. Después de elegir esta opción, se mostrará la distribución de los puntos de datos en un pseudosección. Puede seleccionar una sección del conjunto de datos para invertir mediante el uso de las teclas de flec ha. Todas las instrucciones sobre las teclas a utilizar se muestran en la pantalla. Los puntos de datos seleccionados serán marcadas por cruces o puntos morados, mientras que los puntos de datos restantes son de color negro. En la parte superior de la pantalla, los límites izquierdo y derecho de la sección seleccionados están marcados por líneas verticales de color amarillo. Si hay demasiados puntos de datos en los niveles inferiores, puede reducir el número de datos

17

puntos seleccionando puntos solamente pares o impares en un nivel. Para ello, mueva el marcador amarillo horizontal en el lado izquierdo de la pantalla hacia arriba o hacia abajo para el nivel que desea para diezmar con las teclas de flecha arriba y abajo. A continuación, pulse D para diezmar a los puntos de datos en ese nivel. El programa le permite leer en un archivo que incluye hasta 15.000 electrodos. Usando esta opción, puede seleccionar una sección del conjunto de datos (que cubre hasta un máximo de 500 electrodos y 8000 puntos de datos para un sistema con 256 megabytes de RAM) para invertir. Todo el conjunto de datos a continuación, se puede invertir mediante la inversión de las subsecciones sucesivas del conjunto de datos. Después de seleccionar el conjunto inciso datos que desea invertir, usted debe seleccionar la opción " mi xterminate malos puntos de datos" opción de buscar puntos de datos malos. En general, se recomienda que trate de invertir los datos de todo el conjunto en una sola vez. En la mayoría de los casos, esto se puede hacer fácilmente mediante la adición de más memoria RAM y el espacio libre en el disco duro a tu computador.

R pseudosección Everse - Esta opción da la vuelta al pseudosección horizontalmente de izquierda a derecha. Esto es útil cuando se tiene líneas paralelas de la encuesta, pero las encuestas se iniciaron desde distintos lados.

Figura 4: Ejemplo de un conjunto de datos con unos pocos puntos de d atos malos. Los datos se muestran utilizando el " mi xterminate opción malos puntos de datos".

do ambiar ubicación del primer electrodo - Esto le permite cambiar la ubicación del primer electrodo en la l ínea de levantamiento. Está destinado básicamente para fines de trazado, de manera que la superposición de las líneas de estudio tienen la misma x- ubicaciones para electrodos que coinciden.

18

mi archivo de datos DIT - Al seleccionar esta opción, el editor de texto (Notepad por defecto) se p ondrá en marcha. Para volver al programa RES2D INV, primero debe salir del programa editor de texto.

R ONU otro programa - Esta opción le permite ejecutar otro programa. Está destinado principalmente para los casos en que el programa no había sido instalado y por lo tanto no está en la lista de programas grabados en la opción de inicio de Windows y el programa. Escoger F ONT - En esta opción, puede elegir el tipo de fuente utilizada por el programa cuando se llega a las pseudosecciones. La fuente predeterminada es una fuente Arial. La otra opción es la fuente del sistema. Con la fuente predeterminada Arial el programa debe ser capaz de encajar todo en l a pantalla de la hora de elaborar los pseudosecciones, incluso con el modo 800x600 gráficos. La fuente del sistema por lo general tiene caracteres más audaz cuando podría ser mejor al hacer las copias impresas de las imágenes. Se recomienda que utilice el modo de gráficos de 1024x768 o mayor resolución cuando se utiliza la fuente del sistema.

Res t colores mineral - Esta opción restaura la combinación de colores por defecto utilizado por el programa.

19

9.

Cambio de los ajustes del programa

El programa cuenta con un conjunto de ajustes predefinidos para los factores de amortiguación y otras variables que, generalmente, dan resultados satisfactorios para la mayoría de los conjuntos de datos. Sin embargo, en algunas situaciones, es posible obtener mejores resultados mediante la modificación de los parámetros que controlan el proceso de inversión. Cuando se selecciona la opción " do Hänge Configuración" opción, se muestra la siguiente lista de opciones de menú.

Para cambiar un ajuste en particular, mover el cursor del ratón en la línea correspondiente y luego simplemente haga clic en él. Un mensaje aparecerá entonces lo que le da un poco de información acerca de los valores que puede utilizar, y el tipo de información que necesita para entrar.

Parámetros de amortiguamiento de inversión Las siguientes opciones modifican el uso del factor de amortiguamiento u en la ecuación (1). re factores amping - En esta opción, puede establecer el valor i nicial para el factor de amortiguamiento en la ecuación (1), así como el factor de amortiguación mínima. Si el conjunto de datos es muy ruidoso, se debe utilizar un factor de amortiguación relativamente grande (por ejemplo 0,3). Si el conjunto de datos es menos ruidoso, utilizar un factor de amortiguación inicial más pequeño (por ejemplo 0.1). La subrutina de inversión será generalmente reducir el factor de amortiguamiento en la ecuación (1) después de cada iteración. Sin embargo, un límite mínimo para el factor de amortiguamiento debe ajustarse para estabilizar el proceso de inversión. El valor mínimo por lo general debe establecer a alrededor de un quinto del valor del factor de amortiguación inicial. Chan gramo e de factor de amortiguación con la profundidad - Dado que la resolución del método resistividad disminuye exponencialmente con la profundidad, el factor de a mortiguamiento utiliza en el método de inversión de mínimos cuadrados es normalmente también aumentó con cada capa más profunda. Esto se hac e con el fin de estabilizar el proceso de inversión. Normalmente, el factor de amortiguación se incrementa en 1,05 veces, con cada capa más profunda, pero se puede cambiar. Utilice un valor grande si el modelo muestra oscilaciones no naturales en los valores de resistividad en las secciones inferiores. Esto ayudará a suprimir el

20

oscilaciones. También puede seleccionar la opción de permitir que el programa para determinar el valor para aumentar el factor de amortiguación con la profundidad de forma automática. Esto podría ser una buena opción si el espesor de las capas es mucho más delgado que los valores por defecto, por ejemplo, si ha reducido la distancia entre los electrodos unidad a la mitad en el archivo de datos con el fin de producir un modelo con bloques de modelo más pequeño.

O ptimise factor de amortiguamiento - Si se selecciona esta opción, el programa intentará encontrar el factor de amortiguación óptima u en la ecuación (1) que da el error RMS más baja en cada iteración. Al optimizar el factor de amortiguamiento, el número de iteraciones el programa requiere a converger puede reducirse significativamente. Sin embargo, el tiempo empleado por iteración se incrementará. Para pequeñas y conjuntos de datos de tamaño medio, esto puede reducir significativamente el tiempo total necesidad de ordenador para invertir el conjunto de datos. Para grandes conjuntos de datos con más de 1000 puntos de datos, el tiempo empleado en cada iteración se podría aumentar de manera significativa, ya que es necesaria para resolver los mínimos cuadrados ecuación más de una vez por iteración. Utilice esta opción con cuidado ya que el modelo con el menor error RMS no siempre es el “mejor”

modelo, en particular para los conjuntos de datos muy ruidosos (véase l a sección “Teoría”). En la práctica, la mejora en los resultados de la optimización del factor de amortiguamiento para la mayoría de los conjuntos de datos es más bien pequeño.

L IMIT gama de modelo de resistividad - Al s eleccionar esta opción, se mostrará el siguiente cuadro de diálogo.

Esta opción le permite limitar el rango de valores de resistividad que la subrutina de inversión le dará. En el ejemplo anterior, el límite superior para es 20 veces el valor medio modelo resistividad para la iteración anterior, mientras que el límite inferior es de 0,05 veces (es decir, 1/20 veces). El programa utiliza límites “blandos” que permiten que los valores reales modelo resistividad exceder los límites de a un cierto grado. Sin

embargo, esta opción evitará valores extremadamente grandes o pe queñas modelo de resistividad que son físicamente poco realista. El usuario también puede optar por utilizar el modelo de inversión obtenido en l a primera iteración como el modelo de referencia.

21

V ertical al cociente de filtro de planitud horizontal - Se puede seleccionar la relación entre el factor de amortiguamiento para el filtro de aplanamiento vertical ( fz) al filtro de planeidad horizontal ( fx). Por defecto, el mismo factor de amortiguación se utiliza para ambos. Sin embargo, si las anomalías en la pseudosección son alargadas verticalmente, puede forzar el programa para producir modelos que también son alargadas verticalmente mediante la selección de un valor más alto (por ejemplo 2.0) para la relación de la vertical para filtro de planeidad horizontal. Para anomalías que son alargadas horizontalmente, elija un valor más pequeño (por ejemplo 0,5).

Los parámetros de malla Este conjunto de opciones de cambiar la malla de diferencias finitas o elementos finitos utilizados en la subrutina modelado directo.

F Inite tamaño de la rejilla de malla - Se puede elegir una rejilla de malla utilizada por el programa de modelado adelante para tener 2 o 4 nodos entre electrodos adyacentes. Con 4 nodos por s eparación de los electrodos, los valores de resistividad aparente calculada sería más preciso (en particular para grandes contrastes de resistividad). Sin embargo, el tiempo en la computadora y la memoria requerida son correspondientemente mayor. Por defecto, el programa utilizará la opción 2 nodos si el conjunto de datos incluye más de 90 electrodos.

T SE método de elementos finitos - Este programa le permite utilizar ya sea la finitedifference o método de elementos finitos para calcular los valores de resistividad aparente. Por defecto, el programa utilizará el método de diferencias finitas, que es más rápido, si el conjunto de datos no contiene topografía. Si el conjunto de datos contiene la topografía, la opción por defecto es el método de elementos finitos.

Malla r efinement - Esta opción le permite utilizar una malla más fina (en la dirección vertical) para la diferencia finita o el método de elementos finitos. Los valores de resistividad aparente calculados por cualquiera de los métodos serán más precisas con una malla más fina, pero el tiempo de ordenador y memoria requerida será mayor. El uso de una malla más fina puede dar mejores resultados para los casos en subsuelo resistividad

22

Se espera que 1: contrastes de mayor que 20. Esto es particularmente útil en áreas en l as que una capa de baja resistividad se encuentra por debajo de una capa de alta resistividad.

 Ajustes de progreso de inversión El siguiente conjunto de opciones de control de la ruta de la subrutina inversión tiene durante la inversión de un conjunto de datos.

L ine búsqueda - La rutina de inversión determina el cambio en los parámetros del modelo mediante la resolución de la ecuación (1). Normalmente, el vector de cambio de parámetro re

dará lugar a un modelo con un error RMS inferior. En el caso de que los aumentos de error RMS, tiene dos opciones. Una opción es realizar una búsqueda de línea mediante la interpolación de cuarto grado para encontrar el tamaño óptimo de paso para el cambio en la resistividad de los bloques de modelo. El programa intentará reducir el error RMS, pero también puede ser atrapado en un mínimo local. La alternativa es ignorar el aumento en el error RMS, y esperamos que la próxima iteración dará lugar a un error RMS más pequeño. Esto le podría permitir a saltar de un mínimo local, pero también podría dar lugar un nuevo aumento en el error RMS. Una tercera alternativa es llevar a cabo una búsqueda de línea en cada iteración. Esto suele dar el tamaño óptimo de paso, sino que requerirá al menos un cálculo hacia adelante por i teración. En algunos casos, los cálculos de avance adicionales podrían valer la pena si se reduce el número de iteraciones necesarias para llevar el error RMS hasta un nivel aceptable. Para seleccionar un ajuste en particular, basta con hacer clic en el botón correspondiente. Este ajuste sólo afectará el proceso de inversión para la tercera iteración en adelante. Para las dos primeras iteraciones, donde los mayores cambios en el error RMS generalmente ocurre, el programa llevará siempre a cabo una búsqueda de línea para encontrar el tamaño de paso óptima para reducir aún más el error RMS.

PAG orcentaje cambio para la búsqueda de línea - El método de búsqueda de línea utilizada puede estimar el cambio esperado en el error RMS resistividad aparente. Si la esperada

23

cambio en el error RMS es demasiado pequeña, puede que no sea la pena continuar con la búsqueda de línea para determinar el tamaño óptimo de paso para el cambio de vector de parámetro de modelo. Normalmente se utiliza un valor entre 0,1 y 1,0%.

do límite onvergence - Esto establece el límite inferior para el cambio relativo en el error RMS entre 2 iteraciones. Por defecto, se utiliza un valor de 5%. En este programa, el cambio relativo en el error RMS, en lugar de un valor de RMS abs oluta, se utiliza para dar cabida a diferentes conjuntos de datos con diferentes grados de ruido presente.

límite RMS convergencia - Esto establece el error RMS porcentual en la inversión de los datos de resistividad aparente en el que el programa se detendrá después de que el modelo de productos tiene un error RMS inferior a este límite. Normalmente, un valor de entre 2% y 5% se debe utilizar, dependiendo de la calidad de los datos. norte umber de iteraciones - Esto permite al usuario fijar el número máximo de iteraciones para la rutina de inversión. Por defecto, el número máximo de iteraciones se establece en 5. Para la mayoría de los conjuntos de datos, esto es probablemente suficiente. Cuando la rutina de inversión alcanza este límite máximo, se le pedirá al usuario el número de iteraciones adicionales si desea continuar con el proceso de inversión. Por lo general, no es necesario el uso de más de 10 iteraciones.

METRO valores de resistividad Odel comprobar - El programa mostrará una advertencia si después de una iteración en la inversión del conjunto de datos, un valor de resistividad modelo se convierte en demasiado grandes (más de 20 veces el valor máximo resistividad aparente) o demasiado pequeña (menos de 1/20 el valor de resistividad aparente mínimo). Esta opción le permite desactivar la advertencia.

Datos / Selección de pantalla Esta sección tiene opciones de menor importancia para la pre-procesamiento de archivos de datos y visualización de las secciones durante la inversión.

24

do u límite toff para lecturas de resistencia de pozo - En mediciones transversal del pozo, ciertas combinaciones de los electrodos de corriente y potenciales podría resultar en valores de potencial extremadamente bajas. Esto podría resultar en lecturas con muy ba ja de señal a ruido. Para filtrar tales lecturas potencialmente ruidosos, el programa estima el potencial por una corriente de entrada de la unidad que se mide por cada configuración de electrodo usado en el conjunto de datos. Si el potencial estimado se sitúa por debajo de un cierto límite, en comparación con la combinación que da el máximo potencial, se rechaza la lectura. Por ejemplo, si se elige un valor de 0,001, esto significa configuraciones que darán lugar a valores de potencial que son miles de veces más pequeño que el potencial máximo medido será rechazada. Normalmente se utiliza un valor de entre 0,008 y 0,0007. El programa también le permite introducir las mediciones como valores de resistencia, y llevar a cabo la inversión con la resistencia (en lugar de los valores de resistividad aparente). En este caso, no se utiliza la opción de límite de corte. Para más detalles, por favor refiérase al Apéndice I.

O pción para los intervalos de contorno - Por defecto, el programa utilizará intervalos de curvas logarítmicas para los pseudosecciones y secciones modelo durante la visualización del resultado en la opción “Inversion” en el menú

principal. Esta suele ser la mejor opción para la mayoría de los conjuntos de datos. Sin embargo, puede optar por utilizar el los definidos por el usuario opciones de intervalos de curvas lineales o si lo desea.

Espectáculo pag seudosections durante la inversión - Usted puede optar por mostrar los pseudosecciones durante la inversión de datos, o simplemente mostrar los valores RMS modelo. Se recomienda que usted exhibe las secciones a menos que el controlador de Windows para la tarjeta gráfica es extremadamente lento. A largo plazo, se recomienda que descargue una nueva versión del controlador de Windows.

Guardar los parámetros de inversión - Esta opción guarda los parámetros de inversión en un RES2DINV.IVP archivo. Leer opción de parámetros de inversión -Este vuelve a leer los parámetros almacenados en los archivos RES2DINV.IVP y los utiliza en el programa. RES2DINV_NEW.IVP es un archivo de ejemplo que contiene más parámetros de inversión. Estos archivos también pueden ser utilizados por la o pción 'modo por lotes”. Por favor,

consulte el Apéndice G para el nuevo formato para el archivo de script 'modo por lotes'.

25

opciones de inversión 10. Esta opción le permite llevar a cabo la inversión del conjunto de datos que se ha leído en el uso de la " F ile" opción. También puede mostrar la disposición de los bloques utilizados por el modelo de inversión, así como para cambiar algunos de los parámetros que controlan el proceso de inversión.  Al seleccionar esta opción se mostrará el siguiente menú.

L este cuadrados inversión - Esta opción se iniciará la rutina de inversión mínimos cuadrados. Se le pedirá el nombre del archivo de datos de salida en el que almacenar los resultados, y los intervalos de contorno para los pseudosecciones si hubiera elegido la opción definida por el usuario para las curvas de nivel. Puede detener la inversión pulsando la tecla Q, y luego esperar por un corto tiempo.

Métodos de inversión y Configuración Este conjunto de opciones le permiten seleccionar el tipo de método de inversión regularizada a utilizar.

yo NCLUDE suavización de valores de resistividad modelo - La formulación de mínimos cuadrados se usa en la ecuación (1) se aplica una restricción de suavidad en el vector modelo de perturbación re solamente, y no directamente en los valores de resistividad modelo. En la mayoría de los casos, se producirá un modelo con una variación razonablemente suave de la resistividad

26

valores. En algunos casos, en particular para los conjuntos de datos muy ruidosos, mejores resultados pueden ser obtenidos mediante la aplicación de una restricción de la suavidad en los valores de resistividad modelo también. La ecuación de mínimos cuadrados resultante está dada por ( J T J + u F) d = J T gramo - u fr

(2)

dónde r es un vector que contiene el logaritmo de los valores de modelo de resistividad. Mientras que para los mismos factores de amortiguación esto generalmente producirá un modelo con un error resistividad RMS aparente más grande, esta modificación se asegurará de que el modelo resultante muestra una variación suave en los v alores de resistividad. El uso combinado yo Método nVersión - Esta opción está destinada para su uso en situaciones inusuales, donde los valores de sensibilidad de datos de los bloques de modelo se distorsionan significativamente por las grandes variaciones de resistividad. En algunas situaciones, como por ejemplo una encuesta sobre un cuerpo muy baja resistividad, las trayectorias de corriente podrían estar distorsionados de tal manera que las partes del subsuelo no están bien asignadas y tienen valores muy bajos de sensibilidad de datos en el modelo de inversión. Esto podría conducir a grandes distorsiones justo por debajo del cuerpo de baja resistividad. Mediante la combinación de la Marquardt (o regresión Ridge) y Occam (o suavidad con limitaciones) métodos de inversión, las distorsiones en algunos casos pueden ser reducidos. Esta opción se debe utilizar como último recurso si falla todo lo demás! Parece dar mejores resultados en la resolución de estructuras compactas donde el ancho y el espesor son ligeramente menor que la profundidad, tal como una cueva o mineral-cuerpo cuyo tamaño es sólo un poco menor que su profundidad. En teoría, es posible utilizar sólo la Marquardt (o regresión Ridge) método solo. Sin embargo, se ha encontrado que este método es inestable cuando los datos contienen ruido si se utiliza un modelo con un número relativamente grande de parámetros del modelo.

S electo de la inversión robusta - Esto le permite seleccionar el método de inversión sólida o en bloques. Se debe utilizar cuando se esperan límites precisos a estar presente. Por favor, consulte el Apéndice K para obtener más detalles. Elija l o garithm de resistividad aparente - Por defecto, el programa utilizará el logaritmo de los valores de resistividad aparente como el parámetro de datos al realizar la inversión. Para la mayoría de los casos, esto da los mejores resultados. En algunos casos, por ejemplo, con negativo o resistividad aparente cero, esto no es posible. Esta opción permite que el valor de la resistividad aparente por sí mismo para ser utilizado para este tipo de situaciones.

J cálculo matricial acobian - En este programa tiene tres opciones en el cálculo de la matriz jacobiana J en la ecuación (1). El método más rápido es utilizar el método cuasi-Newton para estimar la matriz jacobiana (Loke y Barker 1996a). Esto podría ser una buena alternativa para usar en el campo donde se utiliza con frecuencia un ordenador portátil más lento con la memoria más limitada. El método más exacto y más lenta es para volver a calcular la matriz jacobiana para todas las i teraciones. Esto requiere un equipo rápido con preferiblemente al menos 256 megabytes de RA M, y alrededor de 512 megabytes de espacio libre en disco duro para muy grandes conjuntos de datos.

27

Una tercera alternativa es para volver a calcular la matriz jacobiana para los primeros 2 sólo iteraciones, y usar el método de actualización cuasi-Newton para las iteraciones posteriores. Los mayores cambios en la matriz jacobiana por lo general se produce en las primeras pocas iteraciones. Por lo tanto, en muchos casos, un nuevo cálculo limitado de la matriz jacobiana da el mejor compromiso entre velocidad y precisión. Por defecto, el programa elegirá la opción de cálculo limitada por la matriz jacobiana. Sin embargo, la configuración por defecto pueden cambiarse por el usuario en esta opción de menú. La opción cuasi-Newton es útil para obtener un rápido vistazo a los resultados, sobre todo en el campo durante el curso de una encuesta utilizando un ordenador portátil. Para el modelo final, lo mejor es utilizar la opción de cálculo limitada en un ordenador de escritorio más rápido con más memoria RAM y el espacio libre en el disco duro.

Tipo de optimi s Método ación - Esta opción le permite elegir dos métodos diferentes para resolver la ecuación de mínimos cuadrados (1). Al seleccionar esta opción de menú, se mostrará el siguiente cuadro de diálogo.

Por defecto, el programa utiliza el método de mínimos cuadrados el 'estándar de Gauss-Newton', particularmente si el número de puntos de datos y / o células modelo es pequeño (menos de unos pocos miles), donde una solución exacta de la ecuación de mínimos cuadrados es calculado. Si el número de puntos de d atos y / o células modelo es grande (más de 2000), el tiempo necesario para resolver la ecuación de mínimos cuadrados podría ser la parte que consume más tiempo del proceso de inversión. Para reducir el tiempo de inversión, un método alternativo que calcula una solución aproximada de la ecuación de mínimos cuadrados utilizando el método de 'incompleto de Gauss-Newton' se puede utilizar. El usuario puede establecer la exactitud de la solución. Para la mayoría de los conjuntos de datos, una precisión de

28

aproximadamente 0,5% (es decir, un límite de convergencia de 0,005 en el cuadro de diálogo anterior) parece proporcionar una solución que es casi el mismo que el obtenido por el método de la 'Standard Gauss-Newton'. El establecimiento de una mayor precisión, por ejemplo 0,1%, en teoría, se dará un resultado que es aún más estrecha con el método de la 'Norma de Gauss-Newton', pero esto es a expensas de un tiempo de inversión más largo. Para grandes conjuntos de datos y modelos (generalmente más de 10.000 puntos de datos / Células modelo), una opción de “compresión de datos” se proporciona con el método de la 'Incompleto Gauss-Newton'. Esto puede reducir

significativamente la memoria del ordenador requerido.

Para los conjuntos de datos de la encuesta líneas muy largas (más de 2000 posiciones de los electrodos), se proporciona una opción de “inversión escasa” que puede reducir en gran medida el tiempo de ordenador y la memoria

requerida. Por ejemplo, la inversión de un conjunto de datos de la encuesta bajo el agua establecidos con casi 7.500 posiciones de los electrodos se redujo por un factor de 10 (a aproximadamente 101 minutos horas en un equipo de 3,2 GHz P4) con la opción de inversión escasa. Este método invierte el conjunto completo de datos y el modelo en una sola vez para producir un modelo continuo y sin problemas.

modelo Discretización Estas opciones le permiten modificar la forma en que el programa se subdivide el subsuelo en las células que se utiliza como modelo de inversión.

re bloques modelo ANTALLA - Esta opción mostrará la distribución de los bloques de modelo y puntos de datos. Los puntos de datos se representan gráficamente en la profundidad mediana de investigación (Edwards 1977) para la matriz utilizada. El programa tiene un límite máximo de 24 las capas del modelo. Si se hubiera elegido la opción 10% para el aumento en el espesor de capas, y la capa inferior está por encima de la fila inferior de puntos de datos, debe cambiar a la opción 25%.

do espesor hange de capas - En esta opción se puede elegir un modelo en el que el espesor de las capas de aumentar en un 10% o en un 25% con cada capa más profunda. Elija la opción 10% si sólo hay un pequeño número (8 o menos) de los niveles de datos. Si hay un gran número de niveles de datos dispersos, la opción del 25% podría ser mejor. Dentro de esta opción, también puede permitir que el número de bloques de modelo a superar el número de puntos de datos. El programa normalmente utiliza un modelo donde

29

la profundidad a la capa más profunda no exceda el pseudodepth máxima en el conjunto de datos . Para utilizar un modelo que abarca un rango d e profundidad más profunda, puede cambiar el factor para aumentar el rango de profundidad del modelo, por ejemplo. 1,0-1,30 para aumentar la gama de modelos de profundidad por 30%.

METRO profundidades odificar a capas - Esta opción le permite cambiar la profundidad de las capas utilizadas por el modelo de inversión. Se puede ajustar la profundidad de manera que algunos de los límites coinciden con profundidades conocidas a partir de datos de las perforaciones y otros. Utilice modelo extendido - Esta opción se extiende a las células del modelo a los bordes de la línea de levantamiento (Figura 5c). UN llow número de bloques modelo exceda de puntos de referencia - Por defecto, el programa tratará de disponer la posición y el tamaño de las células modelo de tal manera que no superen el número de puntos de datos. E sta es probablemente la mejor opción para grandes y medianas conjuntos de datos de tamaño recogidos con más de aproximadamente 50 electrodos, en particular cuando la distribución de los puntos de datos con las separaciones más grande es más escasa. En algunos casos, la anchura de las celdas de modelo en las capas inferiores puede ser más grande que los de las capas superiores. Sin embargo, para los conjuntos de datos más pequeños, que podría ser útil para relajar esta restricción, y permitir que el número de parámetros del modelo para superar el número de puntos de datos. Esto producirá un modelo en el que las células interiores en todas las capas tienen una anchura uniforme que es igual a la más pequeña distancia entre los electrodos (Figura 5b).

 Asegúrese de que los bloques de modelo tienen los mismos anchos de modelo -La muestra en la figura 5b tienen células modelo que son más anchos en los bordes. Esta opción se asegurará de que todas las células tienen la misma anchura que es igual a la distancia entre electrodos unidad. Reducir efecto de bloques laterales - Esta opción afecta el c álculo de los valores de matriz jacobiana para los bloques de modelo situados en los lados y el fondo de la sección de modelo. Normalmente, para un bloque situado en el lado, las contribuciones de todos los elementos de malla asociadas con el bloque modelo se añaden encima de la derecha hasta el borde de la malla. Esto da un mayor peso al bloque lateral en comparación con los bloques interiores. En algunos casos, particularmente cuando se utiliza la opción de inversión robusta, esto puede resultar en un valor inusualmente alta o baja resistividad para el bloque de lado. Esta opción deja fuera la contribución de los elementos de malla fuera de los límites de la línea de encuesta a los valores de matriz  jacobiana para los bloque s secundarios.

Cambiar el ancho de los bloques - Esta opción permite al usuario forzar el programa a usar células modelo que son más anchos que una separación entre electrodos unidad para todas las capas.

Utilice el refinamiento modelo - Por favor consulte el Apéndice P para los detalles.

30

Tipo de cruzada segundo orehole modelo - Hay dos opciones para el tipo de modelo se puede utilizar para la inversión de los datos a través del pozo. Puede optar por utilizar el modelo estándar, donde las dimensiones de los bloques modelo son los mismos que el espacio entre los electrodos en la superficie del suelo y en los pozos. Alternativamente, se puede utilizar un modelo con bloques más pequeños que el tamaño de los bloques es la mitad que el modelo estándar (véase la Figura 20 en el Apéndice I).

Figura 5. Diferentes métodos para subdividir el subsuelo en células rectangulares 2- D. Modelos obtenidos con (a) el algoritmo predeterminado, (b) por lo que el número de bloques de modelo para superar el número de puntos de datos, (c) un modelo que se extiende hasta los bordes de la línea de levantamiento y (d) utilizando los valores de sensibilidad para un modelo de tierra homogénea.

31

Opciones del modelo de sensibilidad

bloques de visualización s ensitivity - Esto mostrará un gráfico de la sensibilidad de los bloques utilizados en el modelo de inversión. El valor de sensibilidad es una medida de la c antidad de información sobre la resistividad de un bloque de modelo contenidos en el conjunto de datos de medición. Cuanto mayor sea el valor de la sensibilidad, más fiable es el valor del modelo de resistividad. En general, los bloques cerca de la superficie por lo general tienen valores de sensibilidad más alto s debido a que la función de sensibilidad tiene valores muy grandes cerca de los electrodos. Los bloques en los lados y el fondo también tienen altos valores de sensibilidad debido al tamaño mucho más grande de estos bloques que se extienden a los bordes de la diferencia finita o malla de elementos finitos. Si usted hubiera llevado a cabo una inversión del conjunto de datos antes de llamar a esta opción, el programa hará uso de la matriz jacobiana de la última iteración. De otra manera,

sens pantalla del subsuelo yo tividad - Esta muestra un gráfico de la sensibilidad de la subsuperficie para bloques de igual tamaño. Esto básicamente elimina el efecto de los cambios en el tamaño del bloque modelo de modo que muestra más claramente el cambio de la sensibilidad del subsuelo con la profundidad y ubicación.

Normalizar los valores de sensibilidad - Por defecto, los valores de sensibilidad calculados se normalizaron dividiendo con el valor promedio de sensibilidad. En esta opción, puede elegir no para normalizar los valores de sensibilidad. Generar bloques modelo - Esta opción genera las células modelo, haciendo uso de los valores de sensibilidad del modelo. Por favor, consulte el Apéndice L para una descripción de esta opción.

32

Opciones IP

Tipo de método de inversión de IP - Para i nvertir una resistividad combinado y conjunto de datos IP, puede optar por invertir tanto conjuntos de datos simultáneamente en cada iteración o para invertir la secuencia de datos IP después de la inversión de los datos de resistividad. IP factor de amortiguación - El factor de amortiguación utilizado en la i nversión de IP es generalmente más pequeño que el usado en la inversión del conjunto de datos de resistividad correspondiente. Si se utiliza un valor de 1,0, a continuación, los mismos factores de amortiguación serán utilizados tanto para la resistividad y las rutinas de inversión IP. Normalmente se utiliza un valor mucho más pequeño de aproximadamente 0,05 a 0,25. Alternativamente, se puede elegir la opción de dejar que el programa calcula automáticamente el factor de amortiguamiento.

segundo Modo atch - En esta opción, se pueden invertir una serie de conjuntos de datos de forma automática. El nombre de los archivos de datos de entrada, y otra información, se proporciona a través de un archivo. Se debe destacar que antes de invertir los conjuntos de datos, usted debe comprobar los datos de malos puntos de datos utilizando el mi opción de DIT datos discutidos anteriormente. Cuando se selecciona la subopción modo por lotes, primero se le pedirá el nombre del archivo de datos que contiene la información requerida. Un archivo de ejemplo proporcionado es RESIS.BTH que permite la inversión de 3 conjuntos de datos utilizando el modo por lotes. El formato de este archivo es la siguiente: -

Línea 1 - Número de archivos de datos a procesar. Líneas 2 a 16 - Este contiene la i nformación para el primer archivo de datos en el siguiente formato. Los valores correspondientes utilizados en el ejemplo RESIS.BTH para el primer archivo de datos también se enumeran a continuación para ayudar a aclarar las cosas.

Tipo de información

Valor en RESIS.BTH

Encabezamiento

ARCHIVO DE DATOS 1

Nombre del archivo de datos de entrada

RATHCRO.DAT RATHCRO.INV

Nombre del archivo de datos de salida  A partir factor de amortiguamiento factor de amortiguación mínimo opción de la línea de búsqueda que se utilizará

0.15 0.03 SIEMPRE

33 límite de convergencia

5.0

Cambio porcentual de búsqueda line

0.4

Número de iteraciones

6 1.00 10 4 NO

Vertical / relación de filtro planitud Horizontal El espesor de las capas del modelo aumento tamaño de la cuadrícula de diferencias finitas

Incluir alisado de modelo de resistividad

Reducir el número de puntos de referencia topográficos SÍ realizar la corrección topográfica

SÍ

Tipo de tendencia lineal para eliminar

DE EXTREMO A EXTREMO

Tipo de cálculo de matriz Jacobiana

RECÁLCULO LIMITED

Si no desea llevar a cabo las correcciones topográficas, basta con introducir “No” en la

línea correspondiente. Para eliminar un lineal de mínimos cuadrados o tendencia media, introduzca lineal o media en el archivo. Incluso si el conjunto de datos no contiene datos topográficos, aún tiene que entrar en el procesamiento de la información topográfica.  Aunque en este caso el procesamiento de la información topográfica será ignorado por el programa, el usuario debe entrarla en algo en las líneas apropiadas o de lo contrario el programa no puede leer en los parámetros para los siguientes archivos de datos en la secuencia correcta. Las otras opciones para el tipo de cálculo matriz jacobiana son “Quasi-Newton” y “Recalcular”. La información para los otros archivos de datos a procesar sigue el mismo formato. Por favor refiérase a la RESIS.BTH ejemplo de archivo.

Un nuevo formato de archivo de script 'modo por lotes también se describe en el Apéndice G q ue hace uso de los archivos de IVP que son creados por el programa RES2DINV. Es más fácil de usar ya que el usuario no tiene que ocuparse de los detalles sobre el formato del archivo que contiene los parámetros de inversión.

11.

modelado topográfico

Si hay relieve topográfico significativo a lo largo de la línea de la encuesta, el efecto de los efectos topográficos puede explicarse si se conocen las coordenadas horizontales y verticales de una serie de puntos a lo largo de la línea de encuesta. Cuando el programa lee en un archivo de datos con datos de topografía (por favor refiérase a la sección 7 para el formato de datos), se seleccionará automáticamente el método de elementos finitos que incorpora la topografía en la malla de modelado utilizado. En este caso, el modelado topográfico se llevará a cabo de forma automática por el programa cuando se invierte el conjunto de datos. Cuando se selecciona el menú “Opción Topografía”, se mostrará el siguiente submenú:

34

re topografía isplay - Esta opción simplemente traza la topografía, tal como en la Figura 6.

S Tipo electo de eliminación de la tendencia - Puede optar por eliminar la elevación media, una tendencia lineal por mínimos cuadrados, o una línea recta que une los primeros y últimos puntos de datos topografía. Si la s uperficie del suelo a lo largo de la línea de encuesta es generalmente plana a excepción de unos pocos puntos, utilice la opción para restar un valor constante de las alturas de los lugares electrodos. Si los puntos de medición están en una pendiente, a elegir entre unos mínimos cuadrados o una línea recta de extremo a extremo para la tendencia lineal. La figura 6 muestra un ejemplo con topografía donde se retiró la tendencia línea recta de extremo a extremo.

T ipo de modelado topográfico - Al seleccionar esta opción mostrará el siguiente cuadro de diálogo donde se puede elegir 3 diferentes métodos de modelado topográficos.

Distorsionada rejilla de elementos finitos con una distorsión uniforme: Este método, y los siguientes dos métodos, utilice una rejilla de elementos finitos deformada de tal manera que los nodos de la superficie de la malla se ajustan a la topografía real. Esto da resultados más precisos que el método más antiguo de la utilización de los factores de corrección para un modelo de tierra homogénea calculado con el método de elementos finitos (Fox et al. 1980) que pueden causar distorsiones en los casos en que grandes v ariaciones de resistividad se producen cerca de la superficie (Tong y Yang 1990, Loke 2000). En esta opción particular, los nodos por debajo de la superficie (y por tanto también las capas del modelo) se desplazan a la misma medida que los nodos de la superficie.

35

Distorsionada rejilla de elementos finitos con una distorsión amortiguado: En esta opción, los nodos del subsuelo se desplazan en menor medida en comparación con los nodos de la superficie, es decir, el efecto de la topografía es “amortiguado” con la profundidad. Esta opción es probablemente una opción razonable si la amplitud de la

curvatura de la topografía es menor que la profundidad de la capa del modelo más profundo. Un factor de amortiguación que controla el grado de amortiguación con la profundidad (figura 7) puede ser modificado por el usuario.

transformación SC con rejilla de elementos finitos distorsionada: Este método uti liza el método de transformación Schwarz-Christoffel para calcular la distorsión en las capas del subsuelo. Es probablemente el mejor método que produce un modelo de la sección más 'natural' que mira en los casos con una gran curvatura comparativamente topografía. En tales casos, el enfoque topografía amortiguada podría producir las capas del modelo inusualmente gruesas bajo áreas en las que las curvas de topografía hacia arriba. En algunos casos inusuales donde la topografía tiene picos agudos muy marcados y puntos de datos dispersos topografía,

la Schwarz-Christoffel método de transformación no podría funcionar. Si esto ocurre, añadir unos cuantos puntos de datos adicionales topografía cerca del pico.

Figura 6. Topografía eliminación de la tendencia lineal para el conjunto de datos montículo Rathcroghan.

36

Figura 7. Diferentes métodos para incorporar la topografía en un modelo de inversión 2-D. (A) Diagrama esquemático de un modelo de inversión 2-D típico sin topografía. malla A-elemento finito con cuatro nodos en la dirección horizontal entre electrodos adyacentes se utiliza normalmente. Las ca pas cercanas a la superficie también se subdividen verticalmente por varias líneas de malla. Los modelos con una rejilla deformada para que coincida con la topografía real donde (b) los nodos del subsuelo se desplazan verticalmente por la misma cantidad que los nodos d e la superficie, (c) el cambio en los nodos del subsuelo se reducen gradualmente con la profundidad o (d) reduce rápidamente con profundidad, y (e) el modelo obtenido con el método de transformación Schwartz-Christoffel inversa.

37

12.

Viendo los resultados de la inversión

En esta opción, se puede leer en un archivo de datos o el archivo de salida producida por la subruti na de inversión y mostrar los pseudosecciones r esistividad aparente medidos y calculados y la sección del modelo. Dentro de esta opción, se puede cambiar el intervalo de contorno usado para dibujar las secciones seudo y modelo, la escala vertical de las secciones, e incluyen la topografía e n la sección del modelo. También puede cambiar la combinación de colores que utiliza el programa. Si hubiera leído anteriormente en un archivo de datos de entrada o se utiliza la opción de inversión en el menú principal, el programa leerá primero el último archivo utilizado. Después de seleccionar esta opción en la barra del menú principal, el programa cambiará a la ventana para la visualización de los pseudosecciones y modelo de la sección. Para volver al menú principal, tendrá que salir de esta ventana. Esta ventana tiene la siguiente barra de menú superior.

Los " F ile opción”hará que aparezca el siguiente submenú.

R ead archivo de datos con resultados de la inversión - se puede leer en el archivo de datos de salida producidos por el subprograma de inversión o un archivo de datos de entrada.

S Ave de datos en formato XYZ - Puede guardar los resultados de la inversión en un archivo de disco con un formato usado por muchos programas de contorno, como Geosoft. Se trata de un programa separado vendido por GEOSOFT Inc., que le permite utilizar técnicas de contorno más sofisticados que los utilizados por el programa RES2DINV.

Guardar los datos en formato Surfer - Esto guarda los datos y las secciones de modelo en formato Surfer. Por favor, consulte el Apéndice N para los detalles. METRO Odel salida en formato LUND - Se trata de un formato de datos utilizado por el sistema de imágenes LUND producido por ABEM. Puede guardar los valores de resistividad modelo en los archivos .RHO y .rms utilizados por el programa MODSEC.

38

Salvar do Los valores ontour - Esta opción guarda los valores utilizados para dibujar los contornos de las secciones y los pseudo modelo.

Figura 8. La resistividad aparente observado pseudosección de la encuesta montículo Rathcroghan junto con la sección del modelo. La topografía se ha incorporado en la sección del modelo.  Al hacer clic en el “ re isplay opción secciones”en la barra de menú superior, aparecerá el siguiente submenú.

39

re ANTALLA datos y secciones modelo - Seleccione esta opción para trazar el modelo y la resistividad aparente pseudosecciones en la pantalla. yo NCLUDE topografía en pantalla modelo - Para resultados de la inversión en el modelado topográfico se llevaron a cabo antes, puede incorporar la topografía en la sección del modelo mediante el cálculo de que la verdadera profundidad y ubicación de los bloques de modelo debajo de la superficie. La figura 8 muestra un ejemplo de un modelo de la sección con la topografía. Los datos de campo para esta encuesta (así como el ejemplo de vertedero en la Figura 3) se midió usando un si stema de imágenes de resistividad Campus Geopulse fabricado por Campus Geophysical Instruments Ltd.

Para ambas opciones de visualización, se le pedirá que introduzca el número de iteración y las curvas de nivel. Después aparece la sección, puede mostrar los resultados para la siguiente iteración pulsando la tecla Av Pág. Del mismo modo, al pulsar la tecla Re Pág mostrará los resultados de la iteración anterior.

do resistividad hoose o la pantalla IP - Si el conjunto de datos tiene datos IP, así, esta opción le permite elegir para mostrar las secciones de resistividad o IP. También puede visualizar tanto la resistividad y las secciones modelo IP al mismo tiempo. T ipo de pantalla polo-dipolo - Esta opción sólo es relevante para los conjuntos de datos polo-dipolo con las mediciones con la “marcha atrás” conjuntos de polo-dipolo “hacia adelante” y (ver Apéndice A). Aquí se

puede elegir para mostrar los valores de resistividad aparente medi dos con la matriz en el“hacia adelante” o configuración “inversa”. Los " do ambiar opción de configuración de la pantalla”le permite cambiar los parámetros que controlan la forma en que se muestran los pseudosecciones de resistividad aparente y la sección del modelo. El siguiente submenú se mostrará al hacer clic en esta opción.

40

re combinación de colores EFAULT - Esto restablecerá la combinación de colores utilizado para colorear las sec ciones a un sistema predeterminado utilizado por el programa. do combinación de colores ustomise - Esto p ermite al usuario cambiar los colores utilizados para las secciones.

Leer personalizar el color s cheme - Puede cambiar el esquema de color a un conjunto de colores de personalizar que anteriormente se ha guardado en un archivo de disco.

norte úmero de secciones para mostrar - Puede optar por mostrar 1, 2 o 3 secciones en la pantalla. Por defecto, se muestran todos los 3 secciones (los pseudosecciones de resistividad aparente observados y calculados y la sección de modelo).

UN pParent diferencia de resistividad - Puede optar por mostrar ya sea la propia calculada pseudosección resistividad aparente, o la diferencia porcentual entre los logaritmos de los valores de resistividad aparente calculados y observados.

V ertical factor de escala de la pantalla - Esta opción le permite especificar la relación de la escala vertical de la escala horizontal, es decir, el factor de exageración vertical, en l a pantalla. Los valores convenientes para usar son 2,0, 1,5 y 1,0. Si introduce un valor

41

de 0.0, el programa utilizará un factor de escala predeterminado para que la pantalla se puede encajar en la pantalla de visualización.

MARIDO escala de trazado orizontal - Esto le permite cambiar la escala horizontal, en términos de número de píxeles por unidad de distancia entre los electrodos. Esta opción es útil cuando se desea representar los resultados de diferentes líneas de estudio con diferentes números de electrodos, con la misma escala.

S cómo ubicaciones de puntos de referencia - De manera predeterminada, la ubicación de los puntos de datos (y los centros de los bloques de modelo) se mostrará por puntos blancos en la pantalla a color pseudosección. Usted puede optar por no mostrar los puntos con esta opción.

R esistivity o de visualización conductividad - secciones Los pseudosecciones de resistividad y de modelo se visualizan normalmente. Sin embargo, puede mostrar las secciones usando valores de conductividad (es decir, el recíproco de la resistividad), que se utiliza más comúnmente en ciertos campos, t ales como hidrogeología.

Monitor pag valores seudodepth - por el pseudosecciones de resis tividad aparente, puede visualizar los valores pseudodepth en la escala vertical, o mostrar la s e paración de los electrodos y los factores “n” (en su caso). Color/ segundo visualización de la falta-Blanco - Por defecto, el programa mostrará l os pseudosecciones y modelo de la sección en color. Esta opción le permite visualizar las secciones en blanco y negro, posiblemente para hacer impresiones en impresoras láser normales. El programa ofrece dos opciones para la pantalla en blanco y negro. La primera opción es el uso de los colores blanco y negro puro junto con un método de tramado para proporcionar los diferentes tonos. Tales imágenes de la pantalla se pueden imprimir en casi cualquier impresora. La segunda opción es utilizar una escala de grises en diferentes tonos de gris, que van desde el blanco puro al negro puro, se utiliza. Es probable que tenga una impresora láser para hacer razonablemente buenas impresiones de las imágenes en escala de grises.

T ipo de pantalla modelo - En esta opción, se pueden visualizar los valores de resistividad en el modelo de la sección en forma de contornos (la opción por defecto), o en forma de bloques rectangulares que constituyen el modelo de inversión. La sección contorneada hace que sea más fácil de visualizar las estructuras geológicas, mientras que las características más finas pueden ser detectados más fácilmente en la sección de bloques rectangulares.

Cuando el " mi DIT datos”está seleccionada, la siguiente‘se visualiza el submenú estadísticas de error RMS’. Esta opción proporciona otro medio de eliminación de puntos de datos malos del conjunto de datos. D e hecho, para transversal pozo de sondeo o ciertos tipos de datos de las matrices no convencionales donde los valores de resistividad aparente no pueden ser convenientemente muestran en la forma de un pseudosección o perfiles, este es el método principal para eliminar los puntos de datos malos.

42

R estadísticas de error de MS - Esta subopción sólo muestra la distribución de la diferencia porcentual entre los logaritmos de los valores de resistividad aparente observadas y calculadas. Sólo se puede utilizar después de llevar a cabo una inversión de prueba del conjunto de datos. Para ahorrar algo de tiempo, se puede utilizar el método cuasi-Newton para llevar a cabo la inversión de ensayo (véase el apartado 10). Al seleccionar esta subopción, la distribución de los errores de RMS se mostrará en forma de un gráfico de barras. Se espera que los “malos” los puntos de datos tendrán relativamente grandes “errores”, por ejemplo por e ncima de 100 por ciento.

Para eliminar los puntos con grandes errores, basta con mover la línea de cursor verde para establecer el límite de error de corte para eliminar los “malos” los puntos de datos.

Los " PAG rint opción”en la barra de menú principal le permite guardar la imagen en pantalla como un archivo de gráficos BMP o PCX, o para hacer una i mpresión directa con una impresora utilizando el programa Windows Paintbrush. Al seleccionar esta opción, se mostrará el siguiente submenú.

Guardar la imagen de pantalla en forma de un archivo de gráficos BMP o PCX le permite editar la imagen con un programa de gráficos de mapa de bits, tales como el programa de Windows Paintbrush, antes de imprimirlo (ver sección 15).

13.

La configuración del programa

Por defecto, el Pincel de Windows se utiliza como el programa para manipular e imprimir el archivo de gráficos BMP o PCX, y el programa Bloc de notas se utiliza como el editor de texto. Puede cambiar estas opciones predeterminadas cambiando las líneas correspondientes en el archivo RES2DINV.INI que debe ser situado en el mismo subdirectorio que el programa RES2DINV.EXE .. Por ejemplo, u n archivo de configuración con el siguiente contenido: -

[programas] pintar = c: \ PSP \ PSP.EXE editor = c: \ gwd \ gwd.exe

43

[archivos] . Extension = * Entrada de lok = c: \ RES2DINV \ = salida c: \ RES2DINV \

utilizará el programa Paint Shop Pro gráficos (que se encuentra en el subdirectorio C: \ PSP) y el editor de texto de la DAS (en el subdirectorio C: \ DAS). También buscará archivos con la extensión LOK cuando se elige la opción de leer en un archivo de datos. El subdi rectorio por defecto con los archivos de datos es C: \ RES2DINV, y el mismo subdirectorio se utiliza para los archivos de salida.

14.

 Ayuda

Si hace clic en la opción “Ayuda”,  se mostrará el archivo de ayuda de Windows para el programa. Al hacer clic en la opción “Info” se mostrará el aviso de copyright.

15.

Impresión de los pseudosecciones

 Además de mostrar las pseudosecciones en la pantalla, se puede hacer una copia impresa de las secciones de modelo. Hay tres métodos para hacer una impresión. El método más simple es utilizar el RES2DINV.EXE para mostrar las secciones pseudosección y modelo que desee imprimir utilizando la opción “Mostrar resultados de la inversión” en el menú principal. Después de que el

programa se ha completado el dibujo pseudosección, haga clic en la opción Imprimir. Esto mostrará las opciones para guardar la pantalla como un archivo BMP, PCX o un archivo para abri r el programa de la brocha para imprimir el archivo. Si necesita imprimir la imagen de la pantalla junto con un poco de texto, que podría ser más conveniente para guardar la imagen de la pantalla como un archivo BMP o PCX. Puede insertar más tarde la imagen directamente en el documento si el procesador de textos, como Microsoft W ord para Windows o Lotus  Amipro o WordPerfect, le permite importar archivos de gráficos.

Otro método consiste en guardar la información en un archivo XYZ (véase la sección“Archivo”), que es un formato general con el apoyo de una serie de gráficos dedicados contorno programa / trazado. Con un programa de este tipo, se puede volver a dibujar las secciones modelo e imprimirlas. Esto es particularmente útil si desea combinar los resultados de diferentes líneas de estudio dentro de la misma área para hacer un pseudo 3-D parcela de las estructuras. Si desea incluir la topografía en los valores del modelo de guardados en el archivo XYZ, tendrá que ejecutar primero la opción “Incluir en la sección

topografía modelo” antes de guardar el archivo XYZ. En este caso, el archivo XYZ contendrá los valores del modelo sin la topografía, seguida de una sección con el modelo después de la incorporación de la topografía.

44

El tercer método es salvar la información del modelo en formato Surfer, y utilizar el surfista de programa de Windows para dibujar e imprimir las secciones. Por favor, consulte el Apéndice N para obtener más información sobre esta opción.

dieciséis

Otros programas útiles

Dos programas shareware que se han encontrado para ser útil son el programa de procesamiento de imágenes Pro taller de pintura, y el programa editor de texto DAS. Ambos s on programas basados en Windows. El programa Paint Shop Pro le permite llevar a cabo una serie de procesamiento más sofisticado para mejorar los gráficos de mapas de bits archivos BMP y PCX. Por ejemplo, puede convertir los archivos a un formato GIF, que ocupa mucho menos espacio en c omparación con los archivos BMP y PCX. El editor de texto DAS es un editor de nivel profesional que le permite realizar muchas operaciones que no son posibles con el programa Bloc de notas. Usted puede hacer uso de estos programas, en lugar del programa Bloc de notas por defecto Microsoft Paintbrush y, mediante la modificación de las líneas correspondientes en el archivo RES2DINV.INI. El archivo RES2DINV.IN2 da un ejemplo.

También hay una serie de programas de dibujo shareware finas y editores de texto que se puede descargar directamente de sitios de shareware de Internet. Además de Paint Shop Pr o, los programas de dibujo que he probado son Lview Pro, VuePrint y artes gráficas Taller para Windows. Todos estos programas le permiten cambiar los archivos BMP y PCX al formato GIF que puede reducir drásticamente el tamaño de los archivos gráficos. Puede descargar estos programas shareware de la mayoría de los sitios de shareware de Internet, tales como www.softseek.com.

El programa de dibujo en 2D surfista de Golden Software Inc. (www.golden.com) tiene una gran variedad de opciones útiles para la elaboración de secciones. El programa RES2DINV tiene una opción para exportar los datos al formato Surfer.

17

Errores en 2-D inversión

Si el programa todavía no da un modelo satisfactorio (suponiendo que los datos de entrada está en el formato correcto), que podría ser debido a uno de los siguientes problemas. Basura dentro basura fuera : Un problema común es malo puntos de datos. Comprobar si los valores inusualmente altos o bajos de resistividad aparente. Estos valores anómalos pueden ser causadas por errores de instrumentación durante el estudio de campo, un mal contacto del electrodo en tierra seca, arena o pedregoso, un cortocircuito de los electrodos debido a condiciones muy húmedas o con objetos metálicos (tales como cercas, tuberías, etc.) o errores como la colocación de electrodos a los conectores equivocados. Usted puede recoger fácilmente estos valores anómalos utilizando la opción "Editar datos" descrito anteriormente. Recuerde, un buen modelo requiere de buenos datos.

45

No unicidad: Es bien conocido en la inversión de los datos de resistividad de resonancia que es posible que los valores de resistividad aparente calculados calculados para ligeramente diferentes modelos a ser casi la misma. Lo mismo es válido para 2-D inversión de los datos de resistividad. Sin embargo, mientras que las características más finas en los diferentes modelos obtenidos para un solo conjunto de datos pueden ser diferentes, las características características principales principales deben ser similares. Además, tener en cuenta que el poder de resolución del método de resistividad resistividad en sí disminuye exponencialmente con la profundidad. Este programa le permite modificar algunos de los parámetros, tales como el factor de amortiguamiento y los filtros de planitud, que influirán en el proceso de inversión. Mediante el uso de diferentes parámetros de inversión, es probable que terminan con un poco diferentes modelos.

condiciones del terreno inusuales: Si la capa superior del suelo está muy sec o, puede que no sea capaz de bombear suficiente corriente en el suelo. En este caso, la técnica de l a resistividad en sí podría no ser adecuado y usted debe considerar métodos geofísicos alternativos (tales c omo EM). Otra situación similar es una capa superior muy conductora. En este caso, la diferencia de tensión medida por los electrodos de potencial podría ser demasiado pequeña (tal vez incluso por debajo del nivel de r uido de fondo). En esta situación, tratar de obtener más actual en el s uelo y comprobar la diferencia de tensión real medida. Si está utilizando la matriz dipolo-dipolo, vuelve a intentarlo con otra matriz como l a Wenner-Schlumberger que da los valores potenciales de mayor tamaño. Si esto n o funciona, intente utilizar otro método geofísico tales como refracción sísmica!

3-D geología: Cabe señalar que una limitación fundamental fundamental de este programa es la suposición de un modelo del subsuelo 2-D. Para las encuestas llevadas a cabo a través de la huelga de las estructuras geológicas alargadas, esta suposición es razonablemente exacta. Sin embargo, si hay variaciones significativas en la resistividad del subsuelo en una dirección perpendicular a la línea de encuesta (es decir, la geología es 3-D), esto podría causar distorsiones en las secciones inferiores del modelo obtenido. Las mediciones realizadas con separaciones de electrodos más grandes no sólo se ven afectados por las secciones más profundas del subsuelo, sino que también se ven afectados por las estructuras a una distancia horizontal más grande de la línea de encuesta. Este efecto es más pronunciado cuando la línea encuesta se coloca cerca de un contacto abrupto con la línea paralela al contacto.

Un programa de inversión de resistividad 3-D (Loke y Barker 1996b), RES3DINV, está ahora disponible. Para invertir un conjunto de datos en 3-D se establece con un número razonable de puntos de datos, se recomienda un microordenador microordenador Pentium rápido con al menos 64 megabytes de RAM! Con 64 megabytes de RAM, el programa de inversión de 3-D puede manejar encuestas rejillas de hasta 32 por 32 electrodos y un máximo de aproximadamente

11.000 puntos de datos. Un 3-D r esistividad libre libre y RES3DMOD programa de modelado directo IP también está disponible. Puede utilizar este programa para calcular el

46

valores de resistividad e IP aparentes para una variedad de estructuras 3-D para diferentes conjuntos de electrodos. Esto le dará una idea de la magnitud de los efectos 3-D. Hay una discusión sobre otras trampas interesantes en las notas libres “Tutorial” en la imagen eléctrica

que está disponible en el sitio web  www.geoelectrical.com.

18

 Algunos ejemplos de campo

 Algunos ejemplos de los conjuntos de datos de campo campo que se invirtieron con este programa son los siguientes: -

Grundfor Línea 2 - Esta encuesta fue llevada a cabo por el Departamento de Ciencias de la Tierra, Universidad de  Aarhus, Dinamarca Dinamarca para mapear la litología litología de los sedimentos sedimentos nearsurface no consolidada y acuíferos en la zona Grundfor (Christensen y Sorensen 1994). Las zonas de baja resistividad en el modelo de inversión (Figura 9) consisten en depósitos glaciares cuaternarios más arcillosos mientras que las áreas de alta resistividad consisten en depósitos más de arena. El modelo obtenido se confirmó por una serie de perforaciones a lo largo de esta línea. Se trata de un conjunto de datos moderadamente grande con 447 puntos de datos y 141 electrodos. En un Pentium 90 MHz que tomó cerca de 1.033 segundos (17,2 minutos) para invertir este conjunto de datos, mientras que en un Pentium 266 MHz II que tomó cerca de 158 segundos (2,6 minutos).

Odarslöv Dyke - El dique provoca una zona prominente alta resistividad aparente cerca de la mitad de la pseudosección (Dahlin 1996) en la parte superior de la figura 10. Esta es un conjunto de datos particularmente difícil de invertir como la anchura del dique de alta resistividad es menor que la profundidad a la sección inferior del dique. Así, la parte inferior del dique no está bien resuelto. Si el factor de amortiguación utilizado es demasiado pequeño, esto dará lugar a oscilaciones en los valores de resistividad en la sección del modelo. Para invertir este conjunto de datos, el peso del filtro de planeidad vertical fue fijado en el doble que la del filtro horizontal desde el dique tiene una dirección predominantemente predominantemente vertical. La opción sólida inversión de datos (ver Apéndice K) se utiliza para reducir el efecto de un par de puntos de datos que podrían ser más r uidoso en los resultados de la inversión. Puesto que hay un contraste agudo y distinto entre la resistividad del dique ígneas y las rocas sedimentarias circundantes, también se utilizó el robusto opción modelo de inversión (Apéndice K). En la sección de modelo, el dique se muestra como un cuerpo de alta resistividad prominente con lados verticales. Este conjunto de datos tiene 701 puntos de datos y 181 electrodos. Las encuestas Dyke Grundfor y Odarslöv se llevaron a cabo utilizando el sistema de imágenes de resistividad automática LUND producido por ABEM Instrument AB.

47

Figura 9: El pseudosección resistividad aparente aparente observada para la encuesta Grundfor Line 2 junto con un modelo de inversión.

Figura 10: El pseudosección resistividad aparente observada para la encuesta dique Odarslöv junto con un modelo de inversión.

48

Figura 11. Cartografía de agujeros en una c apa de arcilla, EE.UU. (a) pseudosección resistividad aparente para la encuesta para mapear los agujeros en la capa de arcilla inferior. (B) Modelo de Inversion y valores (c) sensibilidad de las células modelo utilizado por el programa de inversión.

Los agujeros en la capa de arcilla - USA

Esta encuesta se llevó a cabo con el propósito de agujeros de mapeo en una capa de arcilla que subyace 8 a 20 pies de arena limpia (pers Cromwell. Comm.). Los resultados de la encuesta de imagen eléctrica fueron confirmados posteriormente por las perforaciones. El pseudosección de una línea de este estudio se muestra en la Figura 11a. Los datos de la pseudosección se construyó utilizando datos de la encuesta líneas horizontales superpuestas. Una característica interesante de este estudio es que demuestra la naturaleza engañosa de la pseudosección, sobre todo para la matriz dipolo-dipolo. En el modelo de inversión, se detecta una alta anomalía resistividad por debajo de la ft. Marca 200, que es probablemente un agujero en la capa de arcilla inferior (Figura 11b). Est a característica cae en un área en el pseudosección donde hay una brecha a parente en los datos. Sin embargo,

49

cuerpo tiene valores de sensibilidad más altos (es decir, valores de resistividad modelo más fiable) que en las zonas adyacentes a la misma profundidad con más puntos de datos en la trama pseudosección (Figura 11c). Este fenómeno es debido a la forma de los contornos en la función de sensibilidad de la matriz dipolodipolo, donde las zonas con los valores más altos de sensibilidad están por debajo de los dipolos C1C2 y P1P2, y no en el punto trazado por debajo del centro de la matriz. Si las células modelo se colocan solamente en la ubicación de los puntos de datos, el cuerpo de alta resistividad no podrá contar a partir del modelo de inversión, y no se detectó una característica importante del subsuelo!

encuesta aguas subterráneas Bauchi - Un estudio de las aguas subterráneas mediante el Wenner alfa, beta y Wenner matrices se realizó en la zona de Nigeria Bauchi (Acworth 1987). En esta región, el agua subterránea se encuentra con frecuencia en la capa resistido por encima del lecho de roca cri stalina. La capa de degradado es más gruesa en las zonas con fracturas en la roca de fondo, y por lo tanto dichas fracturas son buenos objetivos para las aguas subterráneas. El conjunto de datos de matriz Wenner Gamma se muestra en la Figura 12 como un ejemplo. Para enfatizar el límite entre la capa de suelo y la roca de fondo, se usó la opción de inversión robusta (Apéndice K). El espesor de la capa inferior resistividad resistido es generalmente de aproximadamente 10 a 20 metros. Hay una zona d e baja resistividad vertical estrecha con una anchura de menos de 20 metros por debajo de la marca de 190 metros que es probablemente una zona de fractura en el lecho de roca. Un pozo de pozo que fue colocado en la marca de 175 metros que se encuentra en el borde de la zona de fractura tenía rendimientos que eran más bajos de lo esperado. La colocación del pozo se basa en gran medida de resistividad y EM de perfiles de mediciones utilizando el Geonics EM34-3, y a ntes de software resistividad inversión microordenador 2-D rápido era ampliamente disponibles.

Figura 12. encuesta array Bauchi Wenner Gamma. (un). pseudosección resistividad aparente. (B) El modelo de inversión con la topografía.

50

19

Solución de problemas

Una breve recopilación de posibles problemas y sus posibles causas y las soluciones se da en esta sección. Memoria insuficiente : Si el programa no puede encontrar suficiente memoria, una posible razón es que se ha reservado para otros programas. En este caso, usted tendrá que quitar los otros programas activos. En Windows, los programas activos generalmente se muestran en la barra de menú Inicio en la parte inferior de la pantalla. Si tiene menos de 512 megabytes d e RAM, se recomienda que usted compra algo de memoria RAM.

Programa se detiene cuando se trata de calcular las derivadas parciales: Los valores de derivadas parciales se almacenan en los archivos GRADWEN.10, GRADTWO.10, GRADDIP.10, GRADWEN.25, GRADTWO.25 y GRADDIP.25. En primer lugar comprobar que estos archivos están en el mismo subdirectorio de su disco duro que los archivos RES2DINV.EXE. Si están allí, un posible problema podría ser debido a sectores defectuosos en el disco duro. E s necesario para escanear el disco duro para comprobar si hay posibles sectores defectuosos sin marcar. Además, verifique que el disco duro tiene al menos 256 megabytes de espacio libre en el que el programa necesita para almacenar archivos temporales para grandes conjuntos de datos.

Programa se bloquea cuando se trata de invertir el conjunto de datos: Esto podría ser debido a los datos dispuestos en un formato incorrecto u otros errores en el conjunto de datos. Cuando el programa lee el archivo de datos, se comprobará si hay algunos errores comunes tales como cero o negativos valores de resistividad aparente. Sin embargo, puede haber algunos nuevos errores inesperados. Una comprobación preliminar que podría hacer es ver los datos utilizando los "Editar datos" opción que se describe en la sección 8. Si usted todavía no puede encontrar el error, por favor envíe una copia del conjunto de datos de la empresa que le vendió este programa (o la autor de este programa) por correo electrónico.

Programa no puede encontrar dongle: Si el programa no puede encontrar la clave de hardware, se ejecutará como un programa de demostración. Si usted hubiera fijado de forma segura el dispositivo de seguridad proporcionado al puerto de impresora y el programa todavía no puede detectarlo, por favor refiérase a la sección 3 para un método para comprobar la fuente del problema. El controlador del sistema de Keylok debe estar instalado para el dispositivo USB. Durante la instalación del controlador, el programa copiará el sistema necesario (SYS) los archivos que necesita Windows para detectar la llave USB. Por favor refiérase a la nota de instalación suministrado con el dispositivo USB para obtener más detalles.

51

Expresiones de gratitud

El Dr. Torleif Dahlin, de la Universidad de Lund, en Suecia proporcionó los conjuntos de datos interesantes y desafiantes Odarslöv Dyke, Gladkvarn y Knivsas. Los conjuntos de datos Grundfor y Cromer fueron amablemente proporcionados por el Dr. Niels B. Christensen de la Universidad de Aarhus en Dinamarca y el Dr. Torleif Dahlin. El conjunto de datos Rathcroghan fue proporcionado amablemente por el Dr. Kevin Barton y el Dr. Colin Brown a partir de datos recopilados por la Unidad de Geofísica Aplicada de la Universidad de Galway, Irlanda. El Dr. Andrew Binley de la Universidad de Lancaster proporcionado amablemente el conjunto de datos de campo crossborehole. El Sr. Ernst Neiderleithinger de Büro für Geophysik Lorenz (Alemania) suministró el conjunto de datos de relleno sa nitario que proporciona la motivación inicial para escribir este programa para manejar grandes conjuntos de datos 2-D. También agradezco al Prof. Douglas Oldenburg de la Universidad de la Columbia Británica en Canadá por señalar algunas de las posibles ventajas del uso de l a ecuación (2). Los datos de Bauchi fue proporcionado por el Dr. Ian Acworth de la Facultad de Ingeniería Civil y Ambiental, la ONU. de New South Wales, Australia. Muchas gracias a Richard Cromwell y Rory Retzlaff de Golder Assoc. (Seattle) para el ejemplo de encuesta para mapear los agujeros en una capa de arcill a. Bill Seaton de ATS, Virginia proporcionado amablemente el conjunto de datos BLUERIDGE.DAT. Dr. Fabrizio Cara del Instituto Nacional de Geofísica y Vulcanología, Italia suministró generosamente el conjunto de datos ITALY_CAVE.DAT. También me gustaría dar las gracias a Ron Kaufman de Technos Inc. para su permiso para utilizar el ejemplo de estudio de campo bajo el agua del lago Lelia. Me gustaría reconocer Michel Ritz del IRD, Francia amablemente para proporcionar el conjunto de datos de ejemplo BOREHOLE_TOPO.DAT. El conjunto de datos de la encuesta de electrodos flotantes a lo largo del río San Francisco en Brasil fue proporcionado amablemente por Rodrigo Machado. Por último, un agradecimiento especial a Rob Barker de la Universidad de Birmingham (Reino Unido) para el Birmingham de lapso de tiempo, la encuesta lecho de roca y la Yeguada Clifton (con Julian Scott) conjuntos de datos, y su asesoramiento a través de los años.

52

referencias  Acworth, RI, 1987. El desarrollo de los acuíferos en un basamento cristalino Trimestral tropical ambiente: diario

Ingenieria

Geología, 20, 265-272.

Barker, R. y Moore, J., 1998. La aplicación de lapso de tiempo eléctrica La tomografía en estudios de aguas subterráneas. El borde de ataque, 17, 1454-

1458.

Christensen y NB Sorensen KI, 1994. El uso integrado de electromagnético

métodos para investigaciones hidrogeológicas. Actas del Simposio sobre la aplicación de Geofísica de Ingeniería y los problemas ambientales, marzo de 1994, Boston, Massachusetts, 163176.

Claerbout, JF y Muir, F., 1973. modelado robusto con datos erráticos. Geofísica, 38, 826-844. Dahlin, T., 1996. 2D estudio de la resistividad para el medio ambiente y la ingeniería aplicaciones. Primer descanso, 14, 275-284.

Dahlin, T. y Zhou, B., 2004. Una comparación numérica de resistividad 2D de formación de imágenes usando diez conjuntos de electrodos. Prospección geofísica, 52, 379-

398.

Dahlin, T. y Zhou, B., 2005. mediciones matriz de gradiente múltiples para multicanal imágenes de resistividad 2D. Cerca de geofísica de superficie, en opinión.

deGroot-Hedlin, C. y Constable, S., 1990, inversión de Occam para generar

suave,

modelos bidimensionales formulario de datos magnetotelúricos. Geofísica, 55,

1613-1624. Edwards, LS, 1977, A pseudosección modificado para la resistividad y induced-

polarización. Geofísica, 42, 1020-1036.

Fox, RC, Hohmann, GW, Killpack, TJ y Rijo, L., 1980, topográfica efectos en la resistividad y encuestas de polarización inducida. Geofísica, 45,

75-93. Griffiths DH y Barker RD 1993, imágenes de resistividad de dos dimensiones y

modelado en áreas de geología compleja. Revista de Geofísica Aplicada, 29, 211-226. Griffiths, DH, Turnbull, J. y Olayinka, AI, 1990, Two-dimensional mapeo de resistividad con una matriz controlada por ordenador. Primer descanso, 8,

121-129.

Hallof, PG, 1990, Reconocimiento y resultados geofísicos detallados, Granito Montaña área del condado de Pershing, Nevada. en Fink, JB, McAlist er,

EO, Sternberg, BK, Wieduwilt, WG y Ward, SH (Eds), 1990, inducida por la polarización: Aplicaciones y historias de casos: Investigaciones en Geofísica No. 4, Soc. Expl. Geophys.

Keller GV y Frischknecht FC 1966. Métodos eléctricos en geofísica prospección. Pergamon Press Inc., Oxford.

53

Loke, MH, 1999. Foto de intervalo inversión de imágenes de resistividad. 5ª Reunión de la

 Actas sección europea de la Sociedad de Ingeniería Ambiental y geofísicos. Em1. Loke, MH, 2000. modelado topográfico de inversión de imágenes de resistividad.

62a Conferencia Eage y resúmenes extendidos exposición técnica, D-2. Loke, MH, 2001. Tutorial: 2-D y 3-D encuestas de imagen eléctricos. Geotomo Software, Malasia. Loke, MH y Barker, RD, 1996a, Rapid mínimos cuadrados inversión de aparente

pseudosecciones resistividad por un método cuasi-Newton. Prospección geofísica, 44, 131-152.

Loke, MH y Barker, RD, 1996b, técnicas prácticas para la resistividad 3D encuestas y datos de inversión. Prospección geofísica, 44, 499-523. Loke, MH y Dahlin, T., 2002. Una comparación de la de Gauss-Newton y métodos-Newton cuasi en inversión de imágenes de resistividad. Revista de Geofísica  Aplicada, 49, 149-162.

Loke, MH, Acworth, I. y Dahlin, T., 2003. Una comparación de suave y métodos de inversión en bloques en las encuestas de imagen 2D eléctrica. Exploración Geofísica, 34, 182-187.

Nelson, GD y Van Voorhis, 1973, Carta al editor en relación con el documento “Espectros resistividad Complejo de la mineralización de cobre porfídico”. Geofísica, 38, 984.

Niederleithinger, E., 1994, El uso de geomagnetismo de alta resolución, 2D -CCgeoeléctrica y la polarización inducida en investigaciones ambientales. Trabajo presentado en SAGEEP '94, Boston, EE.UU..

Olivar, A., de la cal, L. y Sharma, MM, 1990. Un conductividad grano acercarse a arcillosas areniscas. Geofísica, 55, 1347-1356. Robain, H., Albouy, Y., Dabas, M., Descloitres, M., Camerlynck, C., Menchler, P. y Tabbagh, A., 1997, La ubicación de los electrodos infinitos en estudios eléctricos de polo polos: consecuencias para la obtención de imágenes 2D. Revista de Geofísica Aplicada, 41, 313-333. Sasaki, Y., 1992, Resolución de la tomografía de resistivi dad infiere de numérica

simulación. Prospección geofísica, 40, 453-464. Scott, JBT, Barker, RD Peacock, S., 2000. refracción sísmica combinada y proyección de imagen eléctrica. Proc. 6ª Reunión de la Asociación Europea de Medio Ambiente e Ingeniería Geofísica, 3-7 de septiembre de 1997 Bochum, Alemania, EL05.

Seaton, WJ y Burbey, TJ, 2000. Caracterización de acuíferos en el azul Cresta provincia fisiográfica uso de perfiles de resistividad y geofísica de pozo: el análisis geológico. Journal of Environmental & Engineering Geofísica, 5, no. 3, 45-58. Silvester PP y Ferrari RL, 1990, elementos finitos para los ingenieros eléctricos (Segundo. Ed.). Prensa de la Universidad de Cambridge.

54

Slater, L., Binley, AM, Zaidman, MD y el Oeste, LJ, 1997, Investigación de mecanismos de flujo zona vadosa en tiza insaturado utilizando ERT crossborehole. Actas de la Reunión EEGS tercera sección europea, Aarhus, Dinamarca, 17-20.

Verano, JS, 1976, Principios de polarización inducida por geofísica exploración. Elsevier Scientific Publishing Company. Tong, L. y Yang, C., 1990, la incorporación de la topografía en dos inversión de resistividad dimensional. Geofísica, 55, 354-361. Turnbull, J., 1986. Un campo y la técnica de interpretación para levantamientos resistividad más de dos estructuras tridimensionales. Sin publicar. tesis doctoral, Un. de Birmingham, Reino Unido

Van Voorhis, GD Nelson, PH y Drake, TL, 1973, resistividad Complejo Los espectros de la mineralización de cobre porfídico. Geofísica, 38, 49-60. Blanco, RMS, Collins, S. y matrices Loke, P., 2003. resistividad e IP, optimizado para la recolección y la inversión de datos. Exploración Geofísica,

34, 229-232.

Witherly, KE y Vyselaar, J, 1990, Una historia del caso geofísica del álamo Lago depósito de cobre-molibdeno, del área de Houston, Columbia Británica. en Fink, JB, McAlister, EO, Sternberg, BK, Wieduwilt, WG y Ward, SH (Eds), 1990, inducida por la polarización: Aplicaciones y historias de casos: Investigaciones en Geofísica No. 4, Soc. Expl. Geophys.

55

 Apéndice A: encuestas dipolo-dipolo-dipolo y el polo El formato de los puntos de datos para un estudio dipolo-dipolo es ligeramente diferente de la utilizada para los Wenner y el polo polos. La disposición de los electrodos para la encuesta dipolo-dipolo se muestra en la siguiente figura.

Figura 13: (a) disposición de los electrodos para la matriz dipolo-dipolo. (B) Ejemplo de matriz dipolo-dipolo con no entero “ norte ”Valor. Disposición de los electrodos para (c) los “adelante” y (d) “inversa” arrays polo -dipolo. La distancia, un, entre los dos electrodos de corriente, C1 y C2, es la misma que entre los electrodos de potencial. La distancia entre el C1 corriente de electrodo y el electrodo de potencial P1 (el factor de separación dipolo) es generalmente un múltiplo entero, norte, de la distancia entre el par de electrodos actual o potencial. Inicialmente, el factor de separación dipolo ( norte) se establece en 1. A continuación, se aumentó a 2, 3 y así sucesivamente hasta un valor máximo de entre 4 y 6. Cuando se aumenta el factor  de separación dipolo, el potencial medido entre P1 y P2 disminuye rápidamente con el aumento de norte. Por esta razón, no es recomendable utilizar un valor de norte mayor que 6. Uso de valores más grandes de norte podría resultar en datos en lugar ruidoso. En algunos estudios, la distancia ( un) entre la corriente y el potencial

56

pares de electrodos está duplicado o triplicado (Edwards 1977) para aumentar la intensidad de la señal. Así, tanto el factor de separación dipolo ( norte) y la separación entre el par de electrodos de corriente ( un) ha de ser grabada. El uso del archivo BLOCKDIP.DAT como ejemplo, el formato de archivo de datos para la matriz dipolo-dipolo es la siguiente: archivo BLOCKDIP.DAT

comentarios

DOS bloques rectangulares | Título 1.0 3 356 0

| Unidad de separación de los electrodos

| número Array (dipolo-dipolo = 3) | Número de puntos de datos | x-localización dada en términos de primer electrodo | use 1 si se da la ubicación del punto medio

0

| 0 para no IP, utilice 1 si IP presente

0.0 1.0 2.0

1,0 1 30,1 1,0 1 30,1 1,0 1 30,1

. .

| Para cada punto de datos; introducir posición en x "un" | espaciamiento, "n" valor de factor, aparente | resistividad

| |

39,0 1,0 8 30,3 40,0 1,0 8 30,4 0,0,0,0,0

| 2 últimos puntos de datos

|

| Ponga unas Zero al final

El primer valor (en cada línea con los valores de los puntos de datos) da la posición en x del punto de datos. El segundo valor da la separación entre el potencial P1 par de electrodos y P2. Este valor es normalmente un múltiplo entero (de entre 1 a 5 veces) la dist ancia entre electrodos unidad. En algunos estudios de resistividad 2-D, la separación entre el par de electrodos de c orriente es fija. En este caso, este factor tendrá un valor constante (como en el archivo BLOCKDIP.DAT). El tercer valor es el factor de separación dipolo ( norte) con un valor máximo de 8. El cuarto valor es el valor de la resistividad aparente medido.

Los datos presentados BLUERIDGE.DAT es un conjunto ejemplo de datos de campo donde diferentes“ un ”Y“ norte ”Combinaciones (Seaton y Burbey 2000) con el fin de aumentar la profundidad de penetración sin reducir en gran medida la intensidad de la señal.

La matriz de polo-dipolo tiene un número de serie de 6. El formato de los puntos de datos e s similar a la de la matriz dipolo-dipolo. El archivo BLOCKPDI.DAT da un ejemplo de un conjunto de datos polo-dipolo. En comparación con la matriz dipolo-dipolo, la intensidad de la señal medida con el polo-dipolo matriz disminuye menos rápidamente cuando el factor de separación dipolo ( norte) está i ncrementado. A diferencia de las otras matrices comunes, el p olo-dipolo es una matriz no simétrica que da lugar a anomalías algo más complicado en el pseudosección. El punto central de esta matriz se define como el punto medio entre los electrodos C1 y P2. Este programa

57

se supone que la posición en x de electrodo C1 es menor que la del par de electrodos potencial P1 y P 2 (en el lado izquierdo en la figura 13c). Si la ubicación del primer electrodo se da en el archivo de datos de entrada, entonces se usa la localización del electrodo C1.

La matriz de polo-dipolo es una matriz asimétrica, y m ás de estructuras simétricas las anomalías de resistividad aparente en el pseudosección también son asimétricos. En algunas situaciones, la asimetría en los valores de resistividad aparente medidos podría influir en el m odelo obtenido después de la inversión. Un método para eliminar el efecto de esta asim etría es repetir las mediciones con los electrodos dispuestos de manera inversa (Figura 13d). Mediante la combinación de las mediciones con el arrays polo- dipolo “inversa” (Figuras 13C y 13D) “hacia adelante” y, cualquier sesgo en el modelo debido a la naturaleza asimétrica de esta matri z sería eliminado. El PDIPREV.DAT archivo es un conjunto de datos d e ejemplo de campos que combina las mediciones realizadas con el matrices poledipole “inversa” “adelante” y. Para la medición polo -dipolo inversa, la ubicación del primer electrodo está dado por la posición del electrodo P2. Un valor negativo para la

norte factor se utiliza para indicar que se utiliza la disposición polo-dipolo inversa. En este programa, el máximo permitido factor de separación dipolo norte es 32 mientras que la mayor separación entre el par de electrodos P1-P2 es de 16 veces la distancia entre electrodos unidad para la matriz dipolo-dipolo. Para la matriz de polo-dipolo, norte puede tener un valor máximo de 48, mientras que el espaciamiento P1-P2 puede ser de hasta 32 veces la distancia entre electrodos unidad. En la mayoría de los casos, el norte factor d e separación dipolo es un valor entero. Sin embargo, es posible que el norte factor a tener valores no enteros aunque todos los electrodos todavía tienen una constante de separación entre electrodos. Esto puede ocurrir cuando el espaciamiento un entre el p ar C1-C2 es dos veces (o más) la separación de los electrodos unidad. Para el ejemplo mostrado en la figura 13b, la distancia entre electrodos unidad para la línea de encuesta es 1 metro. los un espaciamiento tiene un valor de 2 metros, mientras que la distancia entre los electrodos C1 y P1 es d e 3 metros. En este caso, el norte factor tiene un valor de 1,5. El archivo de datos, DIPOLEN5.DAT, da un ejemplo de un os datos de dipolo-dipolo se establece con la no entero norte valores. Para los casos en que la parte fraccionaria de la norte factor tiene es una serie infinita en notación decimal, el norte valor se debe dar por lo menos a cuatro decimales. Por ejemplo, si el

norte valor es uno y un tercio, se le debe dar como 1,3333.

58

 Apéndice B: La matriz Wenner-Schlumberger En áreas con espacio abierto limitado por una línea de encuesta de largo, la matriz W enner convencional tiene la desventaja de que hay una gran reducción de la cobertura horizontal cuando el espaciamiento de los electrodos se incrementa a fin de lograr una profundidad más profunda de investigación. Por ejemplo, a fin de aumentar la profundidad de investigación por dos veces, el electrodo espaciado un tiene que ser aumentado a 2a (Figura 14). En este caso, la longitud total de la matriz se incrementa desde 3a a 6a.

 Al mismo tiempo, la anchura de la pseudosección disminuye por 3a con cada nivel de medida (Figura 15). En comparación, la anchura de la pseudosección disminuye por 2a para la matriz Wenner-Schlumberger.

Figura 14: Los pasos utilizados por el (a) Wenner y las matrices (b) WennerSchlumberger para aumentar la profundidad de investigación. El valor de resistividad aparente para la matriz Wenner-Schlumberger está dada por  ρ = pi n (n + 1) a R

dónde R es la resistencia medida, un es la separación entre los electrodos P1 y P2 y norte es la relación de las distancias entre el C1-P1 y los electrodos P1-P2. Esta matriz se convierte efectivam ente en la matriz de Schlumberger cuando el norte factor es mayor que 2. Por lo tanto, es realmente una combinación de las matrices de Wenner y Schlumberger adaptados para su uso para una disposición con una línea de electrodos con una separación constante (como normalmente se usa en formación de i mágenes eléctrico 2-D). Además de una mejor cobertura horizontal, la profundidad máxima de penetración de esta matriz es de aproximadamente 15% más grande que la matriz Wenner. Tenga en cuenta que la matriz normal de W enner es en realidad un c aso especial de la matriz Wenner-Schlumberger, donde la norte

factor es igual a 1. En la práctica, ya que el potencial entre los electrodos P1 y P2 disminuye cuando el factor de separación de electrodos norte se incrementa, el valor máximo de norte utilizado en los estudios de campo real es aproximadamente 8. Para aumentar la profundidad de investigación, la separación entre los electrodos P1-P2 se incrementa a 2a y las mediciones

59

se repiten para norte es igual a 1, 2, 3, 4, 5 y 6. A continuación, el espaciamiento P1-P2 se incrementa a 3a, y la misma secuencia de las mediciones se hace. El PIPESCHL.DAT archivo da un ejemplo de un conjunto de datos de campo establecidas usando esta matriz.

Figura 15: Disposición de los puntos de datos en los pseudosecciones para el (a) Wenner y la (b) arrays Wenner-Schlumberger. Mientras que la norte factor de por lo general tiene valores de número entero, el programa también puede aceptar los valores no enteros para este factor. El archivo de datos, WENSCHN5.DAT, es un ejemplo con valores no enteros de la norte factor.

Una ventaja de la matriz normal de Wenner es que tiene la mayor relación señal -tonoise entre las matrices convencionales. Para estar en el lado seguro, debe incluir todas las mediciones normales Wenner cuando se lleva a cabo una encuesta con la matriz de Wenner-Schlumberger. El pseudosección resultante podría parecer un poco desordenado con la superposición de niveles de medición de datos, pero RES2DINV puede manejar un conjunto de datos tales. En zonas muy ruidosas, las mediciones Wenner deberían ser los más fiables.

60

 Apéndice C: Las resistividades de fijación

En algunos casos, la resistividad verdadera de una sección de la sub-superficie podría ser co nocido, por ejemplo de las mediciones de resistividad de pozo. Este programa le permite fijar la resistividad de hasta 256 secciones del subsuelo. La forma de la sección que se va fijado sean rectangular o triangular. Los datos para fijar el modelo de valores de resistividad se introducen en el archivo de datos de entrada después de la sección con la información de la topografía. A modo de ejemplo, parte de la MODELFIX.DAT archivo de datos de ejemplo se enumeran a continuación.

0

| bandera Topografía 2 | Número si las regiones para fijar, poner 0 si no hay R | Tipo de la primera región, R para rectangulares

24,0.7 28,2.3 2.0 2.0

| X y Z coordenadas de la esquina superior izquierda del rectángulo

| X y Z coordenadas de la esquina inferior derecha del rectángulo | valor de la resistividad de la región rectangular

| Amortiguación peso factor T | Tipo de la segunda región, T para triangular 

30,0.0 30,3.0 45,3.0 10.0 2.0

| X y Z coordenadas del primer vértice del triángulo | X y Z coordenadas de segunda esquina del triángulo | Coordenadas de tercera esquina del triángulo | valor de la resistividad de la región triangular | Amortiguación peso factor de

El primer artículo después de la sección de datos de topografía es el número de regiones en las que la resistividad es que se determine. En el ejemplo anterior, 2 regiones se especifican. Si se da un valor de 0 (valor predeterminado), entonces no hay regiones donde la resistividad es especificado por el usuario. A continuación, se le da la forma de la región, R para rectangular o T para triangular. Si se especifica una región rectangular, a continuación, se dan las coordenadas X y Z de las esquinas superior izquierda e inferior derecha del rectángulo, como se muestra en la Figura 16.

Si se elige una región triangular, las coordenadas X y Z de los 3 vértices del triángulo se debe dar en un orden en sentido antihorario. Después se dan las coordenadas de la región que se fije, el siguiente elemento de datos es la resistividad de la región. Después de eso, se necesita el peso factor de amortiguación para la resistencia de la región. Este parámetro permite controlar el grado en el que la subrutina de inversión puede cambiar la resistividad de la región. Por lo general hay un cierto grado de incertidumbre en la resistividad de la región. mediciones de pozos sólo pueden dar la resistividad de una zona muy limitada cerca de la perforación. Por lo tanto, es recomendable que el programa debe permitir (dentro de límites) para cambiar la resistividad de la región. Si se utiliza un peso factor de amortiguación de 1,0,

61

El más grande se utiliza el peso factor de amortiguamiento, menor es el cambio que se permite en la resistividad de la región "fija". Normalmente, se utiliza un valor de aproximadamente 1,5 a 2,5. Si se utiliza un valor relativamente grande, por ejemplo 10.0, el cambio en la resistividad de la región sería muy pequeña durante el proceso de inversión. Tal un valor grande sólo debe utilizarse si la resistividad y la forma de la región se conoce con precisión.

Figura 16: Fijación de la resistividad de regiones rectangulares y triangulares del modelo de inversión.

62

D Apéndice: La matriz Ecuatorial dipolo-dipolo La disposición de los electrodos para esta matriz s e muestra en la siguiente figura.

Figura 17: La disposición de los electrodos para la matriz dipolo-dipolo ecuatorial.

El valor de resistividad aparente para esta matriz está dada por  ρ = 2 pi como R / (s - a) donde

s = (a 2 + segundo 2 ) 0,5

y R es la resistencia medida. Esta matriz se diferencia de las otras mat rices, donde los electrodos están dispuestos en una sola línea. Cada lectura se caracteriza por dos separaciones. La primera separación, a, es la distancia entre los electrodos de corriente par C 1-C1 y el par potencial P1-P2. El segundo espaciamiento, b, es la separación entre los electrodos C1 y C2. Este programa tiene una restricción de que la separación entre el par P1-P2 debe ser la mism a que entre el par C1-C2. Un ejemplo de un archivo de datos para esta matriz es la FAULTEQU.DAT archivo de datos. La parte inicial del archivo se muestra a continuación con los comentarios: -

archivo FAULTEQU.DAT

Ecuatorial matriz dipolo-dipolo 1.00 8 1.50 285 1 0 0.50 1.00 2.74 1.50 1.00 2.74

comentarios

| Título | El espaciamiento de los elec trodos unidad

| número matriz 8 para esta matriz

| La separación b entre el par C1-C2 | El número de puntos de datos | Bandera de punto medio de la matriz | No hay datos IP | x-ubicación, la distancia A, app. resis. | Repita este procedimiento para todos los demás puntos

Tenga en cuenta que para esta matriz, es necesario incluir la separación b entre el par electrodos C1C2 después del número de array.

63

 Apéndice E: estudios submarinos Las posibles situaciones para estudios submarinos se muestran en la siguiente figura.

Figura 18: tres situaciones posibles para estudios submarinos, y una encuesta con electrodos flotantes.

En la primera situación, hay una capa de agua con un espesor constante sobre una superficie del sedimento plana (es decir, no hay topografía). Un ejemplo de tal conjunto de datos se da en el WATER.DAT archivo. La sección inferior de este expediente de observaciones es la siguiente: -

archivo WATER.DAT

39.00 40.00 41.00

1,00 16 1,00 16 1,00 16

51.36 52.09 54.29

comentarios | Los 3 últimos puntos de datos | |

0

| 0 para no topografía

0

| 0 para no hay regiones fijas

1

| 1 para indicar agua capa presente

50.00, -100.00, 200.00

| La resistividad de la capa de agua, izqu ierda y derecha | límites de la capa de agua

4.00

| Espesor de la capa de agua

64

1

| Superficie factor geométrico utilizado para calcular | resistividades aparentes

0

|

Mientras que los límites izquierdo y derecho de la capa de agua se incluyen en el formato de datos, esto es en la actualidad no se utiliza por el programa. Se incluye para su uso futuro para las encuestas con algunos de los electrodos por encima de la capa de agua, por ejemplo, una encuesta a través de un río. Tenga en cuenta que una bandera se incluye para indicar que el factor geométrico normal (tal como se utiliza para un recorrido de superficie) se utiliza para calcular los valores de resistividad aparente. Si se utiliza un valor de 0 para este indicador, el factor geométrico exacta (que incluye el espesor de la capa de agua) debe ser utilizado para calcular los valores de resistividad aparente. Para los estudios de campo reales, se espera que se utilizará el factor geométrico superficie normal.

Si el río o el mar de fondo tiene la topografía, la información topografía es introducida inmediatamente después de la serie de datos como en una encuesta de la superficie terrestre normal. La sección inferior de los archivos es la siguiente: archivo de levantamiento bajo el agua wth topografía

39.00 40.00 41.00

1,00 16 1,00 16 1,00 16

comentarios

51.36 52.09 54.29

2

| Los 3 últimos puntos de datos | |

| 2 para la topografía superficial con x | distancias

7

| Número de puntos de datos de topografía

0.0, 3.0

| El x-ubicación y elevación de cada

10.0, 4.0 15.0, 4.5 25.0, 5.0 30.0, 5.1 35.0, 4.5 45.0, 4.6

| punto de datos topografía

1

| Para indicar 1 st electrodo está en 1 st punto topo

0

| 0 para no hay regiones fijas

1

| 1 para indicar agua capa presente

50.00, -100.00, 200.00

| La resistividad de la capa de agua, izquierda y derecha | límites de la

| | | | |

capa de agua

4.00

| La elevación de la superficie del agua con respecto a | datum

1

local | 1 para indicar factor geométrico superficie normal | usado

El programa calculará automáticamente el espesor de la capa de agua de la elevación de la mar o río abajo y la elevación de la superficie del agua. Tenga en cuenta que este formato permite también la tercera situación posible se muestra en la Figura 18, donde algunos de los electrodos están por encima de la superficie del agua (es decir, en

sesenta y cinco

tierra seca normal) y algunos de los subacuático electrodos. En el ejemplo anterior, los electrodos de entre 0 y 10 metros están por encima de la superficie del agua. Una situación común donde esto ocurre es una e ncuesta que cruza un río. En este caso, el programa as umirá que un electrodo con una elevación que es menor que la elevación de la superficie del agua está bajo el agua. Para el caso en q ue algunos de los electrodos están por encima de la superficie del agua, el factor geométrico superficie debe ser utilizado para el cálculo de los valores de resistividad aparente para todos los puntos de datos.

El modelado topográfico se lleva a cabo automáticamente por el programa cuando se invierte el conjunto de datos. El programa se ajustará automáticamente la superficie de la rejilla de elementos finitos utilizado para modelar el subsuelo para que coincida con la superficie de los sedimentos por debajo de la capa de agua.

Hay una cuarta situación es posible con un cuerpo cerrado de agua que hasta el momento no se admite. Una situación posible es un río o lago con diques. En este caso, algunos puntos con elevaciones por debajo del nivel de la superficie del agua no están bajo el agua. Soporte probablemente se añadirán en el futuro cuando se encuentran con mayor frecuencia este tipo de situaciones.

Muchos estudios submarinos se llevan a cabo arrastrando un cable o streamer detrás de un barco. En este caso, las ubicaciones de los electrodos es improbable que se igualmente espaciados. Esta situación se maneja mejor mediante la introducción de los datos en el formato general de matriz que se describe en el Apéndice L. El LAKELELIA.DAT archivo de datos es un estudio de campo interesante por Technos Inc. con datos en el formato de matriz general (Anexo L) y topografía submarina . Partes del archivo se destacan a continuación para ilustrar el formato de datos utilizado.

archivo LAKELELIA.DAT

comentarios

LAGO LELIA SUBMARINA ENCUESTA

| Título

1.50 11 1

| Unidad de separación de los electrodos

| Formato general matriz

| Nota Wenner tipo sub-array

Tipo de medida (. 0 = aplicación resi stivity.1 = resistencia) | cabecera 0 | De datos está en la resistividad aparente

291 2 0

| Número de puntos de datos

| Superficie x-distancias

| Sin IP

4 82.37,0.0 91.37,0.0 85.37,0.0 88.37,0.0 43,137 4 89.55,0.0 98.55,0.0 92.55,0.0 95.55,0.0 36,599

| Primeros cuatro puntos de datos

4 96.74,0.0 105.74,0.0 99.74,0.0 102.74,0.0 36,458

| elevación es 0,0 aquí

| Nota del electrodo

4 103.92,0.0 112.92,0.0 106.92,0.0 109.92,0.0 37,019 | porque se da más adelante

. .

| Otros puntos de datos | Otros puntos de datos

4 1050.25,0.0 1077.25,0.0 1059.25,0.0 1068.25,0.0 93,072 | Datos de la última dos

4 1056.35,0.0 1083.35,0.0 1065.35,0.0 1074.35,0.0 95,383 | Topografía puntos en | cabecera especial para indicar elevaciones lista separada

66

| se dan por separado del principal | sección de datos

2 38

| 2 para indicar la superficie x-distancias

26.12 999.09 30.48 999.09

| Los dos primeros topografía

. . 1.127,76 999,09 1

| otro topografía

| Número de puntos de datos de topografía | puntos de datos | puntos de datos | Última topografía punto de datos | El primer electrodo coincide con el primer | punto de datos topografía

0

| Poner a cero aquí en lugar de la habitual | topografía de la superficie

0

| No hay regiones fijas

1

| 1 para indicar agua capa presente

58.00, -2000.00,2000.00

| La resistividad de la capa de agua, a la izquierda y | limita el derecho de capa de agua

1000 1

| La elevación d e la superficie del agua | Superficie factor geométrico utilizado para | calcular resistividades aparentes

0,0,0

| Terminar con unos cuantos ceros

La Figura 19 muestra el modelo pseudosección resistividad y la inversión aparente para este conjunto de datos. Los sistemas de topografía móviles se caracterizan por producir líneas de estudio que son muy larga en comparación con la profundidad de investigación del sistema de red de cable utilizado.

Encuestas con electrodos flotantes Otra posible disposición para una encuesta utilizando un barco arrastrando un cable es donde los electrodos están flotando en la superficie del agua. La principal desventaja de esta disposición es que una gran parte de la corriente fluye dentro de la capa de agua, en comparación con la porción de la corriente que fluye en el material por debajo de la parte inferior de agua. La principal ventaja es que es más fácil y rápido de realizar y evita problemas tales como el c able que se está engancharse los obstáculos en el fondo del río. Podría ser un método eficaz en zonas poco profundas, donde el efecto de apantallamiento debido a la capa de agua n o es demasiado grande. Una posible manera de interpretar un conjunto tal de datos es tratarla como un conjunto de datos de superficie estándar, donde la capa de agua se incluye como parte del modelo. Esto probablemente da resultados aceptables donde hay un fuerte contraste entre la capa de agua y la capa superior de los materiales de sub-parte inferior, y donde los datos tiene un bajo nivel de ruido. Sin embargo, si se midió la profundidad de la parte inferior de agua durante el curso de la encuesta (por lo general con una ecosonda), y la resistividad del agua se midió también, el efecto de la capa de agua se pudo calcular e incluido en el modelo. El WATER_FLOAT.DAT archivo es un conjunto de datos sintético que da un ejemplo del formato de datos. La parte inicial del archivo es idéntica a la utilizada para el normal el efecto de la capa de agua se pudo calcular e incluido en el modelo. El WATER_FLOAT.DAT archivo es un conjunto de datos sintético que da un ejemplo del formato de datos. La parte inicial del archivo es idéntica a la utilizada para el normal el efecto de la capa de agua se pudo calcular e incluido en el modelo. El WATER_FLOAT.DAT archivo es un conjunto de datos sintético que da un ejemplo del formato de datos. La parte inicial del archivo es idéntica a la utilizada para el normal

67

encuestas de superficie o encuestas con electrodos en el fondo del agua. La diferencia está en la sección final del archivo donde se dan los parámetros de la capa de agua. archivo WATER_FLOAT.DAT

comentarios

Encuesta sobre la capa de agua y electrodos flotantes | Título

2.0 | Unidad de separación de los electrodos 11 | Formato general matriz 0 | Sin tipo de sub-array Tipo de medición (0 = app.resistivity, 1 = resistencia) | cabecera 0 | De datos está en la resistividad aparente | Número de puntos de datos 942 | distancias horizontales verdaderos 1 0 | Sin IP | Primeros tres puntos de datos 4 2.0,0.0 0.0,0.0 4.0,0.0 6.0,0.0 1,86738 4 4.0,0.0 2.0,0.0 6.0,0.0 8.0,0.0 1,86843 | superficie Nota agua se supone que es 4 6.0,0.0 4.0,0.0 8.0,0.0 10.0,0.0 1,87070 | a 0.0 metros . | Otros puntos de datos . |

4 46.0,0.0 34.0,0.0 94.0,0.0 106.0,0.0 14,91776 | Los dos últimos puntos de datos

4 48.0,0.0 36.0,0.0 96.0,0.0 108.0,0.0 15,51526 | 0 | No hay información topografía

0 3

| No hay regiones fijas

la resistividad del agua

| Encabezamiento

1.5

| La resistividad del agua en ohm.m

la elevación del agua

| Encabezamiento

0.0 18

| La elevación de la superficie del agua

| 3 para indicar encuesta electrodos flotantes

| Número de puntos de profundidad del fondo del agua

0 -1.5 10 -1.5 20 -1.5 26 -1.5 28 -2.0 30 -2.5

| x-ubicación y la elevación del fondo del agua por debajo de la superficie del agua

. .

| Otras líneas de datos de elevación inferior

100 -1.5 108 -1.5

|

1 0,0,0

| Indica primer electrodo está e n primera punto de profundidad

| Segundo punto de profundidad

| | | | | | El último punto de profundidad | Unos pocos ceros para terminar el archivo de datos

Mientras que el ejemplo anterior muestra los datos en el formato de matriz general (Apéndice

L), el programa también puede aceptar datos en el formato basado en el índice. Por ejemplo, el conjunto de datos anterior es para la matriz dipolo-dipolo. Un formato basado índice (ver Apéndice A) dará a la ubicación del punto de datos, la 'a' longitud del dipolo y el factor de separación dipolo 'n'.

68

69

Hay dos métodos que se pueden utilizar para invertir los datos recogidos con electrodos flotantes. El primer método es similar a la utilizada para las encuestas con electrodos en la parte inferior donde la resistividad del agua es fijo y el modelo se extiende por debajo de la parte inferior capa de agua. El segundo método incorpora la capa de agua en el modelo de inversión. Los dos métodos pueden ser seleccionados a través del 'Inversion - métodos de inversión y Configuración - electrodos flotantes método de encuesta inversión' opción.

 Al hacer clic en esta opción, se mostrará el siguiente cuadro de diálogo.

70

El primer método (Uso capa de agua fijo) es más adecuado cuando el espesor máximo de la capa de agua es pequeño (menos de 20%) en comparación con la profundidad máxima de la investigación de la configuración de la encuesta. La Figura 20 muestra la inversión del conjunto de datos WATER_FLOAT.DAT usando este método. Tenga en cuenta que una rejilla de elementos finitos distorsionada se utiliza para modelar la capa de agua. El espesor máximo de la capa de agua en este ejemplo es de 2,5 metros, que es pequeño en comparación con la profundidad máxima de la investigación en alrededor de 15 metros para el conjunto de datos.

Figura 20. Resultados de la inversión de los datos WATER_FLOAT.DAT fijaron usando el primer método de inversión, donde la capa de agua resistividad (representado por la capa superior azul en la sección del modelo) es fijo y el modelo se extiende desde por debajo del fondo del agua.

En los casos en que el espesor de la capa de agua es significativo, un segundo método de inversión (Incorporar la capa de agua en el modelo) es más adecuado. La sección de la capa de agua se incluye en el modelo de inversión. Este método tiene dos posiciones de inversión. En la mayoría de los casos, la resistividad de la capa de agua no varía significativamente a lo largo de la línea de encuesta o con la profundidad. Por lo tanto normalmente se selecciona la opción 'minimizar la variación de la resistividad del agua'. Sin embargo, para acomodar los casos en los que hay v ariaciones significativas en la resistividad del agua (como en una zona de mezcla entre el agua dulce y solución salina), la opción 'Permitir capa de agua resistividad varíe libremente' se puede utilizar.

El material del subsuelo por debajo de la capa de agua co n frecuencia tiene una resistividad mucho mayor que el agua. La opción 'Seleccionar un cambio brusco a través del agua contorno inferior' puede ser utilizado de manera que se permite un fuerte contraste entre la capa de agua y el subsuelo. Esta opción utiliza el método descrito en 'Apéndice T: La incorporación de límites de c apas conocidos' para crear el cambio brusco en el modelo de resistividad. El programa ajusta automáticamente el espesor de las capas del modelo y divide las células modelo en los que pertenecen

71

a la capa de agua y el subsuelo. La figura 21 muestra la disposición de las cél ulas modelo para una de las líneas de levantamiento a lo largo del río San Francisco, divisor de Bahía (BA) y Pernambuco (PE) establece, barrios de la Santa Maria da Boa Vista ciudad (PE), al noreste de Brasil (pers Rodrigo Machado. Comm.). En esta encuesta, un tipo dipolo-dipolo de configuración se utilizó con 5 dipolos del medidor. Sin embargo, las lecturas se tomaron cada metro que da un conjunto de datos con una separación e ntre electrodos unidad nominal de 1 metro. Para reducir el número de células modelo, y tomando en consideración que la resolución de dicha encuesta es muy poco probable que sea mejor que la mitad de la l ongitud del dipolo, la anchura de las células del modelo se fija en 3 veces la separación entre electrodos unidad, es decir, básicamente 3 metros para casi toda la sección. Tenga en cuenta también que la profundidad del agua en algunos lugares es mayor que la profundidad mediana máxima de investigación (como se indica por los marcadores de puntos de datos). Sin embargo, el programa establece la profundidad de las capas de tal manera que hay al menos una capa del modelo por debajo de la parte inferior de agua (representada por la línea azul oscuro). Las células por de bajo de la capa de agua se marca como de color azul claro en la siguiente figura. Una variación de la resistividad agudo está permitido en el método de inversión utilizado entre las dos secciones.

Figura 21. Disposición de células de modelo y el límite inferior de agua. La aparente pseudosección resistividad y el modelo de inversión para este conjunto de datos se muestran en la Figura 22. En esta inversión, las variaciones de resistividad dentro de la capa de agua se reducen al mínimo. Por lo tanto la resistividad capa de agua es generalmente uniforme excepto hacia el extremo izquierdo de la línea de levantamiento donde no hay muchos datos. Hay también una indicación de que entre los 190 y 220 metros marcas, el material de lecho del río cerca de la superficie consiste en sedimentos de baja resistividad.

72

Figura 22. Inversión modelo para la encuesta con electrodos flotantes a lo largo del a lo largo del río San Francisco en Brasil. La línea delgada negro indica el límite inferior de agua.

Encuestas con electrodos flotantes se llevan a cabo con frecuencia usando un cable con un número de electrodos tirados detrás detrás de un barco, por lo general con un sistema de resistividad res istividad metros metros multicanal controlado por un enlace informático a un sistema GPS. A 5 metros separación entre los electrodos es takeouts utilizan comúnmente juntos un tipo dipolo-dipolo de configuración de la matriz (no es necesario simétrico). El barco se está moviendo continuamente y las mediciones se pueden realizar a intervalos arbitrarios arbitrarios aunque la longitud del dipolo es de 5 metros. Con frecuencia, las mediciones son digitalizadas para cada posición 1 metro. Por lo tanto, las posiciones de los electrodos activos son de 1 metro de separación que es mucho menor que la longitud del dipolo de 5 metros. Esto puede resultados en las líneas de levantamiento muy largos, y los modelos de inversión muy grandes con los ajustes por defecto del programa de inversión. inversión.

La resolución de la configuración de la encuesta es poco probable que sea mucho mejor que la longitud del dipolo en particular si hay capa de agua de espesor considerable entre el cable y el subsuelo. El programa normalmente normalmente utiliza un modelo en el que la anchura de las células normales se establece para ser el mismo (o media) como la distancia entre las posiciones de los electrodos a lo largo de la línea. Esto es razonable para los levantamientos topográficos normales normales con una configuración de cable estático. Para efectuar el reconocimiento reconocimiento electrodos flotantes, ya que la distancia entre las posiciones de los electrodos (1 metro) es mucho menor que la longitud del dipolo (5 metros), esto da como resultado es un modelo que es demasiado fino. Para evitar esto, utilice la opción prevista en el 'Inversion - Modelo Discretización - Cambiar el ancho de los bloques' menú. El siguiente cuadro de diálogo se muestra en esta opción. Dado que la resolución horizontal de los datos es poco probable que sea mejor que la mitad de la longitud del dipolo, se recomienda que la anchura de las células modelo fijado en 3 veces la separación entre electrodos unidad unidad (para la situación con 5 metros dipolos y espaciamiento de datos 1 metro ). Esto reducirá en gran medida el número de células modelo en el modelo de inversión.

73

El tiempo de inversión puede reducirse aún más utilizando la 'inversión - Métodos y configuración de Inversión inversión rápida de líneas de levantamiento largos' opción de menú que mostrará el siguiente cuadro de diálogo.

74

La primera opción es usar un método optimizado para el cálculo de la matriz jacobiana. Esto puede reducir significativamente significativamente el tiempo de cálculo sin afectar a la precisión de los valores, y debe dar el mismo m odelo de inversión como el método predeterminado. El segundo método es utilizar un método rápido para calcular la matriz jacobiana. Los valores de matriz jacobiana y modelo de este modo la inversión será ligeramente diferente de la obtenida utilizando el método estándar, pero la diferencia es normalmente pequeño (menos de 5%) y no afecta a la interpretación del modelo. El tercer método, utilizando técnicas de inversión escasa, está destinado a líneas de levantamiento incluso más largos de 1000 electrodos de posiciones o más. Este método todavía producirá un modelo contiguo a lo l argo de toda la línea de encuesta, pero se aprovecha de la escasez de las matrices utilizando en el cálculo de inversión para reducir el número de cálculos necesarios. Una vez más los resultados serán marginalmente diferente (por lo general menos de 5%) forma la obtenida usando el método estándar, pero no afectará a la interpretación de los resultados. El uso de los tres métodos pueden reducir el tiempo de inversión de aproximadamente una cuarta parte de los que requieren el uso de métodos estándar.

La mayoría de los sistemas de estudios con electrodos flotantes digitalizar digitalizar las posiciones de los electrodos a cada 1 metro, mientras que la distancia entre los electrodos takeouts en el cable es de 5 metros. Esta es probablemente probablemente demasiado bien cuando el espesor de la capa de agua es más que un par de metros. A discretización discretización más gruesa, tal como para cada

1,25 o 2,5 metros (es decir, una cuarta parte o l a mitad de la separación entre los electrodos takeouts), es probablemente suficiente. Esto tiene la ventaja de reducir el número de posiciones de los l os electrodos en el modelo de inversión, y por lo tanto reduciendo en gran medida el tiempo de inversión sin afectar significativamente significativamente a los resultados. Si el programa de conversión de datos suministrado con el sistema de encuesta no tiene una opción para cambiar el intervalo de discretización de datos, puede hacerlo cambiando el valor de separación entre electrodos unidad en la línea 2 del archivo de datos ( por ejemplo, de 1,0 a 2.5). En este caso, antes de leer en el archivo de datos, debe seleccionar la opción 'Archivo - Reúne a las posiciones de los electrodos' opción que mostrará el siguiente cuadro de diálogo.

75

 Al seleccionar la opción 'tamizar los electrodos digitalizará las posiciones de las mediciones a la unidad de separación de los electrodos se especifica en el archivo de datos.

Como nota final, se ha observado que, en algunos conjuntos de datos hay falta de coincidencia entre las coordenadas de los electrodos para las mediciones y la topografía del fondo como se mide usando un un sistema GPS de perfiles y. Esto se muestra cuando se utiliza la opción de incorporar la capa de agua en el modelo de inversión por un desplazamiento horizontal (por lo general de 5 a 10 metros) e ntre la topografía del fondo y las estructuras empinadas en el modelo de inversión tales como pequeños barrancos y picos. Esto es probablemente debido a las diferencias en las ubicaciones del sistema de sonda y el cable con los electrodos.

76

 Apéndice F: inversión de datos IP Los datos recogidos de IP (polarización inducida) encuestas consta de dos conjuntos, las mediciones de resistividad aparente normales y una aparente conjunto de datos IP. Por lo tanto, el número de puntos de datos es el doble del número para un estudio de resistividad normal. Un método posible invertir un conjunto de datos tal conjunto es para el tratamiento de los dos conjuntos de datos como completamente conjuntos de datos s eparados. La inversión de los datos de resistividad se completa antes de comenzar en la inversión del conjunto de datos IP. En este caso, el modelo de la resistividad obtenida en la última iteración de la inversión del conjunto de datos de resistividad (que por lo general tiene el error RMS más baja) se utiliza en la inversión del conjunto de datos IP. Este es probablemente el método más eficiente si el modelo de resistividad es una buena representación de la verdadera resistividad del subsuelo. Sin embargo, se ha encontrado que incluso en la inversión de ordenador genera conjuntos de datos, el modelo obtenido puede contener distorsiones en particular si se utiliza un número grande (más de 5) de iteraciones en la inversión del conjunto de datos de resistividad. Esto se debe a que la inversión se basa en un número finito de puntos de datos. La distorsión se vuelve peor cuando los datos contienen ruido. Muchos usuarios tienden a utilizar un gran número de iteraciones  junto con pequeñas factores de amortiguación en un esfuerzo por redu cir el error R MS para el valor más bajo posible. Si un modelo de resistividad con distorsiones se utiliza en la inversión de los da tos IP, entonces las distorsiones afectarán a los modelos de PI para todas las iteraciones de la inversión IP.

El programa ofrece dos métodos para invertir los datos de IP. En el primer método la inversión de los datos de resistividad e IP se llevan a c abo simultáneamente. Inmediatamente después de una iteración de la inversión de los datos de resistividad, una iteración de la inversión de IP se lleva a cabo. En el segundo método, la inversión de los datos IP se lleva a c abo sólo después de la finalización de la inversión de resistividad. En este caso, sólo la resistividad del modelo obtenido en la última iteración de resistividad se guarda.

los

archivos

IPMODEL.DAT,

IPMAGUSI.DAT,

IPSHAN.DAT y

IPKENN.DAT son ejemplos de archivos de datos con datos IP tanto la resistividad y. Este programa es compatible con cuatro tipos diferentes de datos IP; (I) mediciones de cargabilidad dominio del tiempo, (ii) mediciones de efecto ciento dominio de la frecuencia de frecuencia, (iii) las mediciones de ángulo de fase y (iv) los valores de f actor metálico IP. El archivo IPMODEL.DAT tiene los datos IP como valores de cargabilidad. La primera parte de los archivos con los comentarios son los siguientes: -

archivo IPMODEL.DAT

comentarios

Fault y el bloque de modelo

| Título

1.00 3 432 1 1

| Unidad de separación de los electrodos | tipo de matriz | Número de puntos de datos

| Ubicación de centro de matriz dada | 1 para indicar IP presente

77

cargabilidad

| Tipo de datos IP

mseg

| unidad de IP

0.1,1.0

| Retardo, tiempo de integración

1,50 1,00 1 12,04 7,2038 2,50 1,00 1 12,03 7,1983

| x-loc., a, n, app. Res., aprox. IP | Segundo punto de datos

Un valor de “1”, en lugar de “0”, se introduce en l a sexta línea para indicar que los datos I P está presente. Esto

es seguido por las líneas 3 de datos que pr oporciona información sobre la naturaleza de los datos IP. La séptima línea con la palabra cargabilidad muestra que los datos IP se da en términos de capacidad de carga. Muchos instrumentos IP medidos la capacidad de carga sean integrando el área bajo la curva de caída IP. En este caso, la unidad está en mseg (milisegundo). El valor de capacidad de carga obtenido por este método es normalmente calibrado a la Newmont M (

331)

 estándar (Summer 1976), de modo que el valor de cargabilidad en

mseg tiene el mismo valor numérico que la capacidad de carga dada en mV / V. La cargabilidad en mV / V se define como la relación de la tensión secundaria inmediatamente después de la corri ente se corta a la tensión de DC primario.

encuestas IP se han utilizado tradicionalmente en la industria de la exploración de minerales, en particular para los sulfuros de metal, donde se usan generadores que producen altas corrientes del orden de 10 amperios. Los valores de IP aparentes a partir de tales estudios son por lo general menos de 100 mseg. (O mV / V). Un desarrollo reciente es la adición de la capacidad de IP para sistemas basados en baterías utilizados en los estudios de ingeniería y ambientales en los que normalmente corrientes de 1 amperio o menos. Un fenómeno acompañante es la observación de los valores de IP de más de 1.000 mseg. (O menos de -1,000 mseg.) En algunos conjuntos de datos. Tales valores son casi con toda seguridad causados por el ruido debido a una señal muy débil IP. Para comprobar si tales altos valores de IP son reales, comprobar primero la aparente pseudosección resistividad. Si muestra valores inusualmente altos y bajos que varían de un modo irregular, los datos es ruidoso. Si los valores de resistividad aparente son ruidosos, entonces los valores IP aparentes son casi ciertamente poco fiable. A continuación comprobar el pseudosección IP aparente. Si los valores de IP aparentes muestran un patrón irregular (con frecuencia con valores anómalos alineados en diagonal con un vértice en un electrodo de dudoso), entonces los valores de IP son demasiado ruidoso para s er interpretable.

valores de IP que se dan en otras unidades se convierten en valores de cargabilidad internamente por el programa durante el proceso de inversión. Las ecuaciones de c onversión pertinentes pueden encontrarse en los trabajos de Van Voorhis et al. (1973) y Nelson y Van Voorhis (1973).

Para medidas de dominio de frecuencia, el valor IP se da generalmente en Efecto Porcentaje Frecuencia (PFE). Esto se da como PFE = 100 ( ρ

corriente continua -

ρ  AC) / ρ C.A.

dónde ρ corriente continua y ρ C.A. son los valores de resistividad aparente medidos en las frecuencias bajas y al tas. El

IPSHAN.DAT archivo es un ejemplo de un dominio de frecuencia IP

78 conjunto de datos. En la séptima línea, el tipo de datos IP se da como “Porcentaje Frecuencia Efecto”. La línea

de datos novena contiene los valores de las frecuencias bajas y altas utilizadas en la medición. Este co njunto de datos es de un estudio de campo en Birmania, donde las mediciones se hicieron con 20 metros y 40 metros de dipolos (Edwards 1977). Al igual que muchos conjuntos de datos de campo, tiene una distribución muy compleja de los puntos de datos en el pseudosección. Se trata de un conjunto de datos relativamente ruidoso. Es recomendable utilizar la opción “Incluir suavización de valores de resistividad modelo” que se describen en la sección 9 de la inversión de este conjunto de

datos para evitar valores extremadamente grandes y pequeñas modelo de resistividad. El objetivo de la encuesta es un depósito de galena en meta-sedimentos. La geología es bastante complejo con enormes depósitos de galena junto con vetas de cuarzo, mármoles silicificada y zonas de arcilla. Hay una gama muy amplia de valores de resistividad.

valores de IP también se dan como ángulos de fase. En este caso la unidad que se utiliza normalmente es mili-radianes (mrad). Un ejemplo de un conjunto de datos IP establecida con ángulos de fase se da en el archivo IPKENN.DAT con los valores del papel por Hallof (1990). En la séptima línea, el tipo de tipo de datos IP se da como “ángulo de fase”. El segundo parámetro en la línea de datos novena contiene los valores de la

frecuencia a la que se midieron los ángulos de fase. El primer parámetro no se utiliza y normalmente se establece en 0. Muchos instrumentos de encuesta IP moderna puede medir el ángulo de fase para varias frecuencias. En este caso, la inversión de se repite para los datos de á ngulo de fase en cada frecuencia. Esto le dará una idea de la variación del ángulo de fase del subsuelo IP con frecuencia para estudios de resistividad complejas.

datos IP a veces se da como un valor de factor de metal (MF). El valor MF se puede calcular a partir de mediciones de dominio o bien el dominio del tiempo o de frecuencia. En el dominio del tiempo, el factor de metal está valor viene dado por  MF = 1,000 M / ρ

corriente continua

donde la capacidad de carga M se da en térm inos de milisegundo (Witherly y Vyselaar, 1990). En el dominio de la frecuencia, se da como MF = 100000 ( ρ

corriente continua - ρ  AC) / ρ C.A.

2

El archivo es IPMAGUSI.DAT de una encuesta sobre el cuerpo de mineral de Río Magusi (Edwards, 1977), donde la encuesta se realizó con 30,5 metros (100 pi es), 61.0 metros (200 pies) y 91,4 metros (300 pies) de dipolos. De nuevo el pseudosección resultante tiene una distribución muy compleja de los puntos de datos. La resistividad e IP pseudosecciones aparentes medidos,  junto con las secciones de modelo obtenidos se muestran en la Figura 23. El cuerpo de mineral se muestra como un cuerpo de baja resistividad distinta con altos valores de PI cerca de la mitad de l a línea de levantamiento en secciones de modelo. En la inversión de este conjunto de datos, se utilizó el método de inversión robusta (ver Apéndice K) para afilar el límite entre el cuerpo de mineral y las rocas circundantes.

79

El programa de inversión trata de minimizar la diferencia entre los valores de cargabilidad aparentes equivalentes observadas y calculadas. valores de PI en efecto frecuencia ciento y ángulos de fase son directamente proporcionales a los valores de cargabilidad equivalentes. Sin embargo, el valor del factor de metal también implica un valor de resistividad. El modelo que da un error RMS mínimo para los valores de cargabilidad a veces no da un error mínimo RMS para los valores de factor de metal. Sin embargo, el efecto real sobre el modelo de inversión es generalmente muy pequeñas.

Figura 23. río Magusi cuerpo de mineral (a) pseudosección resistividad aparente, (b) sección del modelo de resistividad, (c) factor de aparente de metal pseudosección y la sección (d) Factor de modelo de metal.

80

 Apéndice G: Nuevas opciones de modo de proceso por lotes Esta opción está destinada a usuarios que deseen invertir una gran cantidad de conjuntos de datos en una sola vez por control remoto. Se le permite iniciar el modo por lotes desde otro programa de Windows, o desde una línea de comandos. En la opción modo batch un archivo de comandos, por ejemplo RESIS.BTH, se utiliza para informar al programa de los archivos de datos que deben ser procesados y las diversas opciones de inversión utilizados. En Windows, puede acceder a la opción de línea de comandos haciendo clic en el botón “Inicio”, seguido por la opción “Ejecutar”. Si el programa RES2DINV se encuentra en el directorio C: subdirectorio \

R2DINV, y el archivo de secuencia de comandos por lotes RESIS.BTH se encuentra en el D : subdirectorio \ DATA, el siguiente comando se escribe en el cuadro de diálogo 'Ejecutar”.

C: \ R2DINV \ RES2DINV D: \ DATA \ RESIS.BTH Este comando pondrá en marcha el programa RES2DINV que procesará automáticamente la lista de archivos de datos en el RESIS.BTH. En el archivo RESIS.BTH, se recomienda que se escribe en la unidad y la ruta co mpleta para todos los archivos. Después de que el programa ha procesado todos los archivos de d atos, se apagará automáticamente. Puede procesar hasta 40 archivos de datos utilizando un único archivo de c omandos por lotes.

Desde la versión 3.53, un nuevo y más sencillo formato de archivo para la opción de modo de lote que hace uso de los archivos de parámetros de inversión (PIV) que son creados por el programa RES2D INV sí. El RESIS_IVP.BTH archivo es un ejemplo, con el siguiente formato. archivo RESIS_IVP.BTH

| comentarios

----------------------------------------------------------------------------------------------3

| Número de archivos para invertir

INVERSIÓN DE PARÁMETROS archivos usados

| Cabecera para indicar nuevo modo po r lotes

 ARCHIVO DE DATOS 1

| Encabezado por el primer archivo

C: \ RES2DINV \ RATHCRO.DAT

| archivo de datos para invertir

C: \ RES2DINV \ RATHCRO_L2.INV

| archivo de salida con resultados de la inversión

C: \ RES2DINV \ RES2DINV_L2.IVP

| Archivo que contiene la configuración de inversión

 ARCHIVO DE DATOS 2

| Encabezado por el segundo archivo

C: \ RES2DINV \ RATHCRO.DAT

| archivo de datos para invertir

C: \ RES2DINV \ RATHCRO_L1.INV

| archivo de salida con resultados de la inversión

C: \ RES2DINV \ RES2DINV_L1.IVP

| Presentar con diferentes ajustes de inversión

 ARCHIVO DE DATOS 3

| Encabezado por el tercer archivo

C: \ RES2DINV \ LANDFILL.DAT

| En tercer archivo de datos para invertir

C: \ RES2DINV \ LANDFILL_L2.INV

| archivo de salida con resultados de la inversión

C: \ RES2DINV \ RES2DINV_L2.IVP

| Archivo que contiene la configuración de inversión

Dado que los archivos de IVP son creados por el programa RES2DINV, no es necesario que el usuario lidiar con los detalles sobre el formato del archivo. Para crearlos, el usuario sólo tiene que hacer cli c en las opciones de inversión adecuadas dentro del programa RES2DINV, y luego seleccione la opción “Guardar parámetros de inversión bajo el‘el menú Cambiar configuración’para guardar los parámetros en un archivo de

PIV.

81

 Apéndice H: Velocidad de visualización de gráficos

Como RES2DINV es un programa basado en Windows, las operaciones de gráficos por el programa son manejados a través del controlador de pantalla de Windows para la tarjeta de gráficos por ordenador. Así, el programa funciona con cualquier tarjeta gráfica con tal de que tiene un controlador de Windows. Esto es diferente de los programas basados en DOS donde los problemas de compatibilidad con la alta resolución de los modos gráficos Super VGA son comunes. Una gran ventaja de este programa es que los modos gráficos de alta resolución, tales como los 1024 por 768 píxeles con modo de gráficos de 32 bits, pueden ser utilizadas fácilmente por este programa. Esto es importante para los datos obtenidos de las encuestas con más de varios cientos de electrodos (una encuesta con cerca de 650 electrodos y casi 5200 puntos de datos ahora tiene el récord de encuestas de la tierra) en el que el modo VGA estándar de 640 por 480 píxeles de gráficos es claramente insuficiente.

Dado que las operaciones de gráficos se llevan a cabo a través del controlador de pantalla, la velocidad de la pantalla gráfica depende de la eficiencia del controlador de pantalla. Para algunas tarjetas gráficas, las versiones anteriores de los controladores son relativamente ineficientes que se traduce en operaciones de gráficos dolorosamente lentos como en la elaboración de los pseudosecciones. Si esto sucede en su ordenador, la causa es casi seguro que debido a un controlador de gráficos de versión temprana de la tarjeta de visualización. La mejor manera de superar este problema es reemplazar el controlador de la tarjeta de gráficos de Windows con una nueva versión que debería ser más eficiente. Los controladores se pueden descargar a través de I nternet desde muchos sitios de shareware o en el sitio web del fabricante. Con la última versión del controlador de pantalla puede dar lugar a mejoras significativas en la velocidad de las operaciones de gráficos.

82

 Apéndice I: inversión de datos cruzada pozo La resolución de encuestas eléctricos llevado a cabo con electrodos en la superficie del suelo disminuye exponencialmente con la profundidad. Un método para obtener ra zonablemente buena resolución en profundidad es haciendo mediciones con los electrodos en pozos de sondeo. encuestas de corte del pozo se llevan a cabo con mucha menos frecuencia que las encuestas basadas en la superficie, por lo que a diferencia de las encuestas de superficie, en la actualidad hay falta de normas comunes con respecto a la disposición de los electrodos. La opción transversal pozo de sondeo en este programa está destinado a encuestas poco profundas con la disposición de los electrodos colocados en un cierto patrón general (Sasaki, 1992). La disposición básica de los electrodos se muestra en la Figura 24.

Los electrodos se dividen en 3 grupos con (i) los electrodos de superficie, (ii) los electrodos en la pared del pozo y (ii) electrodos 1 en pozo 2. Se supone que las perforaciones no tienen una carcasa de metal que alteraría gravemente el patrón de el flujo de corriente. El programa divide el subsuelo en un número de bloques de cuadrilátero (Figura 24). Tenga en cuenta que las posiciones de la superficie y los electrodos de pozo controlan la forma en que el subsuelo se divide en los bloques de modelo. El BOREDIFF.DAT archivo es un archivo de datos de ejemplo, con mediciones crossborehole. Una descripción del formato utilizado se presenta a continuación en la sección de comentarios. Tenga en cuenta que en este ejemplo, las mediciones se realizan con la matriz de polo-dipolo, donde sólo se utilizan la C1, P1 y P2 electrodos.

archivo BOREDIFF.DAT

comentarios

electrodos de pozo a diferentes profundidades | Título

1.0 12 840 2 0

| Mas pequeña separación entre electrodos

Los electrodos de superficie

| Cabecera para electrodos de superficie

dieciséis

| Número de electrodos de superficie

0.0, 0.0 1.0, 0.0 2.0, 0.0 3.0, 0.0 4.0, 0.0 5.0, 0.0 6.0, 0.0 7.0, 0.0 8.0, 0.0 9.0, 0.0

| x y z-ubicación de primera electrodo de superficie

10.0, 0.0 11.0, 0.0 12.0, 0.0 13.0, 0.0

|

| número de matriz 12 para encuesta transversal del pozo | Número de puntos de datos | 2 para indicar la ubicación se utiliza el formato XYZ

| 0 para no datos IP

| ubicación del segundo electrodo de superficie | valores de ub icación de las notas se dan para el

| electrodos en orden ascendente de partida | desde la izquierda. | Tenga en cuenta un valor cero para la profundidad de la | electrodos de superficie.

| | | | | |

83

14.0, 0.0 15.0, 0.0

|

Pozo 1 Electrodos

| Encabezado por el primer pozo

10 4,0 1,0 4,0 2,0 4.0 3.0 4.0 4.0 4,0 5,0 4,0 6,0 4,0 7,0 4,0 8,0 4,0 9,0

| Número de electrodos en el primer pozo

|

| x y z-ubicación de primer electrodo | x y z-ubicación de segundo electrodo | Tenga en cuenta la ubicación de lo s electrodos

| se dan desde el electrodo superior | debajo de la superficie. | También tenga en cuenta una posición en x constante para

| pozo vertical. | |

4,0 10,0

|

Los electrodos de pozo 2

| Encabezado por el segundo pozo

10

| Número de electrodos en el segundo pozo

11,0 1,5 11,0 2,5 11,0 3,5 11,0 4,5 11,0 5,5 11,0 6,5 11,0 7,5 11,0 8,5 11,0 9,5

| x y z-ubicación de primer electrodo

11,0 10,5

|

Datos de medición

| Encabezado de sección con las mediciones

| a continuación se dan ubicación de otros electrodos.

| | | | | | |

3 0.00 0.00 1.00 0.00 2.00 0.00 101.5718 | El formato para cada punto de datos es: 0.00 2.00 3 0.00 0.00 3.00 0.00 99.5150 | Número de electrodos utilizados en

0.00 3.00 3 0.00 0.00 4.00 0.00 99.2303 | medición, x y z -ubicación de 0.00 4.00 3 0.00 0.00 5.00 0.00 99.1325 | electrodo C1, x y z-ubicación de 0.00 5.00 3 0.00 0.00 6.00 0.00 101.0616 | electrodo P1, x y z -localización de 0.00 6.00 3 0.00 0.00 7.00 0.00 105.7333 | electrodo P2, resistividad aparente 0.00 7.00 3 0.00 0.00 8.00 0.00 112.6745 | valor. 0.00 8.00 3 0.00 0.00 9.00 0.00 118.5223 | El número de matriz es de 12 para los datos de la encuesta transversal del pozo con las medidas dadas como valores de resistividad aparente. También puede introducir los datos como las mediciones de resistencia, en cuyo caso el número de matriz es de 13. Esto podría ser más conveniente que la mayoría de medidores de resistividad de obtener los valores como un valor de resistencia en ohm, y además el factor geométrico para matrices con electrodos del subsuelo es diferente del factor geométrico para matrices de superficie convencionales. El BORERES.DAT archivo da un ejemplo donde las mediciones se dan como valores de resistencia. establecer un conjunto de datos de campo interesantes donde se utilizó la matriz de polos y polos, y las mediciones también se dan como val ores de resistencia, es BORELUND.DAT de la Universidad de Lund, Suecia.

84

Figura 24. Disposición de los electrodos en una encuesta transversal pozo de sondeo. (A) El valor predeterminado modelo estándar donde las dimensiones de los bloques de modelo son igual a la separación entre los electrodos. (B) Un modelo alternativo con bloques más finas con dimensiones aproximadamente la mitad de la separación entre los electrodos.

Si sólo se utilizan dos electrodos en una medición, sól o la x e z ubicaciones de los electrodos C1 y P1 se dan en el archivo de datos (ver el BOREHIP.DAT archivo

85

que también contiene las mediciones de IP). Sin embargo, si se utilizan los cuatro electrodos, el X y Z-ubicaciones de los electrodos de C1, C2, P1 y P2 se debe dar en este orden. Se pue den combinar las mediciones realizadas con diferente número de electrodos en el mismo archivo de datos.

Cuando los programas se lee en un archivo con valores de resistencia, usted tiene la opción de invertir el conjunto de datos utilizando los valores de resistividad aparente, o directamente usar los valores de resistencia. El uso de valores de resistencia directamente en la inversión tiene la ventaja de que le permite utilizar las lecturas donde no existe el valor de resistividad aparente, o es negativo. Después de leer en un archivo de datos, el programa intentará filtrar las lecturas sospechosas con niveles potencialmente altos de ruido si las medidas se expresan como valores de resistividad aparente, o si había elegido la opción de utilizar valores de resistividad aparente en la inversión. Si decide llevar a cabo la inversión utilizando valores de resistencia, las lecturas no se fil tran.

El factor geométrico para los electrodos del subsuelo es diferente de la utilizada para electrodos de superficie. Como ejemplo, el factor geométrico para una medición en el que se usan sólo dos electrodos es como sigue. Si los electrodos C1 y P1 se encuentran en ( x 1, z 1) y ( x 2, z 2) respectivamente, el factor geométrico k a continuación, está dada por  k = 4 pi / (1 / r

 1 / r

1 +

1' )

dónde

r 1 = sqrt ( dx 2 + dz 2)

r 1 '' = sqrt ( dx

dx = x 1 - x 2

dz = z 1 - z 2

2+

dz 2)

Dz = z 1 + z 2 ecuaciones similares para mediciones con 3 o 4 electrodos se pueden obtener fácilmente mediante la adición de los términos apropiados. La profundidad de los electrodos en Pozo 1 puede ser diferente de los electrodos correspondientes en perforación 2. En la práctica, para obtener los mejores resultados las profundidades no debe diferir demasiado. Este programa requiere que el número de electrodos en ambos pozos debe ser el mismo. En la práctica, esta restricción no es un problema significativo. Puede ser fácilmente superar esta restricción mediante la inserción de electrodos ficticias en el conjunto de datos si el número de electrodos en un pozo de agua es menor (en la sección donde se les da la z-ubicación xy de los electrodos en el pozo de sondeo).

El programa también requiere que un electrodo de superficie se encuentra en la parte superior de ambos pozos. Si no existen físicamente, basta con insertar electrodos ficticias en los lugares apropiados en el archivo de datos. También debe haber por lo menos un par de electrodos de superficie entre los dos pozos de sondeo, así como al menos dos electrodos de superficie a la izquierda de pozo de sondeo 1, y otras dos a la derecha del pozo de sondeo 2. Si tales

86

electrodos no se utilizaron en el estudio de campo, basta con insertar electrodos ficticias en el a rchivo de datos para satisfacer los requisitos de este programa.

Por defecto, el programa subdividir el subsuelo en bloques de cuadrilátero tales que las dimensiones de los bloques son iguales a la separación entre los electrodos (Figura 24a). Sin embargo, también se puede optar por utilizar un modelo en el que las dimensiones de los bloques es de aproximadamente la mitad de la separación entre los electrodos (Figura 24b). Usando el modelo con una discretización más fina podría mejorar significativamente significativamente los resultados de la inversión en algunos casos (Sasaki, 1992). Para seleccionar el tipo de modelo a utilizar, haga clic en las opciones de la 'Configuración de cambio' en la barra superior del menú principal, seguido de 'Tipo de modelo cruzada pozo'. En teoría, es posible subdividir aún más el subsuelo en bloques más pequeños. Sin embargo, dado que la resolución de las mediciones eléctricas disminuye rápidamente con la distancia de los electrodos de medición, esto es poco probable que mejorar significativamente significativamente los resultados. De hecho, que podría resultar en oscilaciones osci laciones no naturales en los valores de resistividad modelo, particularmente particularmente cerca de los electrodos. Esto es como resultado de la compensación conocida entre la resolución de modelo y modelo de variación en la inversión geofísica.

La Figura 25 muestra el modelo obtenido después de

la inversión de

el

BOREHOLE.DAT conjunto de datos. Se trata de un conjunto de datos sintético producido por un programa de modelado hacia adelante con un bloque de alta resistividad entre los dos pozos de sondeo. La opción de utilizar bloques más finos con la mitad del tamaño de la separación entre los electrodos se utilizó en la inversión de este conjunto de datos.

La Figura 26 muestra los resultados de la inversión de un conjunto de datos de campo interesante. Este conjunto de datos es uno de un número que se recogieron por una encuesta para estudiar el flujo de fluidos a través del acuífero UK Chalk en el este de Yorkshire mediante el uso de un trazador de solución salina (Slater et al 1997). Hay una zona de baja resistividad cerca de la superficie, donde la solución salina se irrigó en el suelo, y también prominentes zonas de baja resistividad por debajo de una profundidad de 7 metros debido a la trazador solución salina que había fluido hacia abajo. Otra característica interesante de este conjunto de datos es que no se utilizaron electrodos de superficie.  Así que todos los electrodos de superficie que figuran en el archivo de datos BORELANC.DAT BORELANC.DAT son electrodos ficticias; dos a la izquierda de pozo de sondeo 1 y otros dos a la derecha del pozo de sondeo 2.

Si también se hicieron mediciones de IP, es necesario introducir un valor de 1 en la sexta línea en el archivo de datos, y dos más líneas de datos después de aquella que da el tipo de medición IP utilizada y los dos parámetros asociados a ella (véase el apéndice F). A continuación, introducir el valor de IP después de que el valor de resistividad aparente para cada punto de datos. El BOREHIP.DAT archivo es un archivo de datos a través del pozo de sondeo ejemplo con las mediciones de IP.

87

Figura 25. El modelo obtenido después de la inversión del conjunto de datos BOREHOLE.DAT. BOREHOLE.DAT. Se usó la opción de utilizar bloques más finos con la mitad del tamaño de la separación entre los electrodos.

Figura 26. Modelo obtenido a partir de la inversión de datos de una encuesta transversal pozo para mapear el flujo de un trazador de solución salina en medio de dos pozos de sondeo. Tenga en cuenta las zonas de baja resistividad cerca de la superficie (donde se inyectó el trazador) y por debajo de una profundidad de 7 metros. Las ubicaciones de los electrodos de pozo se muestran con pequeños puntos negros.

88

El archivo de datos de ejemplo boreholes_time_2D.dat muestra el formato de datos para 2-D datos de las perforaciones con lapso de tiempo. El formato es similar a la utilizada para los estudios superficiales normales. La sección de datos con el encabezado de lapso de tiempo se da directamente después de la sección con la información de pozo, e inmediatamente antes de los datos de resistividad aparente. Las opciones de i nversión para los datos de lapso de ti empo de sondeo se establecen del mismo modo que para las encuestas superficiales normales usando l as mismas opciones de menú. La siguiente tabla muestra una parte del archivo de datos con la información en tiempo lapso.

boreholes_time_2D.dat -----------------------------------------------------------------

Comentario

44.000, 64.000

Última electrodo de pozo

datos de la secuencia Time

Cabecera para indicar los datos de lapso de

-------------------------------------

tiempo Número de secciones de tiempo

Encabezamiento

2

Número de serie de tiempo

unidad de tiempo

Encabezamiento

horas

Tipo de unidad de tiempo

intervalo de la sección segunda vez

Encabezamiento

1.0

Tiempo de la segunda serie

Datos de medición

Cabecera de la sección de datos principal

4

4.000, 0.000 0.000, 0.000 0.000 101,4981 101,4983 0.000 12.000,

4

8.000, 0.000 4.000, 0.000 0.000 16.000, 0.000 101,3597 101,3599

......

8.000,

Primera medición, tenga en cuenta dos va lores de resistividad aparente

12.000,

segunda medición Lo mismo para el resto

El archivo boreholes_inclined.dat da un ejemplo del formato de datos con perforaciones no verticales. Tenga en cuenta la adición de una línea de cabecera adicional 'pozos inclinados presente' justo antes de la sección con las posiciones de los pozos de sondeo para indicar que las perforaciones no verticales están presentes. Para pozos inclinados, la subdivisión de la subsuperficie en bloques modelo ya no es ligado directamente a las posiciones de los electrodos. También puede utilizar este formato para pozos verticales son do nde los electrodos correspondientes a los diferentes pozos se encuentran a diferentes profundidades. Esto evita bloques modelo inclinadas o sesgadas cuando hay diferencias significativas en los electrodos profundidades del pozo.

El programa RES2DINV soporta hasta 10 pozos de sondeo (20 perforaciones si el equipo tiene al menos 256 MB de RAM) a lo largo de una línea encuesta. Este límite se podría aumentar en el futuro cuando surja la necesidad. Una posible situación en la que más de 2 “perforaciones” “perforaciones” existe es cuando un electrodo de

corriente se coloca en el extremo de un perforador de pozo de sondeo o penetrómetro. El archivo BOREHOLE3.DAT da un ejemplo del formato de datos con más de 2 pozos. Una descripción del formato es el siguiente.

89

comentarios

archivo BOREHOLE3.DAT

--------------------------------------------------------------------------------------------prueba de pozo polo-dipolo

| Título

1.0 12 908 2 0

| Mas pequeña separación entre electrodos

Los electrodos de superficie

| Cabecera para electrodos de superficie

dieciséis

| Número de electrodos de superficie

0.0, 0.0 1.0, 0.0 2.0, 0.0 3.0, 0.0 4.0, 0.0 5.0, 0.0 6.0, 0.0 7.0, 0.0 8.0, 0.0 9.0, 0.0

| x y z-ubicación de primera electrodo de superficie

10.0, 0.0 11.0, 0.0 12.0, 0.0 13.0, 0.0 14.0, 0.0 15.0, 0.0

|

Número de pozos

| Nueva cabecera durante más d e 2 pozos de sondeo

3

| Número de perforaciones aquí

Pozo 1 Electrodos

| Encabezado por el primer pozo

11

| Número de electrodos en el primer pozo

| número de matriz 12 para encuesta transversal del pozo | Número de puntos de datos | 2 para indicar la ubicación se utiliza el formato XYZ

| 0 para no datos IP

| ubicación del segundo electrodo de superficie | valores de u bicación de las notas se dan para el

| electrodos en orden ascendente de partida | desde la izquierda. | Tenga en cuenta un valor cero para la profundidad de la | electrodos de superficie.

| | | | | | | |

4,0 1,0 4,0 2,0 4.0 3.0 4.0 4.0 4,0 5,0 4,0 6,0 4,0 7,0 4,0 8,0 4,0 9,0

| x y z-ubicación de primer electrodo | x y z-ubicación de otros electrodos

4,0 10,0 4,0 11,0

|

| | | | | | | |

Los electrodos de pozo 2

| Encabezado por el segundo pozo

11

|

11,0 1,0 11,0 2,0 11,0 3,0 11,0 4,0 11,0 5,0 11,0 6,0

| | | | | |

90

11,0 7,0 11,0 8,0 11,0 9,0

|

11,0 10,0 11,0 11,0

|

Pozo 3 Electrodos

| Encabezado por el tercer pozo

11

|

13,0 1,0 13,0 2,0 13,0 3,0 13,0 4,0 13,0 5,0 13,0 6,0 13,0 7,0 13,0 8,0 13,0 9,0

|

13,0 10,0 13,0 11,0

|

Datos de medición

| Los valores de resistividad aparente comienzan después de este

| | |

| | | | | | | | |

3 0.00 0.00 1.00 0.00 2.00 0.00 101.5718 3 0,00 0,00 2,00 0,00 3,00 0,00 99.5150

| Nota 3 electrodos para | una matriz de polos bipole

. .

3 11,00 11,00 11,00 5,00 11,00 4,00 110.7175 3 11,00 11,00 11,00 4,00 11,00 3,00 118.4761 0,0,0,0,0

| Datos de la última dos

| puntos

| Terminar con 5 ceros

Tenga en cuenta que el formato es esencialmente el mismo que el utilizado anteriormente para dos pozos de sondeo, a excepción de la cabecera después de la lista de los electrodos de superficie y el valor para el número de pozos de sondeo.

El BOREHOLE_TOPO.DAT archivo es un ejemplo de un campo de datos de corte del pozo se establece con la topografía. La parte inicial de este archivo de datos con una descripción del formato es el siguiente.

91

archivo BOREHOLE_TOPO.DAT

comentarios

POZO con topografía | Título 1.00000

| Mas pequeña separación entre electrodos

12

| número de matriz 12 para encuesta transversal del pozo

570 2 0

| Número de puntos de datos

topografía actual

| Cabecera para indicar topografía está presente

| 2 para indicar la ubicación se utiliza el formato XYZ

| 0 para no datos IP

modo de información de la topografía (0 = WRT para electrodos de superficie) | cabecera 0 | Indica profundidades de electrodos de pozo se dan WRT | el e lectrodo de superficie en la parte superior del pozo de sondeo

Tipo de x-distancias

| Encabezamiento

2

| 2 para indicar x-distancias son a lo largo de la superficie

Los electrodos de superficie

| Encabezamiento

31 -14.0000 0.00 -13.0000 0.03 -12,0000 0,05 -11.0000 0.08 -10,000 0,11 -9.0000 0.13 . . 15.0000 0.85 16,0000 0,95

| Número de electrodos de superficie

Pozo 1 electrodos

| Encabezado por el primer pozo

17 -14.0000 1.0000 -14.0000 2.0000 . .

| Número de electrodos en pozo

| x-ubicación y la elevación del primer electrodo | x-ubicación y la elevación del segundo electrodo | Lo mismo para el resto

| | | | | | | electrodo de superficie última

| x-ubicación y profundidad de primer electrodo de pozo | x-ubicación y profundidad de primer electrodo de pozo | Lo mismo para el resto | El formato es el mismo que para los archivos de datos de las perforaciones normales

Tenga en cuenta que la información de la topografía se da en las coordenadas z de los electrodos de superficie.

92

electrodos remotos del polo polos y matrices polo-dipolo - Apéndice J En teoría, la matriz polo polos tiene sólo 2 electrodos, el electrodo de C1 corriente positiva y el electrodo de potencial positivo P1. En la práctica una matriz tal no existe ya que todas las mediciones de campo se realizan utilizando una matriz con 4 electrodos (Figura 27).

Figura 27. Un diagrama esquemático de la disposición de los electrodos en una medición polo polos de campo real.

Para minimizar los efectos de la C2 y electrodos P2, generalmente se recomienda que la distancia de estos dos electrodos de los electrodos C1 y P1 debe ser al menos 10 veces (y preferentemente 20 veces) el C1-P1 espaciado r 1. En algunos casos, cuando se utilizan grandes separaciones entre los electrodos C1 y P1, este requisito no se cumple para todas las mediciones. Esto podría conducir a distorsiones en el modelo de inversión (Robian et al. 1997). Para superar este problema, el programa RES2DINV permite al usuario incorporar los efectos de los electrodos C2 y P2 en la inversión.

Para calcular el valor de resistividad aparente medido con la matriz Polepole no id eal, dos factores geométricos posibles pueden ser utilizados. Un método es usar la m isma ecuación que la matriz polo polos ideal. Esto da un factor geométrico aproximada ya que el val or de resistividad aparente resultante no es la misma que la resistividad verdadera para un semiespacio homogéneo. Alternativamente, la exacta

93

factor geométrico que tiene en cuenta las posiciones de los electrodos C2 y P2 se puede utilizar. factor geométrico aproximado (array ideales polo polos) = 2.pi.r1 factor geométrico Exacta = 2 pi / [(1 / r1) - (1 / r2) - (1 / r3) + (1 / r4)] El programa es compatible con ambas convenciones.

Una posible disposición de los electrodos en una encuesta polo polos se m uestra en la Figura 28. En este caso, el C2 y electrodos P2 se fijan a lo largo de la encuesta. Las mediciones s e realizan usando diferentes pares de electrodos a lo largo de la línea de encuesta com o los electrodos C1 y P1. Para una línea encuesta con norte electrodos, hay n (n + 1) / 2 combinaciones posibles. Por ejemplo, si hay 25 electrodos en la línea de encuesta, hay 300 combinaciones posibles. En la práctica, no se hacen todas las posibles mediciones. En lugar de una distancia máxima entre los electrodos C1 y P1 se ajusta en función de la profundidad máxima de investigación necesario (Edwards 1977). En muchos casos, esta distancia máxima es de 6 a 10 veces l a distancia entre electrodos unidad a lo largo de la línea de encuesta. Para evitar la posibilidad de valores de resistividad aparente negativo o cero, la distancia de los electrodos C2 y P2 debe ser al menos 2,5 veces (y preferiblemente 3 veces) el máximo C1 al espaciamiento P1 utilizado. Por ejemplo, si la separación entre electrodos unidad es 1 metro, y el máximo espaciado C1-P1 es de 10 metros, entonces la C2 y electrodos P2 deben estar al menos a 25 metros de la línea de encuesta. Tenga en cuenta que como la distancia de los electrodos C2 y P2 desde la línea de encuesta disminuye, la profundidad máxima efectiva de investigación de la encuesta también disminuye. Puede ser que sea significativamente menor que la profundidad de investigación para la matriz de polo polos ideales (Edwards 1977).

Figura 28. Un diagrama esquemático de la disposición de los electrodos C2 y P2 y los electrodos a lo largo de la línea de levantamiento en una encuesta polo polos. Para incorporar los efectos de los electrodos C2 y P2 en la inversión de datos, su x e y las coordenadas con respecto a la línea de levantamiento deben ser registrados.

94

El POLPOLFX.DAT archivo es un archivo de datos de ejemplo, con mediciones de polos polos donde s e especifican las ubicaciones de los electrodos remotos C2 y P2. Una descripción del formato utilizado se presenta a continuación en la sección de comentarios. A pesar de los Z-localizaciones de los electrodos también se incluyen en el archivo de formato de datos para la integridad, los valores están en la actualidad no se utilizan por el programa. El programa en la actualidad estima las alturas en los electrodos remotos desde las alturas de los electrodos a lo largo de la línea de encuesta. Sin embargo, en el futuro, podrían ser utilizados para dar cabida a los casos extremos en los que la pendiente del terreno en los electrodos a distancia es muy diferente a la largo de la línea de la encuesta.

comentarios

archivo POLPOLFX.DAT bloques Modelo

| Título

2.00 2

| Mas pequeña separación entre electrodos

electrodos de mando a distancia incluidos

| Cabecera para indicar la posición del mando a distancia |

| número de matriz 2 para el poste-poste. electrodos incluidos en el archivo de datos

C2 electrodo remoto X, Y y

| Encabezado por segunda ubicación a ctual del electrodo

Z ubicación

|

-50.000,20.000,0.0 P2 electrodo remoto X, Y y

| X, Y y Z-ubicación de electrodo C2 | Encabezado por el segundo electrodo de potencial

Z ubicación

| ubicación

120.000,0.000,0.0

| X, Y y Z-ubicación de electrodo P2

factor geométrico exacto utilizado

| Tipo de factor geométrico utilizado

295 1

| Número de puntos de datos

0

| 0 para no IP

| 1 x para indicar la localización del centro de la matriz se | dado

1.00

2.00

10.13

| La posición en x, el espaciamiento de electrodos | resistividad aparente para el primer punto de datos

3.00

2.00

10.18

| Segundo punto de datos

El formato de datos es el mismo que para las encuestas normales polo polos, excepto para la sección (después del número de matriz) que contiene información sobre la ubicación de los segundos electrodos actuales y potenciales. La línea de datos “f actor geométrico exacto utilizado” se utiliza para indicar que los

valores de resistividad aparente en el archivo de datos se calcularon utilizaron el factor geométrico exacto. Si se utilizó el factor geométrico aproximado, esta línea de datos se debe administrar como “factor geométrico aproximado utilizado”. El POLDPLFA.DAT archivo es un archivo de datos de ejemplo para una encuesta polo-

dipolo donde se utilizó el factor geométrico aproximado en el cálculo de los valores de resistividad aparente. Tenga en cuenta que en este archivo, se da sólo la ubicación del segundo electrodo de corriente C2.

Para la matriz de polo polos, el efecto del electrodo de C2 es aproximadamente proporcional a la relación de la distancia C1-P1 a la distancia C2-P1. Si los efectos de los electrodos C2 y P2 no se toman en cuenta, la distancia de estos electrodos de la línea de la encuesta debe ser al menos 20 veces el más grande C1-

95

P1 espaciamiento utiliza para garantizar que el error es menos de 5%. En encuestas en los que la separación entre electrodos a lo largo de l a línea de la encuesta es de más de unos pocos metros, puede haber problemas prácticos en la búsqueda de ubicaciones adecuadas para el C2 y P2 electrodos para satisfacer este requisito. Por lo tanto el error en descuidar los efectos de l os electrodos a distancia es mayor para la matriz polo polos. Como guía general, si las distancias de la C2 y electrodos P2 de la línea de encuesta es más de 20 veces el mayor espaciamiento C1-P1 utilizado, la matriz puede ser tratada co mo una matriz Polepole ideal para los propósitos de interpretación. De lo contrario, las coordenadas de los electrodos C2 y P2 deben ser registrados de manera que sus efectos se pueden determinar por este programa. Sin embargo, hay que señalar que a medida que la distancia del electrodo P2 desde la línea de encuesta aumenta, el ruido telúrico recogido por el par P1-P2 también aumenta. Este factor también debe ser tomado en cuenta cuando se col oca el electrodo de P2.

Para la matriz de polo-dipolo, el efecto del electrodo de C2 es aproximadamente proporcional a la cuadrado de relación de la distancia C1-P1 a la distancia C2-P1. Por lo tanto la matriz polo-dipolo se ve menos afectada por el electrodo remoto C2. Si la distancia del electrodo C2 es m ás de 5 veces la mayor distancia C1-P1 utilizado, el error causado por despreciando el efecto del electrodo de C2 es menor que 5% (el error exacto depende también de la localización del electrodo P2 p ara la medición particular). Por lo tanto, como una guía general, si la distancia de electrodo de C2 desde la línea d e encuesta es más de 5 veces el mayor espaciamiento C1-P1, se puede tratar como una matriz polo-dipolo ideal. De lo contrario, introduzca las coordenadas del electrodo de C2 en el archivo de datos.

96

K Apéndice: datos sólidos y el modelo de inversión El método de mínimos cuadrados convencional intentará reducir al mínimo la cuadrado de la diferencia entre los valores de resistividad aparente medidos y calculados. Este método normalmente da resultados razonables si los datos contienen ruido aleatorio o “Gaussian”. Sin embargo, si el conjunto de datos contiene “atípicos” puntos de datos (donde el ruido proviene de fuentes no aleatorias tales errores o problemas con el equipo), este criterio es menos satisfactoria. Tales puntos de datos “atípico” podrían

tener una gran influencia en el modelo de inversión resultante. Para reducir el efecto de tales puntos de datos “atípicos”, un método de inversión donde el absoluto diferencia (o la primero energía) entre los valores

de resistividad aparente medidos y calculados se minimiza puede ser utilizado (Claerbout y M uir, 1973).

En el programa RES2DINV, este método de inversión puede ser seleccionado haciendo clic en la opción “Seleccionar la inversión sólida” en el menú “Cambiar la configuración”. En el cuadro de diálogo “Seleccionar la inversión robusta”, haga clic en los “ datos sólidos restringen” opción. Hay un factor de corte que controla el grado

en el que se utiliza este robusto Restringir datos. Si se utiliza un valor de 0.05, esto significa que el efecto de puntos de datos, donde las diferencias en los valores de resistividad aparente medidos y calculados son mucho mayores que 5 por ciento se reducirá considerablemente.

El método de mínimos cuadrados suavidad limitados convencional (deGrootHedlin y Constable 1990) también intenta minimizar el cuadrado de los cambios en los v alores de modelo de resistividad. Esto producirá un modelo con una variación suave en los v alores de resistividad. Tal modelo es el más adecuado, donde la resistividad del subsuelo también cambia de una manera suave, tal como el ejemplo pluma de contaminación muestra en la Figura 3. Sin embargo, si los cuerpos del subsuelo tienen límites afilados, como la interfaz suelo-roca o cuerpos homogéneos m asivas, la convencional método de los mínimos cuadrados suavidad-restricción tiende a manchar los límites.

Si se utiliza el método de modelo de restricción de inversión robusta en los valores de resistividad modelo, el programa intentará reducir al mínimo la absoluto cambios en los valores de resistividad. Esta restricción tiende a producir modelos con interfaces afilados entre las diferentes regiones con diferentes valores de resistividad, pero dentro de cada región el valor de la resistividad es casi constante. Esto podría ser más adecuado para las áreas donde existe tal situación geológica, tales como la interfaz suelo-roca sobre ciertos tipos de roca de fondo. Como ejemplo, la Figura 29 muestra los resultados de la inversión para un modelo sintético con un bloque de fallo (con una resistividad de 50 ohm.m) y un pequeño bloque rectangular (1 ohm.m) dentro de un medio circundante con una resistividad de 10 ohmios. metro. El modelo producido por el método de mínimos cuadrados estándar tiene un límite gradational para el bloque de fallo. En comparación, el modelo producido por el robusto método de modelo de inversión tiene límites más nítidas y más rectos.

97

Un ejemplo de campo con límites afilados se mostró anteriormente en la Figura 23 para la resistividad río Magusi y conjunto de datos IP. La estructura principal, el cuerpo masivo mineral sulfurado, por naturaleza tiene un contraste de resistividad / IP di stinta y aguda de las rocas ígneas país / metamórficas circundantes.

En el programa RES2DINV, esta opción modelo sólido se puede seleccionar haciendo clic en la opción “robusto modelo de restricción” en el “Seleccionar robusta inversión” cuadro de diálogo. Hay un factor de

corte que controla el grado en el que se utiliza este modelo robusto restricción. Si se utiliza un valor grande, por ejemplo 1,0, el resultado es esencialmente el de la método de inversión mínimos cuadrados suavidad limitados convencional. Si se utiliza un valor muy pequeño, por ejemplo 0.001, el resultado está cerca de la verdadera método de inversión restringida robusto.

Figura 29. Ejemplo de resultados de la inversión utilizando la suavidad-restricción y robustos restricciones modelo de inversión. (A) la resistividad aparente pseudosección (array Wenner) para un modelo de prueba sintético con un bloque de fallo (100 ohm.m) en el lado inferior izquierdo y un pequeño bloque rectangular (2 ohm.m) en el lado derecho con un medio circundante de 10 ohm.m. Los modelos de inversión producidos por (b) el de mínimos cuadrados convencionales método suavidad limitados y (c) el robusto método de inversión.

98

 Apéndice L: matrices no convencionales o generales Esta característica es para atender a disposiciones de electrodos que no están c omprendidos en los tipos habituales de la matriz de electrodos o disposiciones o formas inusuales de la realización de las encuestas. Probablemente hay un número infinito de posibles configuraciones de electrodos que están limitadas sólo por la imaginación del usuario, pero en la mayoría de los casos, es probable que sean ligeras variaciones de las matrices estándar. La Figura 30 muestra cuatro configuraciones no estándar posibles.

Figura 30. Algunas de las posibles matrices no convencionales. (A) fourelectrode no simétrica tipo Wenner-Schlumberger de matriz. (B) matriz dipolo-dipolo con dipolos de tamaño desigual. (C) Una configuración de electrodo posible pero probablemente no viable.

Una posibilidad es una variación no simétrica de una matriz simétrica, tal como el Wenner-Schlumberger o la matriz dipolo-dipolo (Figuras 30a y 30b). Tal configuración podría surgir de encuestas con medidores de resistividad multicanal donde las mediciones se hacen con varios pares de electrodos de potencial para un par común de electrodos de corriente.

Mientras que el programa le da al usuario una mayor libertad en la di sposición de electrodos, algunas matrices que son técnicamente posible, pero tendrían señales potenciales muy bajos deben ser evitados.

Para dar cabida a las diversas posibilidades, un formato de datos que es similar a la utilizada para las encuestas transversal del pozo (Apéndice I) se usa. La ubicación x, así como la elevación de todos los electrodos utilizados en una medición deben ser dados. El MIXED.DAT archivo es un arc hivo de datos de ejemplo, con un formato tal. Esto es e n realidad un conjunto de datos sintética con una mezcla de mediciones usando el Wenner

99

Schlumberger y dipolo-dipolo matrices. La parte inicial de este archivo de datos con comentarios sobre el formato es la siguiente. archivo MIXED.DAT

matriz mixta 1.0 11 0

comentarios | Título | Unidad de separación de los electrodos

| tipo Array (11 para array general)

| tipo Sub-array

Tipo de medición (0 = aplicación. Resistividad, 1 = resistencia)

0 407 1 0

| Encabezamiento

| 0 aquí para indicar valores de resistividad aparente | Número de puntos de datos | Indicador para el tipo de posición en x | No hay valores de IP, 1 si es IP presente

4 0,00 0,00 3,00 0,00 1,00 0,00 2,00 0,00 10.158 | El formato para cada punto de datos es 4 1,00 0,00 4,00 0 ,00 2,00 0,00 3,00 0,00 10.168 | Número de electrodos utilizados en

4 2,00 0,00 5,00 0,00 3,00 0,00 4,00 0,00 10.184 | medición, x y z -ubicación de 0.00 6.00 4 3.00 0.00 4.00 0.00 5.00 0.00 10.225 | electrodo C1, x y z-ubicación de 4 4,00 0,00 7,00 0,00 5,00 0,00 6,00 0,00 10.337 | electrodo C2, x y z-ubicación de 4 5,00 0,00 8,00 0,00 6,00 0,00 7,00 0,00 10.708 | electrodo P1, x y z-localización de 6.00 0.00 9.00 4 0.00 7.00 0.00 8.00 0.00 11.668 | electrodo P2, resistividad aparente 10.00 7.00 0.00 4 0.00 8.00 0.00 9.00 0.00 12.542 | o v alor de la resistencia.

El indicador de tipo de sub-array se utiliza cuando la configuración de electrodos sigue una de las matrices convencionales, por ejemplo la matriz Wenner-Schlumberger. A modo de ejemplo, el MIXEDWS.DAT archivo tiene los datos para sólo en una matriz WennerSchlumberger pero la información se organiza en un formato general matriz. En este caso el número de sub-array es 7. El RATHCMIX.DAT archivo tiene los datos para el montículo Rathcrogan (descrito anteriormente en la Figura 8) para la matriz Wenner pero en un formato de matriz general (por lo que tiene el número de sub-matriz de 1 ) . Para estos conjuntos de datos, los valles de resistividad aparente se pueden mostrar en la forma de un pseudosección. El MIXED.DAT archivo es un conjunto de datos combinados con una mezcla de mediciones con el Wenner-Schlumberger y la matriz dipolo-dipolo. En este caso, no hay ningún tipo de matriz consistente, por lo que el número de sub-array se da como 0.

En cuanto el indicador para el tipo de posición en x, hay dos val ores posibles. Un valor de 1 se utiliza si los valores de x en la localización son los verdaderos distancias horizontales, mientras que un valor de 2 se utiliza si es la distancia a lo largo de l a superficie del suelo. Esto es similar a la convención utilizada para los datos de topografía (Sección 7). En el archivo RATCMIX.DAT, la información que se incluye la topografía como la elevación para cada electrodo en cada línea de datos. En algunos casos, puede ser más conveniente para introducir la información topografía en una lista separada, al igual que para los datos en el formato normal. Un ejemplo con los datos en el formato general matriz, pero con la topografía en una lista separada después de las líneas de datos de resistividad aparente se da en la ITALY_CAVE.DAT archivo. Parte del archivo se muestra a continuación.

100

archivo ITALY_CAVE

comentarios

encuesta cueva Italia

| Título

1.00 11 3

| Unidad de separación de los electrodos

| tipo Array (11 para array general)

| tipo Sub-array (dipolo-dipolo)

Tipo de medición (0 = aplicación. Resistivity.1 = resistencia)

| Encabezamiento

0 584 2

| 0 aquí para indicar valores de resistividad aparente

1

| Nota presencia de datos IP

cargabilidad

| datos IP en cargabilidad

mV / V

| unidad de IP

0.16,0.5

| parámetros de medición IP

| Número de puntos de datos | Indicador para el tipo de posición en x (superficie | distancia)

4 2,00 0,00 0,00 0,00 4,00 0,00 6,00 0,00 126,16 0,23 | 2 primeras líneas de datos

4 2,00 0,00 0,00 0,00 6,00 0,00 8,00 0,00 169,15 -2,80 | Nota 0.0 para la elevación | Otras líneas de datos . . | 88.00 86.00 0.00 4 0.00 92.00 0.00 94.00 0.00 44.73 1.71 | 2 últimas líneas de datos

90.00 88.00 0.00 4 0.00 92.00 0.00 94.00 0.00 42.93 1.37 | Topografía en lista separada | Cabecera para indicar la información de la topografía 2

| está en una lista separada después de las líneas de datos. Formato

85 0,94.158 2,94.269 . . 94,101.948 95,102.062 1 0 0 0

| para los datos de topografía es el mismo que para los datos en | formato normal. 2 primeros puntos de topografía.

| | Otros puntos de datos topografía | | 2 últimos puntos de datos topografía

| | Indica primer electrodo coincide con el primer  | punto de topografía | Termina con algunas de cero

|

Tenga en cuenta la adición de una línea de cabecera “Topografía en lista separada” que indica la información de la

topografía es en una lista separada después de las líneas de datos de resistividad aparente. Tenga en cuenta también en las líneas de datos, la elevación para cada electrodo se da como 0 ,00. El programa calculará la elevación de los electrodos de la información de la topografía.

La medición puede ser dado en términos de resistencia o res istividad aparente. Si no hay topografía, es decir, el Z-ubicaciones son todos iguales a cero o la misma, la fórmula convencional puede ser usado para calcular la resistividad aparente. El factor geométrico k está dada por la siguiente ecuación

k = 2 pi / (1 / r1 - 1 / r2 - 1 / r3 + 1 / r4)

101

donde r1 es la distancia entre los electrodos C1 y P1, R2 es la distancia entre el C1 y electrodos P2, r3 es la distancia entre el C2 y los electrodos P1 y r4 es la distancia entre el C2 y electrodos P2. En el caso cuando no hay topografía, la distancia entre dos electrodos es simplemente la diferencia en el X-ubicaciones de los dos electrodos.

En el caso en que la topografía está presente, no existe una convención ampliamente aceptada. Para ser consistente con el formato utilizado anteriormente, dos convenciones se utilizan diferentes dependiendo del tipo de posición en x dada en el archivo de datos. Si los X-ubicaciones se dan como distancias a lo largo de la superficie (es decir, no verdaderos distancias horizontales), la distancia entre dos electrodos se calcula utilizando la diferencia en los valores de las únicas x-lugares.

Si los X-lugares son verdaderos distancias horizontales, la distancia real entre dos electrodos se utiliza en el cálculo del factor geométrico. Por ejemplo, si los electrodos C1 y P1 están situados en ( x 1, z 1) y ( x 2, z 2) respectivamente, entonces

r1 = sqrt ( dx dónde

2 + dz 2)

dx = x 1 - x 2 dz = z 1 - z 2

En el caso en que la superficie del suelo tiene una pendiente constante, ambas convenciones darán el mismo factor geométrico. Tenga en cuenta que en el formato de datos de la matriz general, el z-ubicación es la elevación del electrodo, que es positivo en la dirección hacia arriba. Esto, desafortunadamente, es diferente del formato para los datos de la pared del pozo establecidos donde los valores z son positivos en la dirección hacia abajo.

La figura 31 muestra una disposición interesante donde los electrodos están dispuestos de tal manera que la horizontal verdadera distancias entre electrodos adyacentes son constantes. Cuando la pendiente de la superficie del suelo no es constante, esto resulta en una disposición en la que la distancia entre electrodos adyacentes, medida a lo largo de la superficie del suelo no es constante sino que depende de la pendiente. En las zonas donde la pendiente es más pronunciada, la distancia entre electrodos adyacentes a lo largo de la superficie del suelo es mayor. Para dar cabida a una disposición de este tipo, el formato general array se puede utilizar. En este caso, el indicador para el tipo de posición en x debe ser de 1 ya que las verdaderas distancias horizontales se dan en los datos. Muy a menudo, la encuesta s e lleva a cabo c on una de las matrices convencionales, tales como la WennerSchlumberger, por lo que el número de sub-array apropiado también debe ser dada.

102

Figura 31. Disposición de los electrodos a lo largo de una línea de encuesta de tal manera que la verdadera distancia horizontal entre electrodos adyacentes es constante independientemente de la pendiente de la superficie del suelo. Tenga en cuenta que la distancia entre electrodos adyacentes a lo largo de la superficie del suelo es mayor en la parte más empinada de la pendiente.

En algunos casos, el conjunto de datos 2-D se genera a partir de una serie de colineales 1-D líneas de sondeo. En tal situación, la distancia entre electrodos adyacentes no es constante, ya sea medida a lo largo de la superficie o en el verdadero sentido horizontal. Para un caso tal, el formato de matriz general (con un tipo de sub-array de 0) se puede utilizar.

Para conjuntos de datos con un tipo sub-array de 0, los datos no se pueden mostrar en la forma de un pseudosección, ni puede editado con la opción "exterminar puntos de datos malo". Para eliminar puntos de datos malos del conjunto de datos, es necesario primero realizar una inversión de todo el conjunto de datos. Sin embargo, es necesario utilizar la "restricción de datos robusta" (ver Apéndice K) para asegurar que los puntos de datos malos no tienen un gran e fecto en los resultados de la inversión. A continuación, vaya a la ventana "Display" y leer en el archivo con los resultados de la inversión. Seleccione la opción "estadísticas de error RMS" sub-opción en el menú "Editar datos". Esto mostrará un histograma donde los puntos de datos se agrupan de acuerdo a la diferencia entre los valores de r esistividad aparente medidos y calculados. Esto le permite eliminar los puntos de datos donde una gran diferencia se produce, por ejemplo, más de 100 por ciento. Después de la eliminación de los puntos de datos más ruidosos, guardar los datos recortados establecidos en un archivo y llevar a cabo la inversión de nuevo con el conjunto de datos recortado.

La sub-opción "Generar bloques de modelo" en la opción "Inversion" en el menú principal asigna la disposición de los bloques en el modelo utilizando los valores de sensibilidad como guía. La figura 32 muestra los valores de sensibilidad de los datos MIXEDWS.DAT establecidos utilizando el modelo de defecto extendido, y para la opción de "generar bloques modelo". Tenga en cuenta que esta opción aumenta la anchura de los bloques cerca de las partes en las partes más profundas de l a sección del modelo con el fin de incrementar los valores de sensibilidad de los bloques. Esta opción también se puede utilizar para

103

normal de datos establece con las matrices estándar. Por defecto, el p rograma utiliza un algoritmo heurístico en parte sobre la base de la posición de los puntos de datos para generar el tamaño y la posición de los bloques de modelo. El "Generar bloques modelo" utiliza un enfoque más cuantitativo basado los valores de se nsibilidad de los bloques de modelo. Sin embargo, en la práctica, el algoritmo más simple y más rápido heurística da unos resultados bastante buenos para la mayoría de los conjuntos de datos.

Figura 32. La disposición de los bloques de modelo utilizado por el programa para la inversión del conjunto de datos MIXEDWS.DAT usando (a) el modelo por defecto extendido y (b) un modelo generado de modo que ninguna de las cuadrículas de la no tiene una sensibilidad relativa valor que es demasiado pequeño.

104

 Apéndice M: tipos de matriz en RES2DINV Una lista de las matrices apoyadas por el p rograma RES2DINV así como los códigos numéricos se dan a continuación. nombre de la matriz

número de código

Wenner (alfa) Polo-Polo

Inline dipolo-dipolo Wenner (beta) Wenner (gamma) Polo-dipolo Wenner-Schlumberger 

Ecuatorial dipolo-dipolo / Matriz general no convencional

1 2 3 4 5 6 7 8 11

encuesta Cross-perforación (valores de resistividad aparente) 12 encuesta Cross-perforación (valores de resistencia)

13

matriz de gradiente (sólo se utiliza como el número de sub-matriz con datos en el formato general array)

15

Tenga en cuenta que la matriz Wenner tiene tres variaciones diferentes (Figura 33). La matriz "normal" Wenner es en realidad la matriz alfa Wenner.

Figura 33. Disposición de los electrodos para las alfa, beta y gamma arrays W enner.

105

La matriz beta Wenner es un caso especial de l a matriz dipolo-dipolo donde el factor "n" es siempre 1. El programa RES2DINV convertirá automáticamente un archivo de datos de matriz beta W enner en una matriz de datos dipolo-dipolo establecidos. La matriz de gamma Wenner no se utiliza con frecuencia, pero tiene un apoyo directo en RES2DINV. El GAMMA.DAT archivo es un a rchivo de datos de ejemplo para esta matriz. Una historia clínica de un estudio de campo en África usando esta matriz se da en la sección 18 de este manual.

106

 Apéndice N - Guardar los resultados en formato SURFER Esta opción permite al usuario guardar el pseudosecciones de resistividad aparente y secciones modelo en el formato XYZ utilizado por el programa trazado gráficos 2-D Surfer. Después de leer en un fichero que contiene el INV en resultados de la inversión en la ventana 'Display', se recomienda que se visualizan los pseudosecciones al menos una vez antes de utilizar la opción de guardar los resultados en formato Surfer. Si el archivo contiene la inversión los datos topográficos,

como

el ejemplo

RATCHRO.INV archivo proporcionado con este paquete de software, tendrá que utilizar seleccione la opción 'Incluir en la topografía exhibición modelo' para mostrar la sección del modelo con la topografía. Después de eso, puede hacer clic en 'Guardar los datos SURFER formato en' opción en el menú Archivo. Para el archivo RATCHRO.INV, se mostrará el siguiente cuadro de diálogo.

En este cuadro de diálogo, puede seleccionar el número de i teración para el que desea guardar los resultados de la inversión. Además, puede optar por guardar los valores de resistividad directamente, o el logaritmo de los valores de resistividad que podría ser más apropiado si los valores de resistividad tienen una amplia gama. SURFER trazará las secciones como una sección rectangular. En muchos casos, especialmente una sección de modelo con la topografía, esto podría no se r apropiada. Para eliminar las secciones fuera de la sección del modelo, SURFER proporciona dos métodos de supresión de parte de la sección contorneada mostrada. El valor por defecto y normalmente mejor método es el uso de la Base

107

función de mapa de superposición en Surfer. En la mitad inferior del cuadro de diálogo, se muestran todos los diversos posibles secciones de datos o modelo. Para un archivo de inv ersión en particular, no todos

el

las posibles secciones están disponibles. Por ejemplo, RATHCRO.INV archivo no tiene datos IP, por lo que los círculos a la derecha del cuadro de diálogo no se comprueban. A modo de ejercicio, haga clic en el 'modelo de resistividad con la topografía' opción para el archivo RATHCRO.INV. Después de hacer clic en el botón Aceptar, el programa guardará e ntonces las

distintas secciones. Para una sección particular, el programa guardará 5 archivos diferentes que se utilizan para diferentes funciones en Surfer. Como un ejemplo, seguir los siguientes pasos para trazar la sección del modelo de resistividad (logaritmo) con topografía para el archivo RATHCRO.INV. Los diferentes archivos se proporcionan en el paquete de software RES2DINV.

1). Crear archivo de contorno. Haga clic en 'Guía' en la barra de menú superior, y luego 'Datos' que debe leer en el archivo RATCHRO_TOPRESLOG.DAT. Después de leer en el archivo de datos, SURFER mostrará el cuadro de diálogo 'Dispersos datos de interpolación'. En la sección 'Línea de cuadrícula Geometría', se muestra el número de líneas que SURFER utiliza para interpolar los valores de datos en una rejill a rectangular. Para este conjunto de datos, los valores predeterminados son probablemente 50 líneas en la dirección X y 6 líneas en la dirección y. En la mayoría de los casos, se debe aumentar los valores por defecto en un factor de 4 a 10 veces para mejorar la precisión de la interpolación. Para este ejemplo, trate de cambiar el número de líneas en la direcciones X e Y a 250 y 30 respectivamente (es decir, aumentar los valores por defecto en un factor de 5). El programa SURFER creará un archivo RATCHRO_TOPRESLOG.GRD correspondiente. A continuación haga clic en 'Mapa' en la barra de menú superior y luego 'Contour'. En la opción 'Contour', leído en el archivo RATCHRO_TOPRESLOG.GRD. A continuación se mostrará el cuadro de diálogo 'Contour Mapa'. En este cuadro de diálogo, haga clic en la opción 'Rellenar Contornos'. A continuación, haga clic en el botón 'Load', y después lee en el archivo RATCHRO_TOPRESLOG.LVL para establecer los valores de contorno y colores. Si desea visualizar la e scala de colores y suavizar los contornos, haga clic en el 'Escala de Color' y las opciones de los contornos suaves. A continuación haga clic en OK. y luego leer en el archivo RATCHRO_TOPRESLOG.LVL para establecer los valores de c ontorno y colores. Si de sea visualizar la escala de colores y suavizar los contornos, haga clic en el 'Escala de Color' y las opciones de los contornos suaves. A continuación haga clic en OK. y luego leer en el archivo RATCHRO_TOPRESLOG.LVL para establecer los valores de contorno y colores. Si desea visualizar la escala de colores y suavizar los contornos, haga clic en el 'Escala de Color' y las opciones de los contornos suaves. A 2). Lo siguiente que necesita para enmascarar las secciones de l a representación del contorno. En el menú 'Mapa', haga clic en 'Cargar Mapa Base'. Elija el RATCHRO_TOPRESLOG.BLN que el archivo de mapa base. En los GS Tapas de opciones de importación ', basta con hacer clic en OK. No utilice las "áreas de Curves opción. A continuación, pulse F2 y, a continuación, haga clic en la opción 'Overlap Mapas' en el menú 'Mapa'. Esto combinar el mapa de contorno y la base de supresión mapa en un único mapa. Seleccione la superposición haciendo clic cerca de una de las esquinas del rectángulo que contiene el mapa de contorno, y luego seleccione el comando 'Modificar superposiciones' en el menú 'Mapa'. En el cuadro de diálogo 'Editar Componentes', haga clic en la 'base' y haga clic en 'Mover al principio' botón. A continuación mantenga pulsada la tecla CTRL del teclado y haga clic en l a superposición de hasta

el mapa base es seleccionada (como se muestra por 'Mapa: Base' el indicador de estado cerca de la esquina inferior izquierda de la pantalla). A continuación haga clic en el menú 'Dibujar' y, a continuación, elija la opción '' atributos de relleno. Seleccionar el color blanco o blanco como color para el relleno

108

atributos. A continuación, haga clic en la 'Línea de Atributos', y también seleccionar el color blanco para que las líneas de superposición no son visibles en la pantalla. 3). Haga clic en el comando 'Post' en el menú 'Mapa', y después lee en el archivo RATCHRO_TOPRESLOG_POST.DAT (que contiene la ubicación de los electrodos) en el cuadro de diálogo 'Open Data'. Después de leer en este archivo, se mostrará el cuadro de diálogo ' Post Mapa'. En este cuadro de diálogo, seleccione la forma y el tamaño de símbolo y otras opciones q ue desee utilizar y, a continuación, haga clic en Aceptar. Después de eso, presione F2 y, a continuación, haga clic en el comando 'Overlay Mapas' en el menú 'Mapa'. A continuación, haga clic en la figura y luego el comando 'Modificar superposiciones' en el menú 'Mapa'. Haga clic en 'Publicar' y el 'Mover al principio', y los electrodos se dibujarán en la parte superior de la sección del modelo.

Después de dibujar la sección en la ventana Surfer, se pueden utilizar otras opciones dentro de SURFER para mejorar el modelo de la sección, como cambiar el tipo de letra, etiquetas, títulos, etc., antes de imprimir los resultados o exportarlo a otro formato. SURFER es compatible con una gran vari edad de formatos, incluyendo formatos gráficos de mapas de bits, como PCX, BMP, GIF, etc .; y así como los formatos basados en vectores tales como el formato DXF de AutoCAD.

Para las secciones de modelo sin la topografía, la etapa (3), donde el mapa de datos posterior se lee no es probable que sea necesario. Sin embargo, es posible que tenga que ajustar las marcas de graduación en la escala vertical en el lado izquierdo de la gráfica de contorno. Las profundidades se dan como valores negativos de modo que SURFER trazará la figura en la orientación correcta con la superficie del suelo en la parte superior. A modo de ejemplo, di bujar la sección de modelo para el ejemplo Sting cueva utilizando la CAVE_MODRESLOG.DAT, CAVE_MODRESLOG.LVL y archivos CAVE_MODRESLOG.BLN proporcionado. Después de dibujar la sección del modelo, es probable que sólo se verá el 0.0 marca en el eje vertical izquierdo de la escala de profundidad. Haga clic en el eje de la izquierda con el botón derecho del ratón y el cuadro de diálogo 'Izquierda Eje' se mostrará. A continuación haga clic en el botón de 'ampliación' en el cuadro de diálogo. En el cuadro de diálogo 'Eje de escala', cambie el valor 'Major Intervalo' a 10 y el 'primer gran Tick' a -40. Haga clic en Aceptar para guardar los valores y el modelo de la sección izquierda del eje y luego debe mostrar marcas a intervalos de profundidad de 10.

Habrá ligeras diferencias en las secciones de contorno dibujado por surfista y RES2DINV debido a las diferencias en los métodos de interpolación utilizados. SURFER no utiliza los puntos de datos originales, sino una cuadrícula rectangular de puntos interpolados para dibujar los contornos. Las diferencias probablemente serán más evidentes en zonas con cambios muy rápidos en l a resistividad. Una forma de reducir las diferencias es aumentar el número de líneas de cuadrícula utilizadas por surfista para la red de interpolación.

109

O Anexo: Modelo refinamiento El programa RES2DINV por defecto utiliza un modelo en el que la anchura de las c élulas modelo interior es la misma que la distancia entre electrodos unidad (por ejemplo, como en la Figura 2). Esto funciona bien en la mayoría de los casos. En algunas situaciones, particularmente donde hay grandes variaciones de resistividad cerca de la superficie del suelo, los mejores resultados se pueden obtener mediante el uso de células modelo más estrechos. Hay dos maneras posibles de reducir la anchura de las células modelo.

La primera es mediante el uso de la opción “Usar el refinamiento del modelo” en el menú 'inversión”. Al hacer clic en

esta opción, se mostrará el siguiente cuadro de diálogo.

Esto le permite elegir las células modelo con anchos d e la mitad de la distancia entre los electrodos unidad. En casi todos los casos, esto da los resultados óptimos. Después de seleccionar esta opción, debe leer en el archivo de datos. Entonces, el programa reducirá automáticamente el electrodo unidad espaciadora se utiliza por medio de la que figura en el archivo de datos. El segundo método consiste en modificar el archivo de datos di rectamente utilizando un editor de texto. El espaciamiento de los electrodos unidad se da en la segunda línea del archivo de datos de resistividad aparente. Por ejemplo, las primeras líneas del archivo PIPESCHL.DAT se reproducen a continuación.

encuesta tubería subterránea

1,00 7 173 1 0

110

1,50 1,00 1 641.1633 En la segunda línea de datos, la separación entre electrodos unidad se da como 1,0 metros que es el espacio real utilizado en esta encuesta. Cuando el RES2DINV lee en este archivo de datos, se fijará la anchura de las celdas modelo a 1,0 metros. El PIPESCHL_HALF.DAT archivo tiene los mismos datos, excepto que la separación entre electrodos unidad se define ahora como 0,5 metros. La parte superior de este archivo se muestra a continuación.

encuesta tubo de metro - media distancia entre electrodos real

0,5 7 173 1 0

1,50 1,00 1 641.1633 Cuando el programa RES2DINV lee este archivo, fijará las células a un modelo 0,5 metros de ancho. La figura 34 muestra los modelos de inversión para este conjunto de datos utilizando células con anchuras de 1,0 y 0,5 metros. El modelo con el ancho de celda 0,5 metros es capaz de modelar con mayor precisión las altas anomalías de resistividad cerca de l a superficie entre las marcas 3 y 12 metros por encima de la tubería de baja resistividad, y por lo tanto tiene un error de resistividad RMS ligeramente inferior aparente. La matriz Wenner-Schlumberger se utilizó para este estudio.

En la mayoría de los casos con grandes variaciones de resistividad cerca de la superficie, usando un modelo con una anchura de la celda de la mitad de la distancia real electrodo de unidad es suficiente. En algunos casos, el uso de un ancho de celda más estrecha de un tercio o un cuarto de la unidad podría estar justificado, en particular si la encuesta se llevó a cabo usando un polo-dipolo o matriz dipolo-dipolo con muy altos factores 'N'. Si la anchura de la celda es un cuarto de la verdadera separación entre electrodos unidad, o menos, se pueden producir oscilaciones en los valores de modelo de resistividad en las primeras capas. Así, el uso de un ancho de celda de menos de una cuarta parte de la verdadera separación entre electrodos unidad no es aconsejable.

Se recomienda que consulte la sección 4.7j de las notas libres “Tutorial: 2D y encuestas de imagen

eléctrico 3-D” que se pueden descargar desde el sitio web www.geoelectrical.com. Esta sección tiene una discusión sobre el efecto del uso de altos fact ores de 'n' con la matriz de polo-dipolo de la sensibilidad de la matriz a anomalías cerca de la superficie.

111

Figura 34. (a) La aparente pseudosección resistividad para el conjunto de datos PIPESCHL.DAT. Los modelos de inversión utilizando (B) células con una anchura de 1,0 metros que es el mismo que el electrodo unidad real, y (c) el uso de células estrechas con una anchura de 0,5 metros. El tubo está marcada por la anomalía de baja resistividad por debajo de los 12 metros. Tenga en cuenta las anomalías de alta resistividad cerca de la superficie entre las marcas 3 y 12 metros. La matriz Wenner-Schlumberger se utilizó para este estudio.

112

 Apéndice P: Combinación de archivos 2-D

El programa RES2DINV tiene 2 opciones para combinar archivos de datos 2-D. En una opción, destinada a una serie de líneas paralelas de la encuesta, se puede combinar una serie de ficheros de datos 2-D en un solo archivo de datos 3-D. Esta es la opción “Intercalar  datos en formato RES3DINV” en el menú “Archivo”. En la segunda opción, destinada a varias jaulas-largo de las líneas de la encuesta, puede combinarlos en un solo archivo de datos 2-D. Esta es la opción “Concatenar  archivos de datos 2-D” en el menú “Archivo”.

Cotejar los datos en formato RES 3DINV

Idealmente, una encuesta 3-D debe llevarse a cabo usando una rejilla rectangular de electrodos con las mediciones en diferentes direcciones (Loke y Barker 1996b). Sin embargo, por diversas razones, a veces esto no es posible. En algunos estudios, se utilizan una serie de líneas de levantamiento 2-D paralelos, con, posiblemente, algunas líneas de unión en una dirección perpendicular. Esta opción permite al usuario al usuario combinar las líneas 2-D con los datos en el formato RES2DINV en un único archivo de datos en el formato utilizado por el programa RES3DINV. Si bien el programa RES3DINV llevará a un verdadero 3-D de inversión (en que se permite que los valores de resistividad a variar en los 3 direcciones simultáneamente durante el proceso de inversión), si el conjunto de datos contiene suficiente información 3-D con el fin de producir una razonablemente exacta modelo 3-D es un asunto diferente.

113

Un archivo de comandos se utiliza para combinar los archivos 2-D. El COLLATE_2D_3D.TXT archivo es un archivo de script de ejemplo. El contenido de este archivo junto con una descripción del formato utilizado es la siguiente.

archivo COLLATE_2D_3D.TXT

| comentarios

La conversión de archivos de datos RES2DINV

| Título

Número de archivos para clasificar

| Cabecera para no. de archivos de datos

3

| Número de archivos de datos

Presentar 1 Parámetros

| Encabezado por el primer archivo

Nombre del archivo de datos en formato RES2DINV | Encabezado por el nombre del archivo d: \ test \ FILE2D_1.DAT

| Nombre completo además de la ruta del archivo

X e Y ubicación del primer electrodo a lo largo de esta línea | Encabezamiento

0.0,0.0

| Coordenadas del primer electrodo

dirección de la línea (0 = X, 1 = Y)

| Encabezamiento

0

| Número dirección de la línea que especifica

signo Line (0 = positivo, 1 = negativo)

| Encabezamiento

0

| Especificar si el electrodo | coordina aumento o disminución | una larga linea

2 presentar parámetros

| Mismo conjunto de parámetros

Nombre del archivo de datos en formato RES2DINV

| para el segundo archivo

d: \ test \ FILE2D_2.DAT

|

X e Y ubicación del primer electrodo a lo largo de esta línea |

0,0, -0,5

|

dirección de la línea (0 = X, 1 = Y)

|

0

|

signo Line (0 = positivo, 1 = negativo)

|

0

|

.... ....

| Lo mismo para el tercer

Nombre de archivo de salida en formato RES3DINV

| Encabezamiento

d: \ test \ FILE_3D.dat

| Nombre de archivo de datos 3-D

Fin del documento

| Cabecera para el final del archivo

| archivo

El formato utilizado permite la posibilidad de que el estudio tenía líneas en 2 direcciones perpendiculares. Una línea encuesta 2-D tiene solamente una dirección horizontal, y la coordenada de un electrodo a lo largo de la línea de encuesta se da como la distancia a l o largo de la línea, a partir del primer electrodo. Para determinar la coordenada del electrodo en el plano XY para una encuesta de 3- D, necesitamos determinar la coordenada del primer electrodo en la línea de encuesta 2-D dentro de la cuadrícula encuesta XY. Esto da un punto de referencia para el que se pueden calcular las coordenadas de los otros el ectrodos a lo largo de la misma línea de encuesta.

114

En segundo lugar, es necesario conocer la orientación de la línea, es decir, si es a lo largo de la dirección X o dirección Y en la cuadrícula XY. Tenga en cuenta que las líneas diagonales no están permitidos.

En una línea de levantamiento 2-D, se supone siempre que la coordenada del electrodo siempre aumenta a partir del primer electrodo (es decir, de izquierda a derecha en una pseudosección típico). Sin embargo, la dirección de aumento del valor a lo largo de la línea de enc uesta 2-D de coordenadas no podría siempre ser el mismo que el X positivo o dirección Y en la c uadrícula encuesta XY 3-D. Los parámetros 'Línea signo' permite la posibilidad de que la dirección de la línea de levantamiento 2-D es opuesta a la dirección de la X positivo o eje Y en la cuadrícula XY.

Si bien es posible combinar una serie de líneas 2-D en un conjunto de datos 3-D, que puede no ser la pena hacerlo. En primer lugar, se recomienda que debe haber al menos 5 líneas paralelas.  Además, la separación entre las líneas no debe ser más del doble de la distancia entre electrodos unidad a lo largo de las líneas. Para más detalles, consulte el capítulo 8 de los libres “Tutorial: encuestas de imagen eléctrica 2-D y 3-D” por MHLoke. Si no tiene una copia, que puede ser descargado desde el sitio web www.geoelectrical.com.

Concatenar los datos en formato RES2DINV Esta opción permite al usuario combinar despliegue a lo largo de los archivos de datos medidos a lo largo de la misma línea en un solo archivo de datos. Todos los electrodos utilizados en las líneas de levantamiento línea a lo largo de una sola línea recta. El archivo CONCATENATE_2D.TXT un archivo de script de ejemplo. El contenido de este archivo junto con una descripción del formato utilizado es la si guiente.

archivo CONCATENATE_2D.TXT

| comentarios

Concatenación de varios archivos de datos en RES2DINV | Puesto Número de archivos para concatenar

| Encabezamiento

3

| Número de archivos de datos

Presentar 1 Parámetros

| Encabezamiento

Nombre del archivo de datos en formato RES2DINV

| Encabezamiento

c: \ data \ file1.dat

| Nombre completo además de la ruta del archivo

X ubicación del primer electrodo a lo largo de esta línea

| Encabezamiento

0.0

| Coordenadas del primer electrodo

signo Line (0 = positivo, 1 = negativo)

| Encabezamiento

0

| Especificar si el electrodo | coordina aumento o disminución | una larga linea

2 presentar parámetros

| Información similar para los días 2 y 3

Nombre del archivo de datos en formato RES2DINV

| archivos siguen

c: \ data \ file2.DAT

|

...... .. ...... ..

|

Nombre de archivo de salida en formato RES2DINV

| Encabezamiento

|

115

c: \ data \ file_123.DAT

| Nombre de archivo para almacenar datos combinados

Fin del documento

|

En algunos sistemas de electrodos múltiples metro de la resistencia, cada nuevo registro se inicia con una ubicación cero para los primeros electrodos de datos. Con el fin de combinar diferentes archivos de datos, los X-ubicación de los puntos de datos en algunos de los archivos tendrán que cambió de manera que todos ellos tienen un origen común. En el formato de archivo de secuencia de comandos anterior, esto se hace mediante la introducción de la coordenada del primer electrodo para la línea de encuesta. En la mayoría de los casos, la dirección de las líneas son las mismas, por ejemplo los X-ubicaciones aumenta desde la izquierda a la derecha. Sin embargo, se proporciona una opción donde se midieron las líneas de estudio en diferentes direcciones.

116

 Apéndice Q: detección de bordes Bedrock En esta opción, el programa intenta elegir la interfaz de lecho de roca de la sección del modelo y guarda el x y y coordenadas en un archivo ASCII. Este método es más exitoso en el lecho de roca tiene una interfaz nítida con el suelo suprayacente. Una situación típica se muestra en la Figura 35. Este estudio se realizó para mapear el un lecho de roca diorita con una intención de extender una cantera existente (Turnbull 1986). La sobrecarga tiene resistividades de menos de 100 ohm.m mientras que la diorita tiene una resistividad de hasta varios cientos de ohm.m.

La función de detección de bordes lecho de roca está disponible d entro de la ventana 'Display', como se muestra a continuación.

 Al hacer clic en esta opción de menú, la sección del modelo primero, se tendrá, a continuación, se muestra el siguiente cuadro de diálogo. Al hacer clic en el cuadro de diálogo 'Detección de bordes' se muestra.

Hay 3 parámetros principales que necesitan que el usuario seleccione. La primera consiste en establecer los límites superior e inferior para la ubicación de la frontera. Por defecto, el programa seleccionará a partir de 0 m. (En la superficie del suelo) hasta que la medida máxima profundidad del modelo de inversión. En la mayoría de los casos, el límite superior se cambia comúnmente, en particular cuando hay gra ndes anomalías de resistividad

117

cerca de la superficie (por ejemplo, cantos rodados). En este ejemplo particular, la profundidad mínima se fija en 15 m. para la interfaz de la roca madre. El segundo parámetro importante es un valor mínimo para el lecho de roca. En este ejemplo, se eligió un valor de 150 ohm.m que probablemente corresponde a degradado parcialmente diorita. El tercer parámetro es el cambio máximo en la profundidad de lecho de roca entre electrodos adyacentes. Esto es para evitar los picos en el perfil de la roca madre.

En la mayoría de los casos, los valores predeterminados para los demás parámetros deberían funcionar bien. Por defecto, se asume que la base tiene una resistividad superior a la sobrecarga. Sin embargo, una disposición está hecha para los casos en que el lecho de roca o de una capa subsuperficial de interés tiene una resistividad inferior.

Para la mayoría de las situaciones, la opción de inversión sólido modelo se debe utilizar (Anexo K). Esto produce normalmente un límite de la roca madre más agudo en el

118

modelo de inversión. Además, se recomienda un modelo en el que las células tienen anchuras de la mitad de la distancia entre electrodos unidad real (Apéndice O). Por ejemplo la encuesta Yeguada se llevó a cabo con una separación de electrodos unidad de 25 m., Pero el modelo de inversión a continuación utiliza anchuras de células más pequeñas.

Figura 35. pseudosección resistividad aparente de la encuesta Haras junto con el modelo de inversión. También se muestra el perfil de la roca madre interpretado. Tenga en cuenta que el lecho de roca en el extremo derecho de la línea está por debajo de la profundidad máxima de la i nvestigación de la línea de levantamiento. levantamiento.

119

 Apéndice I: Time-lapse resistividad inversión invers ión En el estudio de los cambios de la resistividad del subsuelo con el tiempo, las encuestas de imagen bidimensional bidimensional de resistividad se repiten a menudo por la misma línea en diferentes momentos. momentos. Tales estudios incluyen el flujo de agua a través de la zona vadosa, los cambios en la tabla de agua debido a la extracción de agua (Barker y Moore 1998), el flujo de contaminantes químicos y fugas de las presas. Para los experimentos con lapso de tiempo, los cambios en la resistividad del subsuelo se estiman mediante el uso de los cambios en las mediciones de resistividad aparente. Por lo tanto resultados precisos sólo se pueden obtener si la resistividad aparente propios valores son suficientemente precisos. Es mejor evitar el uso del dipolo-dipolo y arrays polo-dipolo en tales experimentos si se usan separaciones más grandes es debido a la señal de relativamente baja a ruido. Puede ser que sea más seguro usar la matriz Wenner o preferiblemente la matriz Wenner-Schlumberger. De lo contrario, los r esultados extraños podrían obtenerse.

El TIMELAPS.DAT TIMELAPS.DAT archivo es un archivo de datos de ejemplo con 3 series de datos de resistividad aparente. La parte superior del archivo junto con comentarios sobre el formato de datos es el siguiente.

archivo TIMELAPS.DAT

comentarios

Ejemplo de datos de lapso de tiempo

| Título

1.00 7

| Unidad de separación de los electrodos

350 1 0

| Número de puntos de referencia

datos de la secuencia Time

| Cabecera para indicar los datos de lapso de tiempo

Número de secciones de tiempo

| Encabezamiento

3

| Número de serie de tiempo

unidad de tiempo

| Encabezamiento

horas

| Unidad utilizada en la medición del tiempo

intervalo de la sección segunda vez

| Encabezamiento

1.0

| diferencia horaria de conjunto de datos por segunda vez

intervalo de la sección segunda vez

| Encabezamiento

2.0

| diferencia horaria de conjunto de datos tercera vez

| tipo de matriz, Wenner-Schlumberger Wenner-Schlumberger en este | ejemplo ejemplo

| se da punto medio de x en la localización de la matriz matriz | No hay datos IP

1.50 1.00 1.0 11.868 11.868 11.868 | El formato para cada punto de referencia es la siguiente: 2,50 1,00 1,0 11.867 11.867 11.867 | posición en x, " un" separación de los electrodos, " norte" factor, 3.50 1.00 1.0 11,863 11,863 11,863 | resistividad aparente para el primer conjunto de datos, aparente

4.50 1.00 1.0 11,859 11,860 11,860 | resistividad de segundo conjunto de datos, 5.50 1.00 1.0 11,855 11,856 11,856 | resistividad aparente para tercer conjunto de datos.

6.50 1.00 1.0 11,851 11,851 11,851 | El formato de datos es en gran medida similar a los datos normales establecidos con sólo un conjunto de lecturas. Para un archivo de datos de lapso de tiempo, la diferencia principal es la sección extra sólo

120

antes de que el bloque de datos con los valores de resistividad aparente. Aquí el encabezado "datos de la secuencia de tiempo" se utiliza para indicar la presencia de los conjuntos de datos adicionales correspondientes a las mediciones de tiempo posteriores. La siguiente parte de esta sección adicional contiene información sobre el número de mediciones de series de tiempo (3 en este ejemplo). El número máximo de conjuntos de datos de series de tiempo permitido es 12 si el equipo tiene al menos 256 MB de RAM, y 6 si la memoria disponible es menor. Para cada punto de datos, el formato para la parte inicial es idéntica a la utilizada para los conjuntos de datos normales. La única adición es que después de que el valor de resistividad aparente para el primer conjunto de datos, se dan los correspondientes valores de resistividad aparente para los posteriores conjuntos de datos de tiempo.

Este programa permite el uso de técnicas de inversión conjuntos utilizando diferentes tipos de restricciones modelo cruzadas. Cuando se selecciona el "método de inversión de lapso de ti empo Seleccionar" en el menú "Cambiar la configuración", se muestra el cuadro de diálogo seguimiento.

En el primer conjunto de opciones en el cuadro de diálogo, varias limitaciones posibles que guía el proceso de inversión para los conjuntos de datos en tiempo múltiples proporcionadas. Con frecuencia, los datos de encuestas de lapso de tiempo realizados en diferentes momentos se invierten de forma independiente. Este método corresponde a la opción 'No hay restricciones" en el cuadro de diálogo anterior. Los cambios en los valores de resistividad del subsuelo se determinan entonces mediante la comparación de los valores de resistividad modelo obtenidos a partir de las inversiones de un conjunto inicial de datos y l os posteriores conjuntos de datos de tiempo. En muchos casos, tales

121

un enfoque ha dado resultados satisfactorios. Sin embargo, en teoría, ya que la inversión de cada conjunto de datos se lleva a cabo de forma independiente, no hay garantía de que las diferencias en los valores de resistividad se deben únicamente a los cambios reales en la resistividad del subsuelo con el tiempo.

En este programa, una técnica de inversión de articulación que utiliza el modelo a partir de los datos iniciales establecidos para limitar la inversión de los posteriores conjuntos de datos de tiempo se utiliza para minimizar las posibles distorsiones (Loke 1999). El conjunto de datos modelo inicial de inversión se usa como un modelo de referencia en la inversión de los posteriores conjuntos de datos de lapso de tiempo. En la opción "Minimizar los cambios", de mínimos cuadrados constriñen para minimizar las diferencias en los valores de modelo de resistividad entre el modelo inicial y el modelo de lapso de tiempo se utiliza un amortiguado simple. Los "mínimos cuadrados suavidad limitan" opción utiliza mínimos cuadrados suavidad Restringir para asegurar que las diferencias en los valores de modelo de resistividad varían de una manera suave. Por último, en la opción "robusto suavidad de restricción",

En la segunda serie de opciones, el usuario puede elegir para llevar a c abo la inversión de todos l os conjuntos de datos de tiempo de forma simultánea o secuencial. En la opción de inversión simultánea, la inversión de la primera o de referencia conjunto de datos es seguido por la inversión de los posteriores conjuntos de datos de tiempo en cada iteración. Tenga en cuenta que incluso si el "no" li mitaciones opciones se elige anteriormente, el modelo obtenido en la iteración final para el primer conjunto de datos se sigue utili zando como modelo de partida para las posteriores series de datos de tiempo.

En la opción de inversión secuencial, la inversión completa del primer conjunto de datos se lleva a cabo primero durante aproximadamente 5 iteraciones (o el número seleccionado de iteraciones). El modelo para el primer conjunto de datos para la iteración final se utiliza entonces como el modelo de referencia para la inversión de los posteriores conjuntos de datos de tiempo. La inversión de los últimos conjuntos de datos de tiempo comienza sólo después de la finalización de la inversión del primer conjunto de datos. Esta inversión secuencial es útil cuando la sección de modelo tiene muy grandes contrastes de resistividad.

Bajo el 'Select tipo de datos de lapso de tiempo de' opción, el usuario puede seleccionar el uso de los valores aparentes de resistividad directamente (bajo 'mediciones directas') o la diferencia en los valores de resistividad aparente entre el tiempo más tarde y primer conjunto de datos (bajo ' Diferencia de datos ') como el parámetro de datos para utilizar utilizado para la inversión de los posteriores conjuntos de datos de tiempo. El uso de la 'diferencia de los datos de opción asegura que si no hay diferencia en los valores de resistividad aparente, no habrá ningún cambio en los valores de modelo.

Por defecto, el modelo obtenido para el primer conjunto de datos se utiliza como el modelo de referencia para todos los posteriores conjuntos de datos de tiempo, es decir, el mismo modelo de referencia se utiliza para

122

todos los otros conjuntos de datos. En el caso donde hay un gran cambio en la resistividad de los modelos de la primera a la última serie de datos, se proporciona una opción donde se utiliza el modelo para el conjunto de datos anterior como el modelo de referencia. En este caso, el modelo para el primer conjunto de datos se utiliza como el modelo de referencia para el segundo conjunto de datos, mientras que el modelo para el segundo conjunto de datos se utiliza como el modelo de referencia para el tercer conjunto de datos, y de manera similar para los otros conjuntos de datos .

La última opción, el peso restricción de tiempo, le da la importancia relativa dada a mantener los modelos de tiempo posteriores lo más similar posible al modelo de referencia. Si se utiliza un peso grande, por ejemplo 10.0, los modelos para los posteriores conjuntos de datos de tiempo se mantienen muy similar a la del modelo de referencia. Esto minimiza eficazmente los cambios a través de l os modelos de tiempo y reduce la posibilidad de ruido en los conjuntos de datos que causan cambios espurias entre los modelos. Sin embargo, esto es a expensas de aumentar el desajuste entre los datos, es decir, la diferencia de RMS entre los datos calculados y medidos. Por defecto, un valor de 3,0 se utiliza para el peso restricción de tiempo. Un valor de 1,0 da la misma importancia a la minimización de los cambios de horario en el modelo y reducir al mínimo el error de datos RMS.

Después de llevar a cabo la inversión del conjunto de datos de lapso de tiempo, por favor, cambie a la ventana “Display” para leer en los resultados en el archivo de la inversión. Se proporcionan varias opciones para mostrar los resultados de los modelos de lapso de tiempo. Las opciones diseñadas específicamente para los modelos de inversión de lapso de tiempo se dan cerca de la parte inferior de la lista de “secciones de visualización” en la ventana “Display”, como se

muestra a continuación.

123

La opción "Display porcentaje de cambio en la resistividad" mostrará el cambio en la resistividad modelo obtenido a partir de la inversión de un conjunto de datos más tarde tiempo en comparación con el modelo de referencia de la inversión de la primera serie de datos. La o pción de "relación de resistividad Display" mostrará la relación de las resistividades del modelo de tiempo más tarde en comparación con el modelo de referencia.

Después de mostrar las secciones con el cambio porcentual en los valores de resistividad modelo, se puede cambiar el conjunto de datos de lapso de tiempo seleccionado utilizando las teclas 'Fin' 'Inicio' y. Por defecto se mostrarán los resultados del último conjunto de datos de lapso de tiempo. Puede mostrar los resultados de un conjunto de datos anteriores ajusta pulsando la tecla 'Inicio'.

Justo debajo de las opciones de menú para mostrar la hora secciones, hay un "cambio en el factor de  Archie Ley" sub-opción. Ley de Archie que proporciona la relación entre la resistividad de una roca porosa y el factor de saturación de fluidos es aplicable para ciertos tipos de rocas y s edimentos, en particular los que tienen un bajo contenido de arcilla. En el caso, la conducción eléctrica s e supone que es a través de los fluidos de llenado de los p oros de la roca. Ley de Archie está dada por

=  uρn ρ

φ

- mw

dónde ρ es la resistividad de la roca, ρ w es resistividad del fluido, φ es la fracción de la roca llena con el fluido, mientras un y metro son dos parámetros empíricos (Keller y Frischknecht 1966). Para la mayoría de las rocas, un es de aproximadamente 1, mientras metro está aproximadamente 2. Bajo ciertas condiciones especiales, la ecuación anterior se puede usar para determinar el cambio en la saturación de fluido o resistividad del fluido con el tiempo.

Si inicialmente el material subsuperficial está saturado con agua, y el cambio en la resistividad es causado por la retirada de agua (es decir una reducción en la fracc ión por unidad de volumen de la roca que se llena de agua), el factor de desaturación (Keller y F rischknecht 1966 ) puede ser calculado a partir del cambio en la resistividad del subsuelo. Para calcular el factor de desaturación, el valor de la metro Se necesita parámetro en la ley de Archie. Por defecto se establece en 2.0 pero este valor puede ser cambiado por el usuario.

En el otro extremo, considere un caso en el que el factor de saturación de fluido no cambia, pero la resistividad de los cambios de fluido (por ejemplo, en experimentos con trazadores con una sal conductora debajo de la tabla de agua). La relación de la resistividad del fluido en el momento más tarde a la resistividad inicial puede calcularse directamente a partir de las resistividades modelo.

Hay que subrayar que las dos situaciones anteriores son casos ideales. Los resultados sólo son válidos si la ley de  Archie se mantiene para el medio bajo la superficie. En muchos casos,

124

la relación entre el contenido de resistividad y agua medio es mucho más complicado (Olivar et al. 1990). Como un ejemplo de un experimento de campo de las mediciones de lapso de tiempo, las Figuras 36 y 34 muestra los resultados de las pruebas en el flujo de agua a través de la zona no saturada realizado en Birmingham (Barker y Moore, 1998). Cuarenta mil litros de agua se vertieron sobre la superficie del suelo cerca de la marca de 24 metros a lo largo de la línea de encuesta (Figura 36) durante un período de 10 horas. El conjunto de datos inicial recogido antes de que el agua de riego se inició junto con se muestra el modelo de inversión correspondiente en las figuras 36a y 36b. El subsuelo se compone de arena y grava es altamente no homogénea. Como comparación, la Figura 36c muestra el modelo de inversión de los datos recogidos después de que se completó la irrigación de la superficie del suelo. La distribución del agua no es muy claro a partir de una comparación directa de los modelos de inversión por sí solos.

Figura 36. (a) La resistividad aparente y (b) una sección modelo de inversión de la encuesta realizada en el inicio del estudio infiltración Birmingham. Esto muestra los resultados de la serie de datos inicial que forma el modelo base en la inversión conjunta con los posteriores conjuntos de datos de tiempo. Como comparación, el modelo obtenido de la inversión del conjunto de datos recogidos después de 10 horas de irrigación se muestra en (c).

125

La distribución de agua se determina más fácilmente a partir de un gráfico de la variación porcentual de la resistividad del subsuelo en los modelos de inversión para los conjuntos de datos tomadas en diferentes momentos (Figura 37) en comparación con el modelo de conjunto de datos inicial. El conjunto de datos recogidos a las 5 horas después del inicio del bombeo muestra una reducción de la resistividad (de hasta más de 50 por ciento) cerca de la superficie del suelo en las proximidades de the24 marca metros. La zona de baja resistividad cerca de la superficie alcanza su máxima extensión y la amplitud máxima después de aproximadamente 10 horas cuando el bombeo se detiene (Figura 37b). Doce horas después de que el bombeo se detiene, el penacho de baja resistividad se ha extendido hacia abajo y ligeramente hacia afuera debido a la infiltración del agua a través de la zona no saturada. Hay una disminución en el máximo porcentaje de reducción en los valores de resistividad cerca de la superficie debido a la migración del agua de la zona cerca de la superficie. Este efecto de la propagación de la pluma se vuelve cada vez más pronunciada después de 24 horas, 36 horas y 52 horas (Figuras 37d, 37d, 37f), debido a una mayor migración del agua. Tenga en cuenta que el límite inferior de la zona con aproximadamente 20 por ciento de reducción en los valores de resistividad tiende a aplanarse a una profundidad de alrededor de 3 metros (Figura 37E), donde la columna de humo desde la superficie se encuentra con la tabla de agua.

Datos perdidos :- Si una medición faltaba de una serie de tiempo, introduzca 0,0 como el valor de la resistividad aparente para que series de tiempo particular. El programa supondrá que un valor cero es una indicación de un punto de datos que faltan.

126

Figura 37. Las secciones que muestran el cambio en los valores de resistividad del subsuelo con el tiempo obtenida a partir de la inversión de los conjuntos de datos recogidos durante la infiltración y de recuperación de las fases del estudio.

127

 Apéndice S: La incorporación de las estimaciones de ruido de datos

 Algunos tipos de instrumentos de resistividad topografía grabar una estimación del ruido aleatorio para cada medición. Esto se hace generalmente mediante la repetición de la misma medición varias veces. La estimación de error de datos se calcula entonces a partir de la desviación estándar en las lecturas. En realidad, esto da una medida de la repetibilidad en las lecturas, es decir, el ruido aleatorio. Tenga en cuenta que esta estimación de error d e datos no incluye fuentes de errores sistemáticos, tales como acoplamiento entre los cables.

Para incluir la estimación del error de datos en un archivo de datos, los valores de error se registran después de que el valor de la resistividad aparente. El MODEL35_N.DAT archivo es un ejemplo de un archivo de datos de una encuesta superficie con los valores de error. Parte del archivo se muestra a continuación con una explicación de los cambios necesarios.

MODEL35_N.DAT

comentarios

Modelo de Bloques

| Título

1.00

| Unidad de separación de los electrodos

7 230 1 0

| array Wenner-Schlumberger

estimación de error para los datos presente

| Cabecera para indicar valores de error |

| Número de puntos de datos

| Tipo de posición en x (centro)

| Sin IP presente

Tipo de estimación de error (0 = misma unidad como datos)

| Encabezado de tipo de error

0

| 0 para indicar que es la misma unidad que los datos | valor

1,500 1,000 1,000 0,1434 10,0145 1,000 1,000 9,9706 2,500 0,1882 3,500 1,000 1,000 0,1761 10,3357

| x-loc., 'a' valor, 'n' valor, aparente | resistividad, valor de error punto de datos

|

Tenga en cuenta que hay dos nuevas incorporaciones en el formato de datos. En primer lugar, justo antes de la sección del archivo con los puntos de datos de resistividad aparente medidos, hay 3 líneas adicionales. 'Estimado de error para los datos actuales' La primera línea es una cabecera para indicar que el archivo de datos contiene las estimaciones de error de datos. La segunda línea es una nueva cabecera para explicar el uso de la siguiente línea. La tercera línea contiene una bandera para indicar al programa de las unidades utilizadas en la estimación de error. En la actuali dad, utilizar un valor de '0' para indicar que el valor de error tiene las mismas unidades que el valor de datos. Por ejemplo, si las mediciones se dan como ohm.m, a continuación, los valores de error deben también estar en ohm.m Las mediciones en la mayoría de los instrumentos se dan como valores de resistencia (es decir ohmios), o tensión más los valores actuales que se pueden convertir a la resistencia valores. Para calcular el valor de resistividad aparente, el valor de resistencia se multiplica por un factor geométrico. El valor de error registrados por el instrumento se encuentra probablemente en la forma de un valor de resistencia (ohm). Por lo tanto la estimación del error registrado por el instrumento también debe multiplicarse por el mismo

128

factor geométrico antes de que se introduce en el formato de archivo DAT utilizada por este programa.

El uso de la bandera para la unidad utilizada por el valor de error hace posible que las futuras versiones del programa para permitir el uso de otras unidades para el valor de error, por ejemplo los datos en ohm.m mientras que el error en ohm. Otra posibilidad es el error d ado como un valor de porcentaje.

La ecuación de mínimos cuadrados en (2) se modifica con la siguiente ecuación

(

TdT

d

+  μ

)

=F g WWJ Tdd FT JWWJ -

 μ 

q

k

d

dónde

W eres una matriz de ponderación diagonal que incorpora el efecto de los errores de datos. Los

puntos de datos con un error menor se les da un mayor peso en el proceso de inversión.

En una situación ideal, la estimación del ruido se debe medir (por ejemplo, por mediciones recíprocos). Sin embargo, en la mayoría de los estudios de campo tales estimaciones de ruido no se miden por falta de tiempo. En tal situación, la mejor opción es probablemente usar la opción 'Robusta restricción de datos' (Anexo K) de manera que el efecto de los valores atípicos en los datos sobre el modelo de inversión se reduce al mínimo. Con est a opción, se le da un peso lineal a la diferencia entre los valores de resistividad aparente medidos y calculados. Esto hace que el resultado de inversión menos sensible a los valores atípicos en comparación con los mínimos cuadrados normales restricción donde se utiliza el cuadrado de la diferencia entre los valores de resistividad aparente medidos y calculados.

129

 Apéndice T: La incorporación de límites de capas conocidos En algunos casos, el perfil de profundidad de un límite entre dos capas con diferentes valores de resistividad se conoce a partir de otras fuentes; tales como datos de las perforaciones o sondeos de refracción sísmica. Con mucha frecuencia, se produce un cambio rápido e intenso a través del límite, por ejemplo, entre una capa de arcilla y arena. Si s e conoce la profundidad a la fr ontera, esta información puede ser incorporado en la inversión de manera que las resistividades de las capas pueden ser más exactamente determinado.

El CLIFLAYER4.DAT archivo es un archivo de datos de ejemplo en el que el límite entre una grava y la capa de arcilla subyacente se conoce a partir de una encuesta de refracción sísmica (Scott et al. 2000). Parte del a rchivo se muestra a continuación con una explicación de los cambios necesarios.

CLIFLAYER4.DAT

comentarios

Clifton

| Título

2 1 392 1 0 3 5 7

| Unidad de separación de los electrodos

| Wenner | Número de puntos de datos

| Tipo de posición en x (centro)

| Sin IP 2 2 2

58.63305 56.77204 54.79331

. . 48 50 0 0 0

| x-loc., 'a' de valor, la resistividad aparente | Segundo punto de datos | Tercer punto de datos | Otros puntos de datos

| 32 32

18.1012

|

17.94035

| El último punto de datos

| sin topografía | No hay regiones fijas

| No capa de agua

capas de Sharp presente

| Cabecera para indicar límites precisos presente

Número de límites

| La información de cabecera

2

| Número de límites, dos en este caso

De contorno 1 Parámetros

| Encabezado por el primer límite

tipo de límite Capa

| Cabecera de tipo de límite

grupo de peso

| Cabecera de peso, para su uso futuro

1.0 Número de puntos de frontera

| Introduzca 1.0 aquí, no se utilizan en la actualidad

dieciséis

| 16 puntos para este límite

X y Z coordenadas de pts 0.0 2.3

| Encabezamiento

18,0 2,3

| x-ubicación y la profundidad de segundo punto

. .

| x-ubicación y profundidad de otros puntos

88,0 3,2

| x-ubicación y profundidad de la segunda último punto

| tipo de capa

| Encabezamiento

| x-ubicación y profundidad de la primera punto

| a lo largo de la frontera

130

92,0 3,2

| x-ubicación y profundidad de último punto

Boundary 2 parámetros

| Encabezado por segundo límite

tipo de límite Capa grupo Peso

| Cabecera de tipo de límite

1.0 Número de puntos de frontera

| Introduzca 1.0 aquí, no se utilizan en la actualidad

3

| 3 puntos por este límite

X y Z coordenadas de puntos 0,0 6,7 45 6,7 92 6,7 0 0 0

| Encabezamiento

| tipo de capa | Cabecera de peso, para su uso futuro | Encabezamiento

| x-ubicación y profundidad de la primera punto | x-ubicación y la profundidad de segundo punto | x-ubicación y profundidad de último punto | terminar con algunos de cero

| |

En el ejemplo anterior, hay dos límites. Si hay más límites, la sección con 'frontera 3 parámetros' hasta la posición en x y la profundidad de la última a lo largo punto el segundo límite repetido. los es

archivo

MODEL_BOUNDARY_GENERAL.DAT da un ejemplo de un archivo de datos con el formato general matriz con un solo límite. La Figura 38 muestra la inversión de los datos CLIFLAYER4.DAT est ablecidos utilizando diferentes opciones de inversión. El primer modelo utilizando el método de inversión suave m uestra una transición gradational entre la arena y capa de arcilla resistividad inferior (Figura 38b) subyacente. La profundidad a la frontera no está bien determinado. El uso del método de inversión de bloques o sólida (Loke et al.

2003) da un límite mucho más agudo (Figura 38c). La adición de la frontera de la encuesta de refracción sísmica ayuda a afinar aún más el límite (Figura 38d).

La mayor limitación de esta opción en la actualidad es que las profundidades a las capas se fijan a lo largo de la inversión. Una opción para permitir a las profundidades para ser cambiados automáticamente probablemente se añadió en una futura actualización.

131

Figura 38. Modelos para la encuesta Clifton de los diferentes métodos de i nversión se muestran conjuntamente. Tenga en cuenta el límite entre la capa de arena de alta resistividad y la capa subyacente de baja resistividad de la arcilla está mucho mejor resuelto en los métodos de inversión límite de bloque y afilados.

132

 Apéndice T: Información especial topografía de encuestas de exploración minera

El software RES2DINV fue escrito originalmente para encuestas de i ngeniería y geofísicos ambientales. Por lo tanto la configuración por defecto en el programa siguen las prácticas en este campo. Una diferencia importante es que en encuestas medioambientales de ingeniería /, l os electrodos se disponen de tal manera que la distancia a lo largo de la superficie del terreno es constante mediante el uso de un cable multiconductor con una separación constante entre los takeouts. Esto ahorra tiempo y el costo de llevar a cabo un levantamiento topográfico detallado antes de que el estudio de resistividad. El programa RES2DINV asume que los X-localizaciones de los electrodos se dan en términos de distancias a lo largo de l a superficie del suelo.

Sin embargo para los estudios de exploración de minerales (con presupuestos más grandes) es común para llevar a cabo un levantamiento topográfico detallado antes de que el estudio de resistividad e IP, y los electrodos se colocan de tal manera que la distancia horizontal real (en lugar de la di stancia de la superficie) entre los electrodos es constante . En áreas con topografía empinada, es posible que la distancia vertical entre dos electrodos adyacentes a exceder la distancia horizontal (pero todavía menor que la distancia de la superficie). Para las encuestas donde se utiliza la verdadera distancia horizontal para el xlocation de los electrodos, se debe usar el formato de datos de la matriz general (véase el Apéndice L). Además, el indicador para el tipo de x-ubicación en el archivo de datos se debe establecer en 1 (véase el ejemplo MIXED.DAT en la página 91, y también la descripción de un caso tal en las páginas 93 y 94).

133

V Apéndice: Trazado de métodos para la matriz de gradiente.  Actualmente no hay ningún método estándar para representar los datos de resistividad aparen te a partir de encuestas con matriz de gradiente. Probablemente hay dos tipos de encuestas con la matriz de gradiente. La primera y tradicional tipo encuesta uso muy pocas posiciones C1C2 diferentes (por lo general sólo una con C1 y C2 situados en los dos extremos de la línea), pero muchas diferentes posiciones P1-P2. Esto se utiliza principalmente en las encuestas de exploración minera.

El segundo y moderno tipo, utilizado principalmente con sistemas de medición de resistividad de múltiples electrodos multicanal, utiliza muchas diferentes posiciones C1-C2 (una matriz de g radiente múltiple). Para dar cabida a las dos posibilidades, dos tipos de datos de las opciones gráficas están disponibles. Como la mayoría de encuestas es probable que se lleven a cabo utilizando el segundo tipo de encuesta con múltiples posiciones C1-C2, esto se describe en más detalle. Los artículos de Dahlin y Zhou (2004, 2005) tienen alguna información sobre este tipo de encuestas.

En primer lugar, una breve descripción del formato de datos para datos de encuestas de matriz de gradiente. El formato de datos de la matriz general (Anexo L) debe ser utilizado con un número submatriz de 15. A continuación es la parte inicial del archivo GRADIENT.DAT ejemplo, con comentarios acerca del formato.

archivo GRADIENT.DAT

comentarios

matriz de gradiente

| Título

1.0 11 15

| Unidad de separación de los electrodos

| tipo Array (11 para array general)

| Sub-array tipo 15

Tipo de medida (. = 0 aplicación de resistividad, 1 = resistencia) | cabecera 0

| 0 aquí para indicar valores de resistividad aparente

1792 1 0

| Número de puntos de datos | Indicador para el tipo de posición en x | No hay valores de IP, 1 si es IP presente

4 0.0,0.0 9.0,0.0 1.0,0.0 2.0,0.0 5,0588 | El formato para cada punto de datos es 4 1.0,0.0 10.0,0.0 2.0,0.0 3.0,0.0 5,0724 | Número de electrodos utilizados en . .

4 46.0,0.0 118.0,0.0 78.0,0.0 86.0,0.0 4,2571 | Datos de la última dos

4 48.0,0.0 120.0,0.0 80.0,0.0 88.0,0.0 4,2258 | puntos 0,0,0,0,0

| Terminar con unos cuantos ceros

Para seleccionar el tipo de trazado método a utilizar, antes de leer el archivo de datos, haga clic en 'Cambiar configuración' y luego '/ Selección de visualización de datos' y, a continuación, 'Tipo de parcela matriz de gradiente' como se muestra en la siguiente figura.

134

 Al seleccionar el 'Tipo de matriz de gradiente' sub-opción, se muestra el siguiente cuadro de diálogo.

El 'Use pseudodepth exacta' elección es la opción por defecto. Sin embargo, el programa cambiará automáticamente al método de pseudodepth aproximada si se encuentra el método pseudodepth exacta no es adecuado para el archivo de datos. Esto ocurre por encue stas con el tipo tradicional de matrices de gradiente, donde sólo se utilizan algunas posiciones C1-C2.

El archivo CROMER02.DAT tiene los datos de una encuesta utilizando una matriz de gradiente múltiple se llevó a cabo por la Universidad de Aarhus de Dinamarca para el camino y la Asociación de la carretera. Una trama de este conjunto de datos en forma de perfiles se muestra en la siguiente figura.

135

Figura 39. Parcela de datos CROMER02.DAT en forma de perfiles. Tenga en cuenta el conjunto de datos cada pseudodepth se divide en dos perfiles que se dibujan en diferentes colores. Esto es porque hay dos posibles configuraciones con el mismo pseudodepth para configuraciones no simétricas de la matriz de gradiente.

Figura 40. dos configuraciones diferentes de la matriz de gradiente con el mismo pseudodepth pero con los electrodos de potencial a diferentes extremos de la matriz. configuración de la matriz de gradiente con (a) la paridad derecha y (b) la paridad izquierdo.

En el ejemplo anterior, la distancia entre los electrodos C1 y C2 son los mismos. El desplazamiento del dipolo P1-P2 desde el electrodo C1 en (a) es el mismo que el desplazamiento del dipolo desde el electrodo C2 en (b). Sin embargo, las dos configuraciones son sensibles a diferentes partes de la sub-superficie (Dahlin y

136

Zhou 2005). Para proporcionar una parcela sin problemas para el perfil, los datos de la 'derecha' y configuraciones de 'izquierda' se representan en diferentes perfiles.

Si el programa encuentra que no es práctico para trazar los perfiles utilizando el pseudodepth exacta, cambiará automáticamente al método de pseudodepth aproximada. Para cada separación C1-C2, los datos se divide en más de 2 perfiles qu e se trazan usando diferentes colores.

Figura 41. Parcela de datos CROMER02.DAT en forma de perfiles utilizando el método pseudodepth aproximada.

137

El modelo de inversión del conjunto de datos CROMER02.DAT que también muestra la gráfica de los datos en pseudosecciones de resistividad aparente se muestra a continuación.

Figura 42. El resultado inversión de los datos CROMER02.DAT establecidos con las secciones de re sistividad y de modelo aparentes.

138

W Apéndice: La matriz polo-dipolo de compensación

Esta matriz se diseñó originalmente para llevar a cabo 3-D encuestas IP rápidamente (White et al 2003), de este modo los datos de estos estudios son normalmente procesados con el RES3DINV programa de inversión de 3-D. El apoyo se añadió a esta matriz en este programa para que los usuarios tienen la opción de invertir cada línea individual. La disposición de los electrodos se muestra en la figura siguiente.

Figura 43. Disposición de electrodos para la matriz polo -dipolo offset. La disposición es similar a la utilizada por la matriz normal de polo-dipolo, excepto que el electrodo de corriente no está en la misma línea que los electrodos potenciales, pero compensado a una distancia “b”

desde la línea de encuesta. Esta disposición se utiliza para reducir el acoplamiento EM entre el electrodo de corriente y el dipolo potencial. El archivo OFF_PDP.DAT da un ejemplo del formato de datos para esta matriz. A continuación se muestra la parte inicial del archivo GRADIENT.DAT ejemplo, con comentarios acerca del formato.

archivo OFF_PDP.DAT

comentarios

-----------------------------------

| ------------------------------------------------- ------

Bloques con IP

| Título

2.00 10 1.00

| Unidad de separación de los electrodos

530

| Número de puntos de datos

1

| Para indicar el punto medio de la matriz dada como x-ubicación

1

| 1 para indicar IP presente

cargabilidad

| Tipo de datos IP

mseg

| unidad de IP

0.10,1.00

| Retardo, tiempo de integración

2,0 2,0 1,0 100,7976 -0,9118 4,0 2,0 1,0 102,0548 -1,8678 6,0 2,0 1,0 104,9487 -4,3391 8,0 2,0 1,0 111,0159 -9,0335

| Para cada punto de datos, tenemos

| tipo de matriz (10 para compensar polo-dipolo) | Desplazamiento del electrodo de corriente desde la progresiva

| x-ubicación, “a” espaciado dipolo, factor “n”, | resistividad aparente, IP aparente |

139

... ..

| Otros puntos de datos

... ..

|

53.0 5.5469

2.0

-10.0

109.3857

0,0,0,0,0

| El último punto de datos | Unos pocos ceros para terminar el archivo

El formato de datos es esencialmente la misma que la matriz normal de polo-dipolo, excepto el número de matriz es de 10 y la distancia de desplazamiento se da inmediatamente después del número de array.

El OFF_PDP_GEN.DAT archivo es un archivo de datos de ejemplo con datos dados en el formato general

la parte inicial de

matriz. A continuación es

el ejemplo

GRADIENT.DAT archivo con el co mments acerca del formato. archivo OFF_PDP_GEN.DAT

comentarios

-----------------------------------

| ------------------------------------------------- ------

Bloques con IP

| Título

2.00 11 10

| Unidad de separación de los electrodos

1.00

| Desplazamiento del electrodo de corriente desde la progresiva

Tipo de medición (0 = aplicación. Resistividad,

| Encabezamiento

| número de matriz 11 para array general | Sub- un tipo RRay (10 para compensar polo-dipolo)

1 = resistencia) 0

| 0 aquí para indicar valores de resistividad aparente

530

| Número de puntos de datos

1

| Para indicar el punto medio de la matriz dada como x-ubicación

1

| 1 para indicar IP presente

cargabilidad

| Tipo de datos IP

mseg

| unidad de IP

0.10,1.00

| Retardo, tiempo de integración

2.0,0.0 4.0,0.0 3 0.0,0.0 100.7976 -0,9118 3 2.0,0.0 4.0,0.0 6.0,0.0 102.0548 -1,8678 6.0,0.0 8.0,0.0 3 4.0,0.0 104.9487 -4,3391 3 6.0,0.0 8.0,0.0 10.0,0.0 111.0159 -9,0335

| Para cada punto de datos, tenemos

... ..

| Otros puntos de datos

... .. 3 64.0,0.0 44.0,0.0 42.0,0.0 109.3857 5.5469 0,0,0,0,0

|

| Número de electrodos usados en la medición,

| x y z-ubicación de electrodo C1, | x y z-ubicación de electrodo P1 | x y z-ubicación de electrodo P2 ,, | resistividad aparente, IP aparente

| El último punto de datos | Unos pocos ceros para terminar el archivo

Tenga en cuenta que el offset “B” debe ser el mismo para todas las mediciones en una sola línea encuesta.

140

 Apéndice X: coordenadas globales o GPS Bidimensional estudios de resistividad se llevan a cabo con los electrodos dispuestos a lo largo de una línea recta. Las coordenadas de las posiciones de los electrodos, y por lo tanto los valores de resistividad modelo del subsuelo, se dan con respecto a los extremos de la línea en forma de coordenadas locales. A fin de incorporar los resultados en un mapa, se requiere que las coordenadas verdaderos (en términos de latitudes y longitudes, o coordenadas X e Y). Las posiciones de los puntos a lo largo de la línea de ahora se pueden obtener fácilmente utilizando el Sistema de Posicionamiento Global (GPS). Con el fin de obtener el los verdaderos valores de coordenadas de los modelos, las coordenadas GPS primero se deben introducir en el archivo de datos leídos por el programa RES2DINV. Los resultados más adelante pueden ser extraídos del archivo XYZ generado a partir del archivo de la inversión de datos.

La información relativa a las coordenadas GPS se da en la sección inmediatamente después de la información de la topografía en el archivo de datos. Un ejemplo de un archivo de datos (sin topografía) utilizando el formato basado en el índice para los datos de resistividad se da en el BETA_GLOBAL.DAT archivo. La sección del a rchivo que contiene la información del GPS se muestra a continuación.

archivo BETA_GLOBAL.DAT

-----------------------------------

comentarios

-------------------------------------------------------

8.750

4.500 10.132

Último

7,500

5,000

9,601

Tres

8,000

5,000

9,601

re puntos ATA

0

bandera Topografía, 0 = no topografía

Global Coordenadas presente

Cabecera para indicar la información del GPS actual

Número de puntos de coordenadas

Encabezamiento

3

Número de puntos de coordenadas

Longitud Latitud local 0,0 20,0 10,0

Encabezamiento Primera coordenada dada como punto

la posición local a lo largo de la línea, GPS Longitud, Latitud GPS

8,0 26,93 14,0 15.5 33.42 17.75

Segundo punto de coordenadas

0

Ceros para indicar otras opciones,

0

tales como regiones fijas etc.

0

no presente

Tercer punto de coordenadas

Un ejemplo de un archivo de datos con la topografía utilizando el formato basado en el índice se da en el RATCHRO_GLOBAL.DAT archivo. La sección del archivo que contiene l a información del GPS se administra después de la sección de topografía se muestra a continuación.

141

RATHCRO_GLOBAL.DAT

comentarios

archivo

-----------------------------------

-------------------------------------------------------

92, -3.207

Último

94, -3.27

topografía tres

96, -3.34

puntos de datos

1

Indica primer punto topografía coincide con el primer electrodo

Global Coordenadas presente

Cabecera para indicar la información del GPS actual

Número de puntos de coordenadas

Encabezamiento

3

Número de puntos de coordenadas

Longitud Latitud local -36,0 80.100,0 200.300,0

Encabezamiento Primera coordenada dada como punto

la posición local a lo largo de la línea, GPS Longitud, Latitud GPS

20,0 80.148,1 200.327,8

Segundo punto de coordenadas

96,0 80.213,3 200.365,4

Tercer punto de coordenadas

0

Ceros para indicar otras opciones,

0

tales como regiones fijas etc.

0

no presente

El archivo IP_MAGUSI_MF_GLOBAL.DAT es un ejemplo de un archivo de datos con los valores de IP, así la información del GPS. El MIXED_GLOBAL.DAT archivo y RATHMIX_GLOBAL.DAT son ejemplos de archivos de datos utilizando el formato de matriz general que contiene las coordenadas GPS. Ambos conjuntos de datos tienen la información de la topografía de cada curva de punto de datos. La sección de la RATHMIX_GLOBAL.DAT con la información del GPS se muestra a continuación.

archivo RATHMIX_GLOBAL.DAT

-----------------------------------

comentarios

--------------------------------------------------

4 60,0 -0,093 96,0 -3,34 72,0 -2,274 84.0 -2.914 1680.0

Última aparente punto de datos de resistividad

0

Ponga cero para la bandera de la topografía

Global Coordenadas presente

Cabecera para indicar la información del GPS actual

Número de puntos de coordenadas

Encabezamiento

3

Número de puntos de coordenadas

Longitud Latitud local -36,0 80.100,0 200.300,0

Encabezamiento

En primer punto de coordenadas local da como posición a lo largo de la línea, GPS Longitud, Latitud GPS

20,0 80.148,1 200.327,8 96,0 80.213,3 200.365,4

Segundo punto de coordenadas

0

Ceros para indicar otras opciones,

0

tales como regiones fijas etc.

0

no presente

Tercer punto de coordenadas

142

El LAKELELIA_GLOBAL.DAT archivo es un ejemplo con datos en el formato general para el conjunto de una encuesta bajo el agua con las coordenadas GPS. También tiene la información de la topografía en una lista separada de las líneas de datos de resistividad aparente. Parte de este archivo se muestra a continuación.

archivo RATHMIX_GLOBAL.DAT

-----------------------------------

comentarios

--------------------------------------------------

4 1056.62,0.0 1083.62,0.0 1065.62,0.0 1074.62,0.0 95.383

El último punto de datos

Topografía en lista separada

Cabecera para obtener información topografía

2

bandera de la topografía

38

Número de puntos de topografía

26.12 999.09

Primer punto de la topografía

... ..

otro topografía

... .. 1.127,76 999,09

puntos de datos

1

Indica primer punto topografía coincide con el primer electrodo

Número de puntos de coordenadas

Encabezamiento

4

Número de puntos de coordenadas

Longitud Latitud local

Encabezamiento

26,0 100.000,0 200.000,0

En primer punto de coordenadas local da como posición a lo largo

El último punto topografía

de la línea, GPS Longitud, Latitud GPS

126,0 100.086,6 200.086,6 500,0 100.410,5 200.237,0 1.092,0 100.923,2 200.533,0

Segundo punto de coordenadas

0

Ceros para indicar otras opciones,

0

tales como regiones fijas etc., no presente

1

Indica capa de agua presente

58.00, -2000.00,2000.00,1000,1

La información de la capa de agua

0,0,0

Terminar con unos cuantos ceros

Tercer punto de coordenadas Último punto de coordenadas

Para hacer uso de las coordenadas de GPS, es necesario ejecutar primero una inversión del conjunto de datos. A continuación lee el a rchivo de la inversión en la ventana 'Pantalla', y luego guardar los valores de modelo en un archivo XYZ a través del 'Archivo - Guardar los d atos XYZ en formato' opción. Los valores de modelo con las coordenadas GPS se mostrarán en la última sección del archivo XYZ.

143

Nota especial OhmMapper Los datos recogidos por el sistema OhmMapper se pueden convertir en el formato utilizado por el programa RES2DINV a través del 'DataMap para OhmMapper' programa proporcionado por Geometrics, Inc. Este es un programa de Windows 95/98/2000 / NT basado. Después de iniciar el programa, por favor haga clic en la opción de exportación para llevar a cabo la conversión de datos. Para saber más detalles, consulte el manual proporcionado por OhmMapper Geometrics.

La figura a continuación muestra un ejemplo de un conjunto de datos recogidos por el sistema OhmMapper  junto con el modelo de inversión producida por el programa RES2DINV. La encuesta se llevó a cabo en un área con granito degradado. El conjunto de datos s e da en el archivo OHMMAPPER.DAT. Las mediciones se realizaron utilizando la matriz dipolo-dipolo con una longitud dipolo "a" espaciamient o de 10 metros y con n valores de 0,5 a 3,5 (Apéndice A). El conjunto de datos tiene 732 puntos de datos mientras que el modelo de inversión consiste en 1240 células. La inversión de este conjunto de datos se llevó unos 15 minutos en un equipo de 300 MHz Pentium II.

Figura 44: (a) pseudosección resistividad aparente medido con un sistema móvil resistividad topografía OhmMapper. sección (b) Modelo obtenido por programa RES2DINV.

144

RENUNCIA Este software se proporciona "tal cual" si n ninguna garantía expresa o implícita, incluyendo su idoneidad para un propósito particular. Ni el autor ni el subagente asumirán la responsabilidad por cualquier daño o pérdida causada por el uso de este programa. En el caso de errores en el programa, se hará todo lo posible para corregir los errores y se le proporcionará un programa de sustitución por un cargo nominal para cubrir los gastos de envío. dispositivos de seguridad defectuosos o discos del programa serán reemplazados de forma gratuita dentro de los primeros 6 meses desde la fecha de compra para usuarios registrados. Los usuarios registrados tienen derecho a actualizaciones gratuitas y soporte dentro del primer año.

 APOYO Por favor, póngase en contacto con su distribuidor si tiene alguna pregunta acerca de este programa. Una fuente frecuente de consultas es cuando los resultados obtenidos no son lo que se espera, o el programa se niega a aceptar el archivo de datos. En tal caso, sólo se proporcionó apoyo si es una copia del archivo de datos se envían  juntos como un archivo adjunto con su correo electrónico, preferentemente comprimido en un archivo ZIP. Las causas más comunes de problemas son: (a) los datos en el archivo DAT no se da en el formato correcto RES2DINV (b) el ruido excesivo en el conjunto de datos (sobre todo las encuestas dipolo-dipolo y / o conjuntos de datos IP) (c) mal funcionamiento del instrumento o de estudio de campo e rrores. Sólo es posible rastrear el origen del problema si se proporciona una copia del archivo de datos.

Tenga en cuenta que el apoyo sólo está disponible para usuarios registrados de la versión completa del programa. No es posible proporcionar apoyo técnico a los usuarios de la versión demo de este software. Para obtener soporte técnico, puede ponerse en contacto con el proveedor de software por correo electrónico a las siguientes direcciones. Por favor adjunte una copia de su archivo de datos a su correo electrónico.

[email protected] [email protected]

La siguiente información será requerida antes de la prestación de asistencia técnica.

1). Su número de ID del dongle. Este número de identificación se muestra cerca de la parte inferior de la caja de la información inicial cuando el programa se pone en marcha RES2DINV si el dispositivo de seguridad está conectado al ordenador.

2). El nombre de su empresa o institución. 3). fecha aproximada de su empresa / institución adquirió el software.

145

Qué hay de nuevo

Versión 3.00 - La mayor actualización a Windows 3.1 y 95. El número máximo de electrodos aumentó a 500. El soporte multi-tarea para Windows 95. Soporte para cualquier tarjeta gráfica compatible y la impresora de Windows.

Versión 3.01 - Mejoras en diferencias finitas subrutina para calcular la matriz Jacobiana que ahora es dos veces más rápido que las versiones anteriores. Versión 3.10 - Soporte para estudios submarinos. La incorporación directa de la topografía en el modelo de inversión utilizando una cuadrícula de elementos finitos distorsionada.

Versión 3.11 - El tiempo tomado por el de elementos finitos subrutina para calcular la matriz jacobiana se reduce. Se añade una opción de utilizar una malla más fina para la finitedifference o de elementos finitos método de modelado adelante.

Versión 3.12 - no enteros valores para el “ norte ”Factor para el dipolo -dipolo, poledipole y Wenner-Schlumberger matrices soportado (consulte el A péndice UN). Los límites máximos para el númer o de niveles de datos y las capas del mo delo aumentó a 52 y 17 respectivamente. Versión 3.13 - La proporción de la mayor distancia entre los electrodos a distancia entre electrodo s unidad para la matriz Wenner y el polo polos se aumentó de 32 a 64. El número máximo de niveles de datos es ahora 64. El número máximo de las capas del modelo también se aumentó a 18 .

Versión 3.20 - Compatibilidad con datos de IP añaden. Versión 3.21 - Se aumentó el número máximo de electrodos a 650. Las mejoras se han hecho en el disco de memoria de intercambio de subrutinas de manera que se incrementa de una determinada cantidad de memoria RAM el número de puntos de referencia el programa puede manejar. Para los sistemas informáticos con más de una unidad de disco duro, el programa seleccionará automáticamente la unidad con más espacio libre en disco para almacenar los archivos de intercambio temporal en disco.

Versión 3.22 - Ligeras mejoras en el uso de l a memoria en la sección de inversión de IP. Se aumentó el tamaño del conjunto de datos IP que puede ser manejado para una cantidad dada de la memoria. Se añadió una opción para optimizar el factor de amortiguamiento de forma automática durante el proceso de inversión. Soporte para la matriz de polo- dipolo “inverso” añadido (véase el Apéndice A). El apoyo a un modo de línea de comandos por lotes incluido (véase el Apéndice I). También se añadió una opción para trazar el modelo de l a sección en forma de bloques rectangulares.

146

Versión 3.30 - Soporte para encuestas de c orte del pozo añaden. Una opción para permitir que el número de parámetros del modelo exceda También se añadió el número de puntos de referencia.

Versión 3.31 - Una opción para un modelo definido utilizado añadió. En este modelo, el usuario especifica el espesor de la primera capa y el factor para aumentar el espesor de cada capa más profunda posterior. El programa también se actualizará automáticamente el directorio utilizado para los archivos de datos de entrada y de salida de los archivos de inversión incluidos en el archivo RES2DINV.INI.

Se añadió una opción para ampliar la subdivisión del subsuelo en bloques a los bordes de la línea de la encuesta - la versión 3.32. También se añadió Opción para mostrar la incertidumbre en los valores de resistividad modelo. Se añadió una opción para incorporar los efectos de los electrodos a distancia utilizados en estudios de resistividad superficial con los arrays de polos polos y el polo-dipolo - Versión 3.33. Versión 3.34 - Unas correcciones de errores. El usuario ahora puede escalar las profundidades de las capas del modelo de manera que la profundidad a la última capa puede ser mucho mayor que el límite predeterminado permitido por el programa. Cuando el programa guarda los valores de modelo en el formato XYZ, que ahora también guardar las coordenadas de las esquinas de los bloques en el modelo. Una opción para utilizar una malla muy fina en la dirección vertical para contrastes de resistividad de más de 250: 1 se añadió, así como una opción para 6 nodos en la dirección horizontal entre los electrodos adyacentes.

Versión 3.35 - El usuario puede ahora establecer el número máximo de electrodos, desde 150 a 1500, a través del programa JACOBWIN.EXE. Dibujo de las secciones de contorno de c olor es ahora mucho más rápido en la mayoría de los ordenadores. Versión 3.36 - Se añadió una opción de mínimos cuadrados inversión sólida. Para los datos muy ruidosos con '' valores atípicos, el modelo resultante será menos sensible a tales puntos de referencia cuando se selecciona el método de inversión de datos robusta. Para las áreas donde la geología del subsuelo tiene interfaces afilados, el robusto modelo de método de inversión dará mejores resultados.

Versión 3.40 - Soporte para las matrices no convencionales para estudios de resistividad. Versión 3.41 - Dos nuevos métodos para el modelado topográfico utilizando una cuadrícula distorsionada amortiguado y la transformación inversa de Schwartz-Christoffel añadió.

Versión 3.42 - Soporte para electrodos remotos para encuestas IP.

Versión 3.43 - Soporte para la matriz Wenner Gamma, y estudios submarinos con matrices no convencionales. número máximo de electrodos aumentó a 2000. El programa también ha sido optimizado para los conjuntos de datos donde la unidad

147

separación de los electrodos se ha reducido a la mitad del valor real con el fin de obtener un modelo en el que la anchura de los bloques es la mitad del tamaño usual. Esto ayuda en casos en los que hay grandes variaciones de resistividad lateral cerca de la superficie.

Versión 3.44 - Soporte para las encuestas IP con matrices no convencionales. Esta característica es útil en algunos casos donde los datos de una serie de superposición colineal 1D líneas de levantamiento de resonancia se pueden combinar en un solo dato 2-D SET para obtener un modelo 2-D.

Versión 3.45 - Opción de guardar los resultados en formato SURFER añaden. También soporte para encuestas con algunos electrodos bajo el agua y algunos electrodos por encima del nivel de la superficie del agua añadida.

Versión 3.46 - adiciones menores para los datos en formato general array. Opción de utilizar un modelo con la mitad de la distancia entre electrodos unidad para los datos en formato general array se añade (ver la RATCMIX2.DAT archivo para un ejemplo). Ahora es posible llevar a cabo la inversión de los datos en el formato general matriz con la versión de demostración, pero los resultados sólo se mostrará temporalmente en la pantalla durante la inversión. El método de elementos finitos rejilla amortiguada distorsionada se establece como el método por defecto para el modelado topográfico.

Se añadió una opción para llevar a cabo la i nversión de la secuencia de datos IP - Versión 3.47. Ligeros cambios en la estructura del menú. Versión 3.48 - Soporte para un máximo de 4 pozos de sondeo en la opción cruzada pozo.

Versión 3.49 - Incompleto opción de inversión de Gauss-Newton añadió. número máximo de electrodos aumentó a 4.000, y el número máximo de puntos de datos a 20000. opción de combinar una serie de ficheros de datos 2-D en formato RES2DINV en un único archivo de datos 3-D en el formato RES3DINV añadió. Versión 3.50 - detección de bordes roca de fondo y de lapso de tiempo opciones añadidas. Versión 3.51 - Opción para incluir estimaciones de ruido de datos incluidos.

Versión 3.52 - Opción de encuestas con electrodos flotantes añadió. El apoyo a la topografía en las encuestas transversal del pozo incluido. Versión 3.53 - opción de inversión escasa para las l íneas de campaña muy larga (2.000 a 16.000 posiciones de los electrodos) añadió. Este método invierte el conjunto completo de datos y el modelo en una sola vez para producir un modelo continuo y si n problemas. Nuevo formato de archivo de comandos por lotes opción de modo que hace uso de los archivos que contienen los parámetros de inversión producidos por el programa RES2DINV.