BENCHMARKS
Benchmarks
Resultados de Converge comparados con soluciones de referencia publicadas, con las ecuaciones publicadas y con cálculos a mano. Cada caso da la referencia, la tolerancia y el resultado, y dice con claridad lo que aún no se ha validado. Los casos se ejecutan como pruebas con nombre en la batería de pruebas de Converge.
RESUMEN
Dónde están los motores
Una auditoría independiente el 15 de octubre de 2026 recalculó 30 resultados de Converge a mano o con código independiente. En los casos sintéticos, los datos se construyeron a partir de respuestas conocidas, así que la verdad es exacta. 25 se aprobaron a la primera. Los cinco fallos estaban en valores por defecto, puntos de referencia y bandas de incertidumbre, no en la física central. Cuatro están corregidos y ahora se aprueban. El quinto es una corrección de método cuyo único caso de prueba sigue en el límite, y lo decimos más abajo.
- Casos de validación
- 30
- Aprobados en la batería de pruebas
- 29
- En el límite, declarado
- 1
- Tolerancia típica
- 0.05–3%
Los valores están en unidades de campo, como en las referencias publicadas. Estado al 15 de octubre de 2026.
SIMULACIÓN DE YACIMIENTOS
SPE1: el primer proyecto de soluciones comparativas de la SPE
Odeh, A. S. (1981), «Comparison of Solutions to a Three-Dimensional Black-Oil Reservoir Simulation Problem», Journal of Petroleum Technology 33(1), SPE 9723. Una malla de 10 × 10 × 3, un inyector de gas y un productor durante diez años. El caso viene integrado en Converge.
| Magnitud | Comportamiento publicado | Converge | Tolerancia | Estado |
|---|---|---|---|---|
| Meseta de caudal de petróleo | 20,000 STB/d | 20,000 STB/d | 1 STB/d | Aprobado |
| El productor pasa de control por caudal a control por presión | Cerca del día 1.000 | Día 973 | Día 800 a 1.300 | Aprobado |
| Pico de la presión media del yacimiento | Cerca de 6.400 psia | 6,423 psia | 6.100 a 6.700 psia | Aprobado |
| Caudal de petróleo del campo, RGP, presiones de fondo de los pozos, presiones de las celdas | Corrida de referencia publicada de un simulador estándar de la industria, 11 tiempos de reporte hasta el día 3.650 | 1% · 2.5% · 2% · 3% | Mayor diferencia indicada | Aprobado |
| Cierre del balance de masa | Masa conservada en cada paso | 9.3 × 10⁻⁸ | Por debajo de 10⁻⁶ del volumen poroso | Aprobado |
La corrida tiene 132 pasos de tiempo y tarda unos 1,4 segundos en un solo núcleo de portátil. Un caso Buckley-Leverett unidimensional sitúa el frente de agua a menos de tres celdas de la construcción analítica de Welge.
V1–V30
Los 30 casos de validación auditados
Cada caso es ahora una prueba con nombre. «Corregido» significa que el caso falló en la auditoría, se corrigió la causa y la prueba ahora se aprueba.
| Caso | Cálculo | Referencia | Tolerancia | Estado |
|---|---|---|---|---|
| V1 | Punto de burbuja de Standing, 35 °API, gravedad del gas 0,75, 650 scf/STB, 200 °F | 2.748 psia (cálculo manual) | ±1.6 psi | Aprobado |
| V2 | Punto de burbuja de Standing, 32 °API, 0,72, 520 scf/STB, 205 °F | 2,601 psia | ±1.6 psi | Aprobado |
| V3 | Factor volumétrico del petróleo de Standing en el punto de burbuja, ambos fluidos | 1,375 y 1,3037 rb/STB | 0.0005 rb/STB | Aprobado |
| V4 | Viscosidad del petróleo vivo de Beggs-Robinson en el punto de burbuja, ambos fluidos | 0,4927 y 0,6185 cp | 0.1% | Aprobado |
| V5 | Densidad del petróleo en el punto de burbuja | 695.0 kg/m³ | 0.1% | Aprobado |
| V6 | Factor Z, Dranchuk-Abou-Kassem con Sutton, a cuatro presiones | Ver la tabla del factor Z | 0.05% | Aprobado |
| V7 | Corrección por presión de la viscosidad del agua de McCain a 5.000 psia | +28% | ±0,5 puntos | Aprobado |
| V8 | La división de la fracción pesada conserva el peso molecular del C7+ | 218.0 g/mol | 10⁻⁶ g/mol | Aprobado |
| V9 | Factor volumétrico del gas a 3.000 psia | 0.000983 rb/scf | 10⁻⁶ rb/scf | Aprobado |
| V10 | Petróleo original in situ por balance de materiales, yacimiento volumétrico construido a partir de un N conocido, sin acuífero | 45.00 MMSTB | 1% | Aprobado |
| V11 | El mismo yacimiento con un acuífero que no tiene | No debe dar un N erróneo | Ajuste bloqueado, con los motivos | Corregido |
| V12 | Ajuste histórico por balance de materiales en el yacimiento de ejemplo | 150 MMSTB | 3% | Aprobado |
| V13 | La banda P90–P10 del ajuste histórico contiene la verdad | 150 MMSTB dentro de la banda | Cobertura de al menos el 80 % en 10 historiales sintéticos | Corregido |
| V14 | Petróleo original in situ por balance de materiales dinámico | Verdad sintética conocida | 1% | Aprobado |
| V15 | Ajuste de Arps a un pozo sintético, qi 3.500 STB/d, Di 60 %/año, b 0,70 | Los parámetros generadores | Datos exactos: 3 decimales. Ruido del 3 %: 5 % en qi | Aprobado |
| V16 | EUR de Arps con un límite de 30 años | 4.78 MMSTB | 0.005 MMSTB | Aprobado |
| V17 | EUR exponencial y armónico | Formas cerradas: 1.980 y 9.210 STB | 0,5 STB y 5 STB | Aprobado |
| V18 | La banda P90–P10 del bootstrap de declinación contiene la verdad | Cobertura de al menos el 80 % en 12 pozos sintéticos | Observado: 10 de 12 | En el límite |
| V19 | Pendiente y permeabilidad de Horner, restauración sintética con k = 180 md | 26.43 psi/cycle, 180 md | 3% | Aprobado |
| V20 | Presión extrapolada de Horner p* | 3,650 psia | 2 psi | Aprobado |
| V21 | Daño (skin) de la misma restauración, Horner y log-log | +4.5 | ±0.3 | Corregido |
| V22 | Permeabilidad de flujo radial (log-log) | 180 md | 3% | Aprobado |
| V23 | Gradiente hidrostático en el modelo de pozo | 18,06 lbm/ft³ da 0,1254 psi/ft | 10⁻⁴ psi/ft | Aprobado |
| V24 | Conversión de unidades al importar, 40 m³/d y 25 barg | 251,59 bbl/d y 377,29 psia | Exacto con 2 decimales | Aprobado |
| V25 | Regalía e impuesto sobre la renta, cada año de un caso de 10 años | Cálculo manual independiente | 10⁻⁶ | Aprobado |
| V26 | VPN al 10 %, descuento a mitad de periodo | Suma independiente de los flujos de caja descontados | 10⁻³ | Aprobado |
| V27 | TIR | Raíz de la ecuación del VPN | 10⁻⁹ | Aprobado |
| V28 | SPE1 (Odeh 1981) | Comportamiento publicado | Ver la tabla de SPE1 | Aprobado |
| V29 | Balance de masa de la red de superficie en cada nodo | Masa entrante igual a masa saliente | 10⁻⁹ | Aprobado |
| V30 | Contabilidad del gas de la red con gas lift | Gas de formación más gas de levantamiento igual a gas que llega; gas de venta igual a gas de formación | 10⁻⁶ | Corregido |
V18: el método bootstrap se corrigió y cumple su objetivo de cobertura del 80 % en 12 pozos sintéticos. En el único pozo de ejemplo que usó la auditoría, el EUR verdadero sigue justo por encima del P10, algo que una banda del 80 % hará aproximadamente una vez de cada diez. No ajustamos la semilla ni ensanchamos la banda para forzarlo a entrar. V25–V27: los datos del caso del auditor no se publicaron, así que en su lugar se verifica un caso equivalente de 10 años construido a mano.
FLUIDOS
Correlaciones PVT frente a las ecuaciones publicadas
Cada correlación se codifica a partir del artículo original y se verifica con un cálculo a mano de la misma ecuación. Resultados para el fluido V1: 35 °API, gravedad del gas 0,75, 650 scf/STB, 200 °F.
| Propiedad | Fuente publicada | Cálculo manual | Converge | Estado |
|---|---|---|---|---|
| Punto de burbuja | Standing (1947) | 2,748 psia | 2,747.8 psia | Aprobado |
| Bo en el punto de burbuja | Standing (1947) | 1.375 rb/STB | 1.3749 rb/STB | Aprobado |
| Viscosidad del petróleo vivo | Beggs y Robinson (1975) | 0.4927 cp | 0.4924 cp | Aprobado |
| Punto de burbuja y Bo | Vasquez y Beggs (1980), Glasø (1980), Petrosky y Farshad (1993), Al-Marhoun (1988) | Cálculo manual para cada conjunto | Coincide; el punto de burbuja devuelve la RGP de entrada | Aprobado |
| Corrección por presión de la viscosidad del agua | McCain (1990) | +28% | +28% | Aprobado |
| Corrección por gas ácido | Wichert y Aziz (1972) | Cálculo manual de ε | Coincide | Aprobado |
El motor de ecuación de estado se verifica por separado. Los puntos de burbuja y de rocío del petróleo SPE5 y de un gas condensado, en Peng-Robinson y SRK, coinciden con una implementación independiente de código abierto dentro del 0,01 %, y las composiciones de flash bifásico dentro de 1 × 10⁻⁵. Las presiones de vapor de componentes puros quedan a menos del 3 % de los datos del NIST.
GAS
Factor Z: Dranchuk y Abou-Kassem
Dranchuk, P. M. y Abou-Kassem, J. H. (1975), Journal of Canadian Petroleum Technology 14(3), con las propiedades pseudocríticas de Sutton (1985). Gravedad del gas 0,7, 200 °F.
| Presión | Cálculo manual | Converge | Diferencia | Estado |
|---|---|---|---|---|
| 1,000 psia | 0.9215 | 0.9213 | 0.02% | Aprobado |
| 3,000 psia | 0.8883 | 0.8880 | 0.03% | Aprobado |
| 5,000 psia | 0.9999 | 0.9997 | 0.02% | Aprobado |
| 10,000 psia | 1.4308 | 1.4310 | 0.01% | Aprobado |
Tolerancia del 0,05 %. El factor volumétrico del gas a 3.000 psia es 0,000983 rb/scf, como en V9.
PRUEBA DE POZO
Horner y log-log en una restauración sintética
Un periodo de flujo de 48 horas a 800 STB/d y luego 72 horas de cierre, generados con la solución de fuente lineal y superposición. Respuesta conocida: k = 180 md, daño +4,5, presión inicial 3.650 psia. Espesor neto 32,8 ft, porosidad 0,20, compresibilidad total 1 × 10⁻⁵ 1/psi, viscosidad 1 cp, factor volumétrico 1,2 rb/STB.
| Resultado | Respuesta conocida | Converge | Tolerancia | Estado |
|---|---|---|---|---|
| Pendiente de Horner | 26.43 psi/cycle | 26.46 psi/cycle | 3% | Aprobado |
| Permeabilidad de Horner | 180 md | 179.8 md | 3% | Aprobado |
| Presión extrapolada p* | 3,650 psia | 3,650.0 psia | 2 psi | Aprobado |
| Daño (skin), Horner | +4.5 | +4.47 | ±0.3 | Aprobado |
| Permeabilidad y daño, derivada log-log | 180 md, +4.5 | Dentro de la tolerancia | 3%, ±0.3 | Aprobado |
| Deconvolución de una prueba multicaudal, von Schroeter, Hollaender y Gringarten (2004) | Respuesta unitaria conocida | Recuperados | 5 %, o 10 % con ruido de sensor de 0,05 psi | Aprobado |
La pendiente de referencia es m = 162,6 qBμ / kh. La biblioteca de modelos también se verifica frente a los límites analíticos publicados: Theis (1935), van Everdingen y Hurst (1949), Gringarten, Ramey y Raghavan (1974), Cinco-Ley y Samaniego (1981) y Warren y Root (1963). El ejemplo 5.1 de Earlougher (1977) da k = 12,8 md y un daño de 8,6; Converge reproduce ambos dentro del 5 %.
DECLINACIÓN
Declinación de Arps
Arps, J. J. (1945), «Analysis of Decline Curves», Transactions of the AIME 160. Un historial sintético de 36 meses a partir de qi = 3.500 STB/d, Di = 60 %/año, b = 0,70.
| Resultado | Respuesta conocida | Converge | Tolerancia | Estado |
|---|---|---|---|---|
| qi ajustado, datos exactos | 3,500 STB/d | 3,500 STB/d | 1 STB/d | Aprobado |
| Di y b ajustados, datos exactos | 60%/yr, 0.70 | 60.0%/yr, 0.700 | 3 decimales | Aprobado |
| qi, Di y b ajustados con 3 % de ruido | Los valores generadores | Recuperados | 5 %, 8 puntos, 0,15 | Aprobado |
| EUR con un límite de 30 años | 4.78 MMSTB | 4.78 MMSTB | 0.005 MMSTB | Aprobado |
| Caudal y acumulada, Arps b de 0 a 1,5, hiperbólico modificado, Duong, exponencial estirada | Biblioteca independiente de código abierto de análisis de declinación | Coincide | 10⁻⁶ relativo | Aprobado |
BALANCE DE MATERIALES
Balance de materiales en un yacimiento con N conocido
Havlena, D. y Odeh, A. S. (1963), «The Material Balance as an Equation of a Straight Line», Journal of Petroleum Technology 15(8). Un yacimiento de petróleo volumétrico cuyo historial de producción y presión se construyó a partir de N = 45 MMSTB.
| Caso | Referencia | Converge | Tolerancia | Estado |
|---|---|---|---|---|
| Yacimiento volumétrico, sin acuífero | N = 45.00 MMSTB | 45.03 MMSTB | 1% | Aprobado |
| Mismo yacimiento, acuífero erróneo activado | No se muestra ninguna respuesta errónea | Ajuste bloqueado por no ser fiable | — | Corregido |
| Entrada de agua de un acuífero radial infinito | Ahmed (2006), ejemplo 10-4 | Constante y entrada de agua a 6 y 12 meses | 0.5% | Aprobado |
| Extracción subterránea F, cada fila | Ahmed (2006), ejemplo 11-1 | Reproducido | 0.01 Mbbl | Aprobado |
| Balance de materiales dinámico de petróleo | Solución analítica de círculo cerrado | N y b_pss | 2% | Aprobado |
| Balance de materiales dinámico de gas | Yacimiento de gas en régimen pseudoestacionario, G conocido | 0.04% | 1% | Aprobado |
ECONOMÍA
VPN y TIR
El flujo de caja descontado se verifica frente a sumas a mano y a una biblioteca financiera independiente de código abierto. El caso más sencillo para comprobarlo usted mismo: −100 y luego 30 al año durante cinco años, descuento al final del periodo al 10 %.
| Resultado | Cálculo manual | Converge | Tolerancia | Estado |
|---|---|---|---|---|
| VPN al 10 % | 12.476 | 12.476 | 10⁻³ | Aprobado |
| TIR | 15.24% | 15.24% | 10⁻⁹ | Aprobado |
| Filas de regalía e impuesto, caso de 10 años | Bucle independiente, cada año | Coincide | 10⁻⁶ | Aprobado |
| Aritmética de reparto de producción | Hoja manual: tope de cost oil, tramo del factor R, impuesto, participación del Estado | Coincide | 10⁻⁹ | Aprobado |
Las plantillas fiscales de la PIA 2021 de Nigeria son valores de evaluación preliminar codificados a partir de los anexos de la ley. No son asesoramiento fiscal y no forman parte de este benchmark.
POZOS
Levantamiento vertical y aporte del yacimiento
Las correlaciones de flujo multifásico se verifican frente a tres implementaciones independientes de código abierto. No se encontró ningún ejemplo resuelto de libro de texto con respuesta publicada para estas correlaciones, así que la concordancia es con código independiente.
| Elemento | Referencia | Concordancia |
|---|---|---|
| Beggs y Brill (1973) | Dos implementaciones independientes, 14 estados de flujo | 0.1% |
| Gray (1974) | Dos implementaciones independientes, 14 estados de flujo | 0.2% |
| Duns y Ros (1963) | Implementación independiente, flujo burbuja, tapón y niebla | 1% |
| Perfil completo de un pozo de gas, 10.000 ft, seco y húmedo | Implementación independiente | 0.02–0.18% |
| Perfil monofásico de líquido, pozo desviado | Hidrostática exacta más fricción de Colebrook | 10⁻⁹ |
LÍMITES
Lo que aún no está validado
Estos puntos siguen abiertos. Ninguno se presenta como validado en ninguna parte del producto.
- La envolvente de fases EOS completa frente a un paquete PVT comercial independiente. Los puntos de saturación y los flashes están verificados; la envolvente completa no.
- La red de superficie frente a herramientas independientes de modelado de redes. El balance de masa y la contabilidad del gas cierran, pero no se ha hecho ninguna comparación externa.
- El pronóstico integrado frente a herramientas independientes de modelado integrado de activos. Coincide con el propio motor de balance de materiales de Converge dentro de 0,1 psi.
- SPE9 es solo una concordancia parcial: el caudal de petróleo sigue a la referencia, pero la producción de agua no. SPE3 no se ha intentado. Ninguno de los dos se reivindica.
- Curvas tipo de transitorio de caudal para gas, modelos de parafinas y asfaltenos, y varios ejemplos impresos de pruebas de pozo de libros de texto.
REPRODUCIR
Ejecute los casos usted mismo
Cada caso anterior puede volver a ejecutarse en su propia cuenta. Abra Integraciones y use el editor de scripts, o llame a las mismas funciones mediante la API REST o MCP. El fluido del Sample field es el del caso V1, así que su página de PVT muestra de inmediato el punto de burbuja de 2.748 psia.
SPE1
const r = converge.sim.run({ preset: 'spe1' });
table(r.field, 'SPE1 field results');PVT y factor Z
converge.pvt.blackOilTable({
fluid: { kind: 'black-oil', api: 35, gasSG: 0.75,
gor: 650, reservoirT: 200 } });
converge.pvt.gasProperties({
psia: 3000, TF: 200, gasSG: 0.7 });Pruebas de pozo
const s = converge.welltest.synthetic({
reservoir: { k: 180, h: 32.808, phi: 0.2, ct: 1e-5, mu: 1,
B: 1.2, rw: 0.3, skin: 4.5, C: 0.001 },
pi: 3650, q: 800, tpHr: 48, shutInHr: 72,
pointsPerPeriod: 70, boundary: { kind: 'infinite' } });
const b = s.buildup;
converge.welltest.horner({ dtHr: b.dtHr, pws: b.pws,
tpHr: b.tpHr, pwf: b.pwf, window: [5, 70],
fluid: { q: 800, B: 1.2, mu: 1, h: 32.808,
phi: 0.2, ct: 1e-5, rw: 0.3 } });Declinación, balance de materiales, VPN
// points: [{ t: days, q: STB/d }, ...]
converge.dca.fit({ kind: 'arps', points });
converge.mbal.analyse(converge.example('mbal.analyse'));
converge.econ.metrics({
cashFlows: [-100, 30, 30, 30, 30, 30],
discountRates: [0.1], discounting: 'end-period' });Las entradas están en unidades de campo por defecto. Llame a converge.setUnits('metric') para trabajar en bara, °C y Sm³.