ISSN-e: 3006-9467; ISSN: 0016-7975 / 1011-9565
Petrofísica/Petrophysics/Petrofísica
Alfonso Quaglia
IngºGeóº, PhD., Inter-Rock USA. Correo-e: quagliaa@inter-rock-ca.com
Juan Carlos Porras
IngºGeóº, MSc., Inter Rock Panamerican. Correo-e: porrasjc@inter-rock-ca.com
Paola Almedo
IngºGeóº, Inter-Rock, C. A. Correo-e: almedopaola@gmail.com
Rafael Panesso
IngºGeóº, MSc., Inter Rock Panamerican Colombia. Correo-e: panessor@inter-rock-ca.com
Marlene Villalba
Ing°Petró°, PhD, Inter-Rock, C. A. Correo-e: Villalbam@inter-rock-ca.com
Recibido: 14-10-24; Aprobado: 18-11-24
Abstract
This work focuses on developing a workflow for defining net pay thicknesses in conventional reservoirs and comparing it with other relevant methods. The background research provided a fundamental understanding of potential variations in net pay thickness estimations resulting from different methods applied across the studied areas. This was achieved by developing petrophysical models for each area and their respective key wells designated as Well-1, Well-2, and Well-3, which included conventional well logs and core analysis data. Petrophysical properties derived from the developed petrophysical models include: volume of clay (Vclay), porosity (Phi), water saturation (Sw), and permeability (K), which were calibrated with conventional core analysis. The formations selected for this study were designated as X, Y, and Z, to which different cutoff parameter methods, recognized as the most relevant in net pay determination, were applied. The method yielding the highest net pay thicknesses was Shepherd's method, described in his 2009 book Oil Field Production Geology (AAPG Memoir 91, p. 189-193). However, this does not necessarily make it the most recommendable or reliable method. A comparative study with other existing methods is required, considering the specific characteristics of each reservoir under study. In Shepherd's method, the hydrocarbon pore volume thickness equals the gross thickness multiplied by the net-to-gross ratio, porosity, and oil saturation. While Shepherd's work doesn't explicitly define cutoff parameters, focusing primarily on net pay determination, these parameters were determined in this study using PhieSo crossplots. PhieSo represents the Oil Pore Volume, defined as the percentage of porosity occupied by oil and calculated as PHIE*(1-BVW) (Inter-Rock, C. A., 2011).
El objetivo de este trabajo consistió en el desarrollo de un flujo de trabajo en yacimientos convencionales para la definición de espesores y su comparación con algunos métodos relevantes. Los antecedentes de esta investigación representaron los pilares fundamentales para entender las eventuales variaciones de los resultados de los métodos empleados en las áreas de estudio para la determinación de los espesores útiles. Esto fue posible luego de desarrollar los modelos petrofísicos para cada área y sus respectivos pozos clave designados como Pozo-1, Pozo-2 y Pozo-3, los cuales contaban con registros convencionales y análisis de núcleos. Entre las propiedades petrofísicas derivadas de los modelos petrofísicos desarrollados están: volumen de arcilla (Vclay), porosidad (Phi), saturación de agua (Sw) y permeabilidad (K), las cuales se calibraron con análisis convencionales de núcleos. Las formaciones que sirvieron de objeto para este estudio se designaron como X, Y, y Z a las cuales se les aplicaron distintos métodos de parámetros de corte conocidos como los más relevantes en materia de determinación de espesores útiles. Se observó que el método que logro resultados con valores más altos de espesores útiles fue el método de Shepherd descrito en su libro publicado en 2009. Oil field production geology. AAPG Memoir 91, p. 189-193. Esto no lo hace, en absoluto, el método más recomendable o de menor certidumbre. Para ello se requiere de un estudio comparativo con el resto de los métodos existentes y tener presente las características propias de cada yacimiento estudiado. En este método el espesor del volumen poroso hidrocarburífero es igual al espesor bruto x la relación del espesor neto con respecto al bruto x la porosidad x saturación de petróleo. En su trabajo, Shepherd no menciona cómo determinar los parámetros de corte, mayormente se enfoca en la determinación de los espesores útiles. Sin embargo, para este método los parámetros de corte se obtuvieron mediante los gráficos de la relación PhieSo, la cual corresponde al Volumen Poroso de petróleo y está definido como el porcentaje de la porosidad que se encuentra ocupado por petróleo y es igual a PHIE*(1-BVW) (Inter-Rock, C. A., 2011).
Este trabalho teve como objetivo o desenvolvimento de um fluxo de trabalho em reservatórios convencionais para a definição de espessuras líquidas e sua comparação com alguns métodos relevantes. O embasamento desta pesquisa representou o pilar fundamental para entender as eventuais variações nos resultados das espessuras líquidas calculadas pelos métodos empregados nas áreas de estudo. Isso foi possível após o desenvolvimento dos modelos petrofísicos para cada área e seus respectivos poços chave designados como Poço-1, Poço-2 e Poço-3, os quais continham perfis convencionais e análises de núcleos. Entre as propriedades petrofísicas derivadas dos modelos petrofísicos desenvolvidos estão: volume de argila (Varg), porosidade (Phi), saturação de água (Sw) e permeabilidade (K), as quais foram calibradas com análises convencionais de núcleos. As formações que serviram de objeto para este estudo foram designadas como X, Y e Z, às quais foram aplicados distintos métodos de parâmetros de corte, reconhecidos como os mais relevantes em termos de determinação de espessuras líquidas. Observou-se que o método que obteve resultados com maiores valores de espessuras líquidas foi o método de Shepherd descrito em seu livro publicado em 2009, Oil Field Production Geology (AAPG Memoir 91, p. 189-193). Isso não o torna, de forma alguma, o método mais recomendável ou de menor incerteza. Para isso, se faz necessário um estudo comparativo com os demais métodos existentes e levar em consideração as características específicas de cada reservatório estudado. Neste método, a espessura do volume poroso de hidrocarboneto é igual à espessura bruta multiplicada pela relação entre a espessura líquida e a bruta, pela porosidade e pela saturação de óleo. No trabalho de Shepherd, não há menção de como determinar os parâmetros de corte, o foco principal é a determinação das espessuras líquidas. Entretanto, para este método os parâmetros de corte foram obtidos através dos gráficos da relação PhieSo, que corresponde ao Volume Poroso de Óleo e é definido como a porcentagem da porosidade que se encontra ocupada por óleo, sendo igual a PHIE*(1-BVW) (Inter-Rock, C. A., 2011).
Palabras clave/Keywords/Palabras-chave:
Cutoff parameters, espessura líquida, espesor neto, modelos petrofísicos, net pay, parâmetros de corte, parámetros de corte, porosidade, porosidade, porosity, saturação de água, saturación de agua, water saturation.
Citar así/Cite like this/Citação assim: Quaglia et al. (2024) o (Quaglia et al, 2024).
Referenciar así/Reference like this/Referência como esta:
Quaglia, A., Porras, J. C., Almedo, P., Panesso, R., Villalba, M. (2024). Determinación y análisis comparativos de espesores útiles en caracterización de yacimientos convencionales empleando métodos existentes relevantes. Casos de estudio en norte y sur américa. Geominas 52(93). 61-84.
Introducción
Las metodologías consideradas como universales o estándares para la determinación de espesores útiles en ocasiones presentan incertidumbre, una de las razones es la cantidad de factores que se deben considerar al elegir un valor de corte. De hecho, el problema con frecuencia sobrepasa cualquier enfoque sistemático y analítico y, con frecuencia, la solución debe gestionarse según las particularidades de cada caso. Esto a su vez implica que siempre habrá un grado considerable de subjetividad e incertidumbre en este proceso.
Existen diversos métodos para determinar espesores útiles, sin embargo, ninguno de ellos podría considerarse como definitivo, sobre todo tomando en cuenta los diferentes propósitos para los cuales estos se han determinado. Podría argumentarse que el uso más frecuente de los espesores útiles es el cálculo de hidrocarburo original en sitio (POES o GOES).
Usualmente, los espesores útiles para esta aplicación están relacionados con la determinación de un factor de recobro en el yacimiento. Otro uso de los espesores útiles tiene que ver con la definición del volumen total de hidrocarburos en el yacimiento con la finalidad de determinar las reservas y el potencial comercial del mismo. Una tercera razón es para evaluar la factibilidad de implementar recobros secundarios en los yacimientos.
Para determinar espesores útiles se requiere definir parámetros de corte, los cuales deben fijarse según las características de cada yacimiento y del propósito para el cual se establecen. Es decir, tienen un carácter variable y dinámico, debiendo ser ajustados al estado actual del yacimiento estudiado, y a las condiciones económicas y tecnológicas al momento de ser definidos. La necesidad de crear flujos de trabajo para la aplicación, análisis y resolución de ecuaciones y fórmulas matemáticas es cada vez más común en empresas que requieren de tecnología para ejecutar diversos estudios. La finalidad primordial en la aplicación de tecnología y los mencionados flujos de trabajo se fundamenta en recortar los tiempos de ejecución de los proyectos y en optimizar los costos sin sacrificar la calidad de los resultados.
La industria petrolera tiene como propósito determinar los espesores útiles y las propiedades petrofísicas de los yacimientos para la estimación cada vez más precisa de los volúmenes de hidrocarburo. Es por ello que es de suma importancia para los geólogos e ingenieros la optimización de flujos de trabajo que permitan definir la calidad y el comportamiento de producción de los yacimientos petroleros, disminuyendo en lo posible los niveles de incertidumbre.
Finalmente, este trabajo plantea el desarrollo de un flujo de trabajo para determinar espesores útiles mediante la optimización del proceso de selección de parámetros de corte y la comparación entre las metodologías existentes más relevantes.
Generalidades del Área de Estudio
Los datos empleados en este trabajo están sujetos a pautas de confidencialidad, por lo que no se proporcionará información detallada sobre la ubicación geográfica exacta o los nombres reales de los pozos investigados.
Ubicación geográfica: Los yacimientos convencionales a estudiar se ubican geográficamente en Norte y Sur América. Para el siguiente trabajo fueron utilizados tres (3) pozos. En la tabla 1 se identifican los pozos cargados a la base de datos, su ubicación relativa, las edades geológicas, topes y bases de las formaciones de las zonas de estudio.
Tabla 1. Ubicación geográfica de los pozos, edad geológica y topes de las formaciones.
Descripción geológica: con respecto a la estratigrafía del “POZO-1”, este atraviesa a la formación “X” de edad Cretácico (Aptiano medio-Albiano inferior). Consta de areniscas cuarzosas, grano-decrecientes hacia el tope, generalmente gruesas, blancas, poco consolidadas, con estratificaciones cruzadas y rizaduras ocasionales. Posee intercalaciones de lutitas limolíticas y lutitas bituminosas de color oscuro. La mayor parte de la Formación “X” se depositó en un ambiente continental fluvial anastomosado. Sin embargo, a fines de la depositación invadieron condiciones marinas, con aumento de lutitas, presencia de glauconita y unos pocos fósiles marinos (Tschopp, 1953) (Figura 1).
Figura 1. Columna estratigráfica POZO-1, Formación “X”, Cretácico inferior.
Con respecto a la Estratigrafía del “POZO-2”, este atraviesa la formación designada con la letra “Y”, la cual representa una columna litológica de origen Jurásico Medio, comprende un espesor hasta de 1.000 metros y está compuesta por areniscas finas a gruesas hasta conglomeráticas, apretadas (tight), intercaladas con arcilitas, dividida en tres miembros con permeabilidades que van de 1 a 100 mD (miliDarcy). Esta formación fue depositada en un ambiente fluvio-deltaico, dominado con acción subordinada de mareas y olas, con asociación de facies correspondientes a planicies deltaicas distales, canales distributarios y frentes deltaicos (Wagner, Alencastre, Terrasanta y Masiero, 2016) (Figura 2).
Figura 2. Columna estratigráfica POZO-2, Formación “Y”, Jurásico Medio.
Figura 3. Columna estratigráfica POZO-3, Formación “Z”, Pensilvánico superior.
Por último, la estratigrafía del “POZO-3”, el cual atraviesa una columna litológica de areniscas ubicada en Norte América, está conformada por la formación designada arbitrariamente con la letra “Z” y cuya secuencia fue depositada en el periodo Pensilvánico superior. La zona productiva está compuesta por areniscas de grano fino a muy fino, consolidadas, con un rango de porosidad promedio que está en el orden de 8 a 10%, con pequeñas cantidades de feldespato, intercalado con lutitas grises. (Branch, G.A., y Drennan, K. M, 1991) (Figura 3).
Antecedentes
Inter-Rock, C. A. (2018), en su trabajo “CARACTERIZACIÓN PETROFÍSICA DEL YACIMIENTO LINDERO ATRAVESADO, GRUPO CUYO, CUENCA NEUQUINA” utiliza la combinación de RMN y Cluster Analysis, junto con una modificación de la metodología de Scribner, para definir tipos de roca y espesores útiles.
Inter-Rock, C. A. (2015), en su trabajo “CARACTERIZACIÓN PETROFÍSICA INTEGRADA DE LA FORMACIÓN LAJAS EN EL YACIMIENTO CUPEN MAHUIDA”, propone la incorporación de perfiles de Resonancia Magnética Nuclear (RMN) y su combinación con la metodología de Cluster Analysis para determinar tipos de roca a partir de los bins de porosidad de RMN.
Shepherd, M. (2009) “Volumetrics, oil field production geology– AAPG MEMOIR 91, p. 189-193”, define unidades de medida de hidrocarburos, y presenta como estimar petróleo (STOIIP) o gas en el sitio (GIIP), considerando factores importantes antes de comenzar la producción en el campo. Menciona tres métodos principales para el cálculo de volúmenes de reservorio y cómo validar esos volúmenes.
Snyder, R. H. (1971) “A review of the concepts and methodology of determining net pay” hace referencia a los métodos utilizados para determinar la producción neta en los yacimientos que contienen hidrocarburos. Entre otras cosas, define dos tipos de intervalos, uno descrito como “intervalo bruto” entre el tope y la base del horizonte geológico que se origina por el geólogo en la “correlación” y el otro que se define como “intervalo bruto de reservorio” y que está representado por la distancia (α-ῼ) existente entre el tope del primer
horizonte probable de almacenar fluido hasta la base del último, todos estos horizontes ubicados dentro del “intervalo bruto” descrito anteriormente. El “intervalo bruto de reservorio” contiene tanto roca reservorio, como roca no reservorio. Esto hace que el método propuesto por Snyder, empleado aún y desde hace muchos años, tienda a confundir a profesionales de otras ramas distintas a las geociencias. Como aspecto positivo, el método de Snyder aborda de manera específica conceptos de net pay (espesores netos hidrocarburíferos), las razones por las cuales es necesario calcularlo, y los pasos a seguir para obtener resultados más precisos al momento de definirlo, como por ejemplo análisis de núcleos y técnicas de mapeo para así luego calcular su volumen.
Worthington, P. (2005) “The applications of cutoffs in integrated reservoir studies”, abarca metodologías aplicables para net reservoir y métodos para determinar cortes petrofísicos en caso de agotamiento de fuentes primarias y recuperación secundaria por inyección de agua. Tomando en consideración que los cut-offs están condicionados dinámicamente, es decir están combinados a una referencia de parámetros claves de diagnóstico (diámetro de poro circular equivalente y puntos terminales de permeabilidades relativas), además de parámetros convencionales con el fin de obtener resultados más reales y certeros de un yacimiento. Presenta un flujo de trabajo para la aplicación de los parámetros de corte adecuados a los fines de la integración de estudios de reservorios, que no cubre todas las situaciones, pero sí ilustra cómo deben utilizarse los cortes en estudios de campo tanto determinísticos como estocásticos.
Worthington, P y Cosentino, L. (2005) “The role of cutoffs in integrated reservoir studies”, abordaron los diferentes enfoques para cuantificar los cut-offs, pero sin establecer uno de forma general para los yacimientos, ya que varían dependiendo del objetivo a estudiar en el mismo, tomando en cuenta el almacenamiento y el flujo de hidrocarburos, así como las permeabilidades relativas y la movilidad de los hidrocarburos. Se presta atención en la relación entre la identificación de los cut-offs y aquellos elementos claves de la práctica sistémica contemporánea de los estudios integrados de yacimiento, además de demostrar los principales beneficios de un conjunto de cortes petrofísicos debidamente condicionados, de modo que una empresa energética pueda realizar con mayor eficacia el valor del activo.
Metodología
A continuación, se presenta un flujograma de trabajo, el cual contempla cada uno de los pasos seguidos para el estudio realizado (Figura 4).
Figura 4. Flujograma de trabajo.
Tabla 2. Descripción de metodologías aplicadas más relevantes de acuerdo con diferentes autores.
Para el diseño del flujograma de trabajo se analizaron las metodologías más relevantes en la determinación de parámetros de corte y espesores útiles, además de los antecedentes relacionados con proyectos anteriores donde fue necesario determinar espesores útiles. Los antecedentes seleccionados para este trabajo fueron analizados con el objetivo de entender las distintas metodologías para el establecimiento de espesores útiles. Para ello se realizó una tabla comparativa que hace mención de los factores relevantes de cada método de acuerdo a diferentes autores (Tabla 2).
Para la incorporación de los pozos de estudio, registros y análisis de núcleos, se recopiló toda la información disponible de los pozos, de los distintos tipos de registros (eléctricos, acústicos, nucleares y litológicos), análisis de núcleos tanto convencionales como especiales (porosidad, permeabilidad, densidad de grano, saturación de fluidos, entre otros) y curvas evaluadas de las propiedades petrofísicas obtenidas (arcillosidad, porosidad, saturación de agua y permeabilidad con el fin de generar una base de datos lo más completa posible y con la menor incertidumbre.
Adecuación y QC de los datos: Antes de iniciar la interpretación de la información obtenida a partir de los registros, se debe realizar control de calidad para verificar los datos. En primera instancia, este control de calidad consiste en la validación del registro de campo y su comparación con datos digitales, para asegurar de que sean los mismos, haciendo énfasis en la sección donde se hará la interpretación. Posteriormente, se realizan varios pasos para la edición y adecuación de datos, entre los cuales podemos mencionar: ajustar las escalas de las curvas a utilizar, ajustar profundidades en las curvas, verificar si existen secciones repetidas en los registros, realizar correcciones ambientales, y eliminar y corregir valores nulos y anómalos.
Estimación de las propiedades petrofísicas mediante registros y calibración con análisis de laboratorio: para la realización de una evaluación cuantitativa de los perfiles de pozos es necesario conocer los parámetros petrofísicos de la formación, entre ellos: exponente de cementación (m), factor de tortuosidad (a) y exponente de saturación (n). Para la obtención de dichos parámetros, se utilizaron las mediciones resultantes de las propiedades eléctricas de núcleos en las formaciones “X”, “Y” y “Z”, específicamente las pruebas de factor de formación (FF) e índice de resistividad (IR).
Otro factor importante es la determinación de una función de temperatura en función de la profundidad, para lo cual se comenzó por determinar el gradiente geotérmico de las distintas formaciones consideradas para este estudio, a partir de la siguiente fórmula:
(1)
Donde:
GG: Gradiente geotérmico
Tf: temperatura de fondo (°F)
Ts: temperatura de superficie (°F)
Pf: profundidad final (pies)
Luego se estimó la temperatura de las formaciones estudiadas empleando la siguiente ecuación:
(2)
Donde:
Tform: temperatura de formación (°F)
GG: Gradiente geotérmico
Pform: profundidad de formación (pies)
Ts: temperatura de superficie (°F)
Propiedades petrofísicas: Arcillosidad (Vclay), porosidad (Φ), saturación de agua (Sw) y permeabilidad (k). Como parte del flujo de trabajo que realizó, la interpretación petrofísica se realizó de la siguiente manera:
Se calculó la arcillosidad, porosidad y saturación de agua para cada uno de los pozos. Para ello se escogió los modelos a utilizar para cada propiedad. Por ejemplo:
- La arcillosidad o volumen de arcilla (Vclay) se determinó a partir de la curva de Gamma Ray.
- La Porosidad (ɸ) se calculó utilizando el registro de Densidad.
- La saturación de agua (Sw) se determinó a partir de la ecuación de Archie.
Algunas de las formaciones consideradas para este estudio se contaban con datos fidedignos de salinidad. Para la determinación de Rw en donde no se contaba con información de salinidades, se utilizó el gráfico de “Pickett”, empleando las curvas de resistividad para definir las zonas de agua.
La permeabilidad (K) se determinó empleando la ecuación de Timur:
(3)
Este procedimiento se realizó para todos los pozos, no sin antes proceder a la calibración y verificación de los modelos petrofísicos. Para la calibración de los modelos petrofísicos se realizó una comparación y ajuste de los valores de las propiedades petrofísicas obtenidas mediante los registros con los datos de núcleos disponibles de dichas formaciones. Es importante mencionar que la calibración fue realizada con los datos obtenidos a condiciones de sobrecarga (NOBP).
Determinación de espesores útiles en función de los parámetros de corte: realizada la evaluación petrofísica, es necesario obtener los valores de los parámetros de corte para determinar las zonas correspondientes al Espesor Neto de Reservorio (Net Reservoir) y zonas que corresponden al Espesor Neto de Hidrocarburos (Net Pay). Dichos parámetros de corte para arcillosidad, porosidad y saturación de agua fueron calculados para cada uno de los métodos seleccionados en este estudio. (Tabla 2).
Snyder (1971): Snyder menciona en su estudio que normalmente se elige un indicador de permeabilidad y/o un valor de porosidad como valor de corte basado en un juicio intuitivo. Generalmente, este método es suficiente cuando hay una diferencia marcada entre los valores más altos y más bajos de las propiedades petrofísicas que definen lo que se consideraría reservorio y no reservorio, generalmente definido por contrastes obvios entre las curvas de registros que dan origen a propiedades como arcillosidad y porosidad principalmente. Se graficó la curva de arcillosidad con el fin de obtener el corte de Vclay utilizando el criterio mencionado de que todo lo que deflecte hacia la izquierda son arenas/reservorio y lo que deflecte a la derecha son arcillas/no reservorio. Con esto se obtiene el valor más aproximado de espesor bruto o también conocido como “arena total” y se excluyen las rocas con alto contenido de arcilla, como por ejemplo lutitas o algunas limolitas. Luego se visualizaron, mediante gráficos cruzados o “x-plots” e histogramas, los valores de porosidad correspondientes a la “arena total” obtenida con el corte de Vclay. Se tomó un valor mínimo siguiendo las recomendaciones de Snyder y se graficó este resultado como espesor neto de arena, es decir, las arenas que tengan propiedades de reservorio. Para el caso del espesor neto hidrocarburífero, aplicamos un corte de saturación de agua que habrían sido correlacionados con las calidades de los espesores de reservorio y las pruebas de pozos disponibles. Lo ideal en este caso sería poder soportar estos parámetros con datos de laboratorio. Este procedimiento se realizó para cada uno de los pozos del estudio y se tomó como referencia la figura 5, modificada del estudio original de Snyder.
Figura 5. Definición de espesores útiles (Modificado de: Snyder, 1971).
Figura 7. Histograma y “x-plots” para el cálculo del corte de phie. POZO-1.
Figura 8. Gráfico Vclay vs Vclay acumulado. POZO-1.
Figura 9. Gráfico Phie vs Phie acumulado. Scribner. POZO-1.
En las figuras 6 y 7 se muestra un resumen de los gráficos de este método aplicado al Pozo-1, a manera de ejemplo de lo que se realizó también en los otros pozos del estudio.
Figura 6. Gráficos para la obtención de parámetros de corte. POZO-1.
Scribner (1984): este método, originalmente creado para establecer parámetros de la función de permeabilidad, fue posteriormente adaptado y ampliamente utilizado para graficar valores de arcillosidad y porosidad contra el acumulado de sí mismas, obteniéndose un gráfico en el que el valor del parámetro de corte corresponde a la intersección entre las rectas tangentes a las dos primeras tendencias que se originan en dicho gráfico.
En este caso particular se utilizan las figuras 8 y 9 para mostrar los resultados obtenidos en el Pozo-1 a manera de ejemplo de lo que se hizo también en los otros pozos del estudio. Para realizar los cálculos de los parámetros de corte en este método se exportaron las curvas de arcillosidad y porosidad efectiva en formato. LAS para cada uno de los pozos, luego se cargaron en Excel con la finalidad de calcular su acumulado y realizar los gráficos de Vclay vs Vclay acumulado y Phie vs Phie acumulado.
En este método de Scribner, o más bien adaptación del método de Scribner para generar parámetros de corte de permeabilidad, únicamente se obtienen cortes de arcillosidad y porosidad efectiva, ya que su trabajo tampoco indica el uso de esta metodología para obtener valores de corte de saturación de agua. Sin embargo, como valor agregado para este método y para realizar una comparación completa de todos los espesores útiles se utilizaron los valores de corte de saturación de agua proporcionados por un método similar al del método de Snyder adaptado por la empresa Inter-Rock C. A. y así obtener el net pay y no solamente el net reservoir.
Worthington y Cosentino (2005): En su diagrama de flujo mencionan que para determinar los parámetros de corte debe realizarse de una forma sinérgica y deben estar condicionados dinámicamente. Este método cubre la doble publicación hecha por Worthington en 2005, haciendo énfasis en el trabajo compartido con Luca Cosentino donde el condicionamiento dinámico incluye una variante representada por el concepto de “movilidad”. Para ello identifican un límite que debe estar vinculado a un parámetro dinámico, como la permeabilidad absoluta; el diámetro de poro circular equivalente; la movilidad, que está relacionada con la viscosidad del fluido.
En caso de contar con data especial de laboratorio, también podría relacionarse con la presión capilar; la saturación de agua irreducible; o la permeabilidad relativa extrapolada a partir de los “end points”, esto con el fin de identificar un punto de cruce ubicado en los que se consideraría la vida útil del reservorio y que puedan expresarse en términos de un valor de corte. Una vez identificado este límite, se puede relacionar con y quizás.
En la figura 10, se observa que el límite relacionado con el parámetro dinámico a utilizar en este ejemplo está representado por la movilidad, que de acuerdo al trabajo realizado por los autores se define como la relación entre la permeabilidad absoluta y la viscosidad.
Figura 10. Generación esquemática de cortes sinérgicos V_(sh,) ϕ y S_w para el agotamiento primario del yacimiento, donde la movilidad y la saturación de agua están distribuidas normalmente. (Worthington y Cosentino, 2005).
Worthington y Cosentino mencionan en su trabajo que para dar entrada a las reservas de petróleo y de gas se propone que los parámetros de corte se vinculen a un valor límite de movilidad. La viscosidad de un gas es al menos un orden de magnitud menor que la de petróleo ligero y la magnitud de esta disparidad probablemente explica la diferencia entre los valores tradicionales de corte de permeabilidad de 0,1 mD para el gas y 1 mD para el petróleo. Dicho lo anterior, es importante mencionar que para este método los parámetros de corte se determinaron según la naturaleza de cada yacimiento. Los autores de este método advierten que, aunque no existe una definición universal de espesores útiles, este podría ser un método adaptable a casi todas las naturalezas de reservorios, por lo que ofrecen una serie de términos, que ayudan a tener una idea más clara de cómo ilustrar el mecanismo de obtención de los espesores útiles. Figura 11.
Figura 11. Interrelación esquemática secuencial de parámetros de corte para la determinación de espesores útiles. (Worthington y Cosentino, 2005).
a. Espesor Total de Roca (G): incluye toda la roca dentro del intervalo a evaluar.
b. Espesor Neto (S): se refiere a aquella roca que pudiera tener características de reservorio. El criterio para determinarlo se expresa en función de un valor de corte de Volumen de Arcilla (Vclay).
c. Espesor Neto de Reservorio (R): se refiere a aquella parte del Espesor Neto que tiene propiedades útiles de reservorio. La condición para determinarlo se expresa en función de un valor de corte de porosidad (Phi).
d. Espesor Neto de Hidrocarburos (P): se refiere a los intervalos netos de reservorio que contienen cantidades significativas de hidrocarburos. Se determina usando como criterio un corte de saturación de agua (Sw).
Para el POZO-1, a pesar de que contiene petróleo, no se utilizó la permeabilidad de 1 mD, por ser este un yacimiento que contiene hidrocarburo pesado, por lo que el valor de permeabilidad tuvo que ser mucho mayor. Para obtener este valor según los lineamientos de los autores de este método, se verificaron los valores y análisis de PVT, donde la viscosidad era de 23 cp y la movilidad de 0,65 mD/cp., dando como resultado una permeabilidad de corte de 15 mD. La fórmula para realizar el cálculo de dicha permeabilidad fue la siguiente:
(4)
Donde:
M: movilidad (mD/cp)
k: permeabilidad de corte (mD)
μ: viscosidad (cp)
Es importante mencionar que para el POZO-2, por ser un yacimiento que contiene gas, se asumió que el valor de permeabilidad es de 0,1 mD, mientras que para el POZO-3 se utilizó un límite de permeabilidad de 1 mD, que según los autores corresponde a un yacimiento de petróleo liviano. A continuación, los gráficos del Pozo-1 mostrando la determinación de parámetros de corte según el Método de Worthington y Cosentino, (Figuras 12,13 y 14). Estos gráficos fueron también generados para los Pozos 2 y 3.
Figura 12. “X-plot” permeabilidad vs porosidad. POZO-1.
Figura 13. “X-plot” Vclay vs porosidad. POZO-1.
Figura 14. “X-plot” Sw vs porosidad. POZO-1.
Figura 15. Mapa del espesor del volumen poroso de hidrocarburos que muestra la distribución de hidrocarburos en un reservorio según método de Shepherd.
Shepherd (2009): en su publicación “VOLUMETRICS, OIL FIELD PRODUCTION GEOLOGY. AAPG MEMOIR 91”, menciona que los volúmenes o espesores intermedios pueden relacionarse con las propiedades medias de la roca en una zona petrolífera de la siguiente manera:
- Volumen neto de la roca (NRV) = Volumen bruto de la roca (GRV) x promedio del net to gross
- Volumen neto del poro (NPV) = Volumen neto de la roca (NRV) x porosidad promedio.
- Volumen poroso de hidrocarburos (HPV)= Volumen neto del poro (NPV) x saturación de petróleo promedio.
- Petróleo original en sitio (STOIIP) = Volumen poroso de hidrocarburos (HPV) x 1/Bo (factor de volumen de formación para petróleo).
A su vez, estas propiedades de las rocas y de los fluidos pueden calcularse para el tramo de hidrocarburos de la siguiente manera:
- Net to gross = Volumen neto de la roca (NRV) Volumen bruto de la roca (GRV).
- Porosidad = Volumen neto del poro (NPV) Volumen neto de la roca (NRV).
- Saturación de petróleo = Volumen poroso de hidrocarburos (HPV) Volumen neto del poro (NPV).
- Factor volumen de formación = Volumen poroso de hidrocarburos (HPV) Petróleo original en el sitio (STOIIP).
Con ayuda de un mapa del espesor hidrocarburífero se visualiza el método. (Figura 15).
La figura 15 muestra que el espesor del volumen poroso de hidrocarburo es igual al espesor bruto x net-to-gross x porosidad x saturación de petróleo (Shepherd, 2009).
Shepherd, en su trabajo, no menciona cómo determinar los parámetros de corte, mayormente se enfoca en la determinación de los espesores útiles. Sin embargo, para este método los parámetros de corte se obtuvieron siguiendo los conceptos fundamentales por él esgrimidos en relación con los espesores porosos saturados de hidrocarburos mediante los gráficos que tradicionalmente emplean
la variable PhieSo para la determinación de los parámetros de corte convencionales.
Para obtener el corte de arcillosidad se realizaron gráficos de PhieSo vs Vclay. PhieSo corresponde al Volumen Poroso ocupado por petróleo y está definido como el porcentaje de la porosidad que se encuentra saturada por petróleo y que es igual a PHIE*(1-BVW) (Inter-Rock, C. A., 2011).
Luego de realizados los cálculos de PhieSo se emplearon gráficos cruzados XY para determinar los parámetros de corte en los pozos destinados para este estudio. En estos gráficos se muestra la variable PhieSo en el eje “y” y el parámetro de corte a determinar en el eje “x”. El procedimiento se lleva a cabo trazando una línea envolvente que separa una nube de puntos considerados como parte del espacio poroso que contiene hidrocarburos de aquellos que no lo contiene, la cual intercepta el eje X y así obtener el corte de la propiedad que se quiere determinar. La figura 16 muestra la determinación del corte de Vcl para el Pozo-1 el cual resulto en 27 % y a partir de ese valor, la figura 17 muestra la determinación del parámetro de corte de porosidad, el cual resultó en 10 %.
Figura 16. Gráfico PhieSo vs Vclay para obtener corte de arcillosidad. POZO-1.
Figura 17. Gráfico de Phie vs Vclay para obtener corte de porosidad. POZO-1.
Los valores de corte para saturación de agua en este método fueron determinados de la misma manera; relacionando el espacio poroso saturado por hidrocarburo versus la saturación de agua (Sw), de manera de poder visualizar el valor máximo de esta variable diferenciando el espacio poroso con hidrocarburo del que contenía altas saturaciones de agua, el cual resulto en 60 % de Swmax para el POZO-1.
Finalmente, para la determinación de espesores útiles y sumarios de evaluación petrofísica, se completaron los valores de parámetros de corte los cuales se ingresan en el módulo de sumarios de evaluación de la plataforma o software disponible para hacer la corrida correspondiente a los pozos objeto de estudio y obtener los espesores netos en una tabla que permita su visualización, análisis y comparación.
Para ello, se escogen las curvas anteriormente obtenidas en la evaluación petrofísica sobre las cuales se aplicarán los parámetros de corte según los criterios que priven para cada caso o método utilizado. Estos criterios se establecen para diferenciar los distintos tipos de espesores, como por ejemplo los espesores de reservorio neto o “net Reservoir” y los espesores netos hidrocarburíferos o “net pay”. Los parámetros de corte se aplicarán según los símbolos matemáticos comparativos , que se identifican con “< =” (menor o igual que), “> =” (mayor o igual que) o “=” (igual que), indicando el método de parámetros de corte utilizado, los cuales representaran valores promedio en porcentaje (%).
Las curvas seleccionadas para la aplicación de los parámetros de corte fueron: Vclay (volumen de arcilla), Phie (porosidad) y Sw (saturación de agua).
A continuación, se muestra la tabla 3 con los parámetros de corte calculados para cada uno de los métodos seleccionados en este trabajo:
Tabla 3. Valores determinados de parámetros de corte por métodos utilizados.
Tabla 4. Métodos de determinación de espesores útiles según su Autor, fecha de publicación y propósito.
Como se observa en la tabla 3, los valores de los parámetros de corte varían según sea el método empleado, sin embargo, se puede notar cierta consistencia en la mayoría de ellos, estos valores que se mantienen dentro de un margen o rango similar para cada uno de los pozos.
Discusión de resultados
De acuerdo con la metodología descrita anteriormente, en este capítulo se mostrarán y se discutirán los resultados obtenidos de los espesores útiles mediante análisis comparativos de los distintos escenarios considerados.
Los Cálculos y algoritmos involucrados en los análisis de parámetros de corte para el establecimiento de espesores útiles se alimentaron principalmente de los siguientes datos de entrada: Vclay (arcillosidad), Phie (porosidad) y Sw (saturación de agua). Cada Investigador dentro de su método emplea algoritmos, gráficos y consideraciones distintas que fundamentan su razonamiento para la determinación de espesores útiles. (Tabla 4).
En la tabla 4 se muestran los métodos considerados como los más relevantes o más comúnmente empleados en la determinación de espesores útiles. Entre los autores de dichos métodos tenemos a: Snyder, Scribner, Worthington & Cosentino y Shepherd. Esta tabla proporciona información del año de publicación, el propósito y los parámetros de corte a utilizar para cada método.
Snyder (1971) fue uno de los primeros en mencionar que para seleccionar y definir los parámetros de corte a utilizar en un yacimiento era importante tomar en consideración el propósito de los espesores útiles a calcular, sin embargo, analizando otros antecedentes no siempre era el caso, en algunas ocasiones utilizaban criterios simples, haciendo que dicho punto pasara desapercibido.
Otro punto relevante que se observó en las metodologías estudiadas es que para cada autor los espesores e intervalos en un yacimiento tenían distintas nomenclaturas, así como también, distintas definiciones. A continuación, en la tabla 5, modificada de Worthington y Cosentino (2005) donde se puede observar lo mencionado anteriormente.
Tabla 5. Comparación de las nomenclaturas utilizadas originalmente por los autores de los métodos más relevantes de determinación de espesores útiles. (Modificada de Worthington & Cosentino – 2005).
Con respecto a las curvas evaluadas, los tres pozos presentaron curvas similares y suficientes para realizar los cálculos necesarios en este estudio, así como también los datos de núcleos disponibles en los tres pozos para la calibración de las propiedades petrofísicas obtenidas, como lo son los datos de porosidad y permeabilidad de los yacimientos. Dicho esto, los valores obtenidos de los parámetros petrofísicos en este estudio se observan que para el POZO-2 y POZO-3 son areniscas consolidadas, sin embargo, para el POZO-1 los parámetros petrofísicos obedecían a una arenisca menos consolidada. Luego de proceder a la Calibración y verificación de las propiedades petrofísicas, observamos diversos rangos de las mismas en los pozos de estudio.
En el POZO-1, se observa que los valores de porosidad son mayores de 20 % y tiene permeabilidades máximas aproximadas a los 4.000 mD, estos valores se consideran altos en comparación a los otros yacimientos evaluados. A partir de los 11.497 pies presenta una caída gradual en la resistividad y la saturación de agua llega a valores cercanos a 100 %, quizás a consecuencia de una zona de transición.
El POZO-2 por ser un yacimiento de areniscas apretadas sus valores son bajos en comparación a los pozos 1 y 3. Los valores de porosidad no son mayores del 12 % y sus permeabilidades oscilan entre 0,01 a 0,2 mD.
En el POZO-3 la zona de mayor interés de este yacimiento está representada por porosidades de hasta a un 10 % y las permeabilidades entre < 1 mD a 10 mD.
A continuación, a manera de ejemplo, se muestra un gráfico de espesores útiles (Pozo-1/ método Snyder), realizado en cada uno de los métodos/Pozos de estudio. Figura 18.
Figura 18. Gráfico de espesores útiles. Método Snyder. POZO-1.
En la tabla 6 se tiene un resumen de los resultados obtenidos de los espesores útiles de todos los métodos en todos los pozos.
Tabla 6. Resultados finales de espesores útiles a partir de los métodos desarrollados en este estudio.
De acuerdo a los resultados obtenidos que se muestran en la tabla 6, se procedió a realizar una comparación de todos los métodos. A grandes rasgos, se observó que el método propuesto por Shepherd obtuvo resultados más favorables en cuanto a los espesores útiles en comparación con el resto de los métodos. Sin embargo, la mayoría de los valores eran consistente entre los métodos aplicados a cada pozo y no variaron significativamente con contadas excepciones. A continuación, se realiza un análisis más detallado de los resultados para cada uno de los yacimientos/pozos:
POZO-1/Yacimiento “X”: se observó que, en este caso, no hubo una variación significativa en el espesor de reservorio. El valor más alto para reservorio neto fue de 264.5 pies, el cual se obtuvo por el método de Worthington & Cosentino, así mismo se comparó el valor de reservorio neto hidrocarburífero, donde se logra observar que, en este caso, también el método de Worthington y Cosentino obtuvo el mayor valor de 175.25 pies de espesor para el yacimiento “X”. Muy probablemente, esto se debe a que este método funciona muy bien para yacimientos con permeabilidades altas.
POZO-2/Yacimiento “Y”: los espesores útiles tanto para el espesor neto de reservorio y espesor neto hidrocarburífero presentaron diferencias considerables. Como se mencionó anteriormente, este yacimiento es de areniscas apretadas. Se observa que el método que presentó resultados bastante favorables en comparación de los demás fue el método de Shepherd con espesores de reservorio neto e hidrocarburífero de 799 y 780.5 pies respectivamente, mientras que el método menos favorable en este caso fue el de Worthington y Cosentino con 218,5 pies de reservorio neto y 208 pies de reservorio neto con hidrocarburo; de esta manera se confirmó que este último método quizás no sea el más adecuado para pozos con permeabilidades y porosidades bajas.
POZO-3/Yacimiento “Z”: En este caso, dos métodos resultaron ser los más favorables, tanto para los espesores netos como los espesores netos hidrocarburíferos, al coincidir con valores de 40 pies y 38,5 respectivamente cuando se aplicaron los cortes al yacimiento “Z”. Estos métodos resultaron ser los de Snyder y Shepherd, lo cual demuestra “consistencia” en los resultados y modelos petrofísicos robustos. El método menos favorable en este caso fue el de Worthington y Cosentino con tan solo 7,5 pies de reservorio neto y 7 pies de espesor neto hidrocarburífero, quizás por estar en presencia de un yacimiento con arenas consolidadas de poca permeabilidad, aún en presencia de hidrocarburos gaseosos o petróleos livianos.
Es importante recalcar que la aplicación de parámetros de corte “fijos” favorece la omisión de algunos intervalos cuya saturación de agua es variable debido a la variación de calidad de la roca almacén, en la cual y en muchos casos, la Sw representa el agua irreducible o no movible. Esto deja abierta la discusión para que en futuros estudios se consideren parámetros de corte “variables” en cuanto a la Sw ya que está íntimamente relacionada a la calidad de los reservorios o los “Tipos de Roca”. Existen varios tipos de roca, rocas buenas, regulares y de mala calidad, por ello la aplicación de un corte de Sw variable lograría aportar resultados más realistas y más favorables que los otros métodos los cuales utilizan un corte de Sw fijo.
Conclusiones
1. Se generó una base de datos de tres pozos seleccionados para este trabajo, como resultado de la recopilación, organización y clasificación de los datos de registros de pozos y análisis de núcleos, tanto convencionales como especiales, que sirvió de base para la verificación de propiedades petrofísicas y para el cálculo de parámetros de corte de las formaciones “X”, “Y” y “Z” designadas para este estudio.
2. La revisión previa de los antecedentes referentes al tema de selección de parámetros de corte para el establecimiento de espesores útiles, permitió facilitar la metodología para determinar los parámetros de corte a partir de los métodos relevantes existentes seleccionados para el estudio.
3. Se observó que a pesar de existe suficiente información para determinar los espesores útiles, no es el caso para el caso de determinar parámetros de corte; cada método presenta un enfoque distinto y tiende a ser relativamente subjetivo, ya que depende del propósito del estudio y de las características del yacimiento a evaluar.
4. En líneas generales, los métodos de Snyder y Shepherd generaron valores más altos de espesores útiles en comparación al resto de los métodos con la excepción del yacimiento “X” del POZO-1 donde el método de Worthington y Cosentino superó el conteo de pies de espesores netos y netos hidrocarburíferos de los otros métodos; muy probablemente debido a la consideración de variables de permeabilidad y movilidad, las cuales aportan un peso significativo en yacimientos de alta calidad de almacenamiento y flujo. Se observó que para cada yacimiento no todos los métodos son favorables. Las características de los yacimientos generalmente varían y no es recomendable que los parámetros de corte sean fijos, como generalmente sucede en la práctica. Sobre todo, con relación a la saturación de agua (Sw).
5. Cada metodología para determinar espesores útiles presenta sus ventajas. Snyder, aunque subjetivo, ofrece resultados gráficos aceptables; Scribner utiliza propiedades acumuladas para calcular parámetros de corte; Shepherd considera el volumen poroso de petróleo apoyado por los gráficos PhiSo; Worthington y Cosentino incorporan la permeabilidad y movilidad como propiedades que relacionan el carácter dinámico de los métodos para determinar los parámetros de corte.
6. Este trabajo logra aportar algo con lo que no se cuenta comúnmente en la práctica, que es la manera de implementar no solo diversos métodos de espesores útiles sino también una guía para determinar los parámetros de corte de cada método en función de las características de los yacimientos.
Recomendaciones
1. Asegurarse de preparar y adecuar los datos de los pozos objeto de estudio para garantizar la calidad de los mismos y con ello de los resultados.
2. Aprovechar la disponibilidad de datos de núcleo para calibrar las propiedades petrofísicas y mejorar la representatividad de las propiedades de los yacimientos y disminuir la incertidumbre al momento de calcular los parámetros de corte y espesores útiles.
3. Definir con antelación el propósito de calcular los parámetros de corte, ya que los métodos y los parámetros a determinar dependen en gran medida de los objetivos del estudio. Por ejemplo: cálculo de POES o evaluación del potencial de hidrocarburos remanentes para recobros secundarios.
4. El uso de parámetros de corte fijos no garantiza la óptima determinación de espesores útiles, por lo que se recomienda la implementación de cortes variables sobre todo de la saturación de agua (Sw) en virtud de sus variaciones según el tipo de roca o calidad de reservorio.
5. Procurar la disponibilidad de análisis especiales de núcleos (SCAL) para poder aplicar cortes de tipo variable implementando el carácter dinámico como lo recomiendan Worthington & Cosentino en su trabajo.
6. Definir tipos de roca existentes en los yacimientos, ya que nos permite priorizar la calidad de los reservorios presentes en función de los valores de garganta poral más representativos facilitando la eventual aplicación de los parámetros de corte variables.
Referencias bibliográficas
Scribner, M.K. & & Roy Engineer. HP-41CV Analyzes Porosity, Permeability Data. Petroleum Engineer International. 1984.
Baker Hughes (2004). BAKER ATLAS. LOG INTERPRETATION CHARTS. Baker Hughes Incorporated.
Barbato, R., Quaglia, A., Porras, J., Panesso, R., & Soria, J. / Inter-Rock C.A. (2006). LOG DATA PREPARATION. A REAL OPPORTUNITY TO SAVE MONEY AND TIME! Cancún, México. Congreso Mexicano del Petróleo, SPE Paper 104072.
Biguelow, E. (1995). INTRODUCTION TO WIRELINE LOG ANALYSIS. Houston, Texas. Editorial Western Atlas.
Crain, E. (2018). CRAIN’S PETROPHYSICAL HANDBOOK. 26 de febrero de 2020, [https://www.spec2000.net/00-fineprint.htm].
Ellis, D. & Singer, J. (2007). WELL LOGGING FOR EARTH SCIENTISTS. Dordrecht, Holanda. Editorial Springer.
García, R. & Valera Y. (2014). INTÉRPRETE DE PROGRAMAS DE USUARIO PARA EL CÁLCULO DE PARÁMETROS PETROFÍSICOS. 25 de febrero de 2020, [http://scielo.sld.cu/scielo.php?script=sci_arttext&pid=S2227-18992014000200008].
Greve, J. (2020) NEURALOG- AUTOMATED WELL LOG DIGITIZING SOFTWARE. [https://www.neuralog.com/pages/NeuraLog].
Krygowsky, D. (2003). GUIDE TO PETROPHYSICAL INTERPRETATION. Austin, Texas.
López, M. (2013). PETROFÍSICA Y EFECTOS AMBIENTALES EN REGISTROS GEOFÍSICOS DE POZOS. Ciudad Universitaria, México. Universidad Nacional Autónoma de México.
Rider, M. (1996). THE GEOLOGICAL INTERPRETATION OF WELL LOGS (2DA EDICION). Sutherland, Escocia. Editorial Rider-French Consulting Ltd.
Schlumberger. (1999). PRINCIPIOS/APLICACIONES DE LA INTERPRETACIÓN DE REGISTROS. Schlumberger Educational Services.
Schlumberger. (1972). INTERPRETACIÓN DE PERFILES. Nueva York: Schlumberger Limited, 19-80.
Theys, P. (1999). LOG DATA ACQUISITION AND QUALITY CONTROL (2DA EDICION). Paris, Francia. Editorial Technip.
Inter-Rock, C. A (2018), “CARACTERIZACIÓN PETROFÍSICA DEL YACIMIENTO LINDERO ATRAVESADO, GRUPO CUYO, CUENCA NEUQUINA”.
Inter-Rock, C. A (2015), “CARACTERIZACIÓN PETROFÍSICA INTEGRADA DE LA FORMACIÓN LAJAS EN EL YACIMIENTO CUPEN MAHUIDA”.
Inter Rock, C. A. WORKSHOP “Determinación de espesores Útiles”. Jornadas Técnicas Repsol. Puerto la Cruz. 2006.
Shepherd, M. (2009) “VOLUMETRICS, OIL FIELD PRODUCTION GEOLOGY– AAPG MEMOIR 91, p. 189-193”.
Snyder, R. H. (1971) “A REVIEW OF THE CONCEPTS AND METHODOLOGY OF DETERMINING NET PAY”.
Worthington, P. (2005) “THE APPLICATIONS OF CUTOFFS IN INTEGRATED RESERVOIR STUDIES”, SPE. Annual Technical Conference. Dallas.
Worthington, P., Cosentino, L. (2005) “THE ROLE OF CUTOFFS IN INTEGRATED RESERVOIR STUDIES”, SPE. Annual Conference. Denver.
The author(s) declare(s) that she/he/they has/have no conflict of interest related to hers/his/their publication(s), furthermore, the research reported in the article was carried out following ethical standards, likewise, the data used in the studies can be requested from the author(s), in the same way, all authors have contributed equally to this work, finally, we have read and understood the Declaration of Ethics and Malpractices.